Standing shocks in a two-fluid solar wind

NASA Technical Reports Server (NTRS)

Habbal, Shadia R.; Hu, You Qiu; Esser, Ruth

1994-01-01

We present a numerical study of the formation of standing shocks in the solar wind using a two-fluid time-dependent model in the presence of Alfven waves. Included in this model is the adiabatic cooling and thermal conduction of both electrons and protons. In this study, standing shocks develop in the flow when additional critical points form as a result of either localized momentum addition or rapid expansion of the flow tube below the existing sonic point. While the flow speed and density exhibit the same characteristics as found in earlier studies of the formation of standing shocks, the inclusion of electron and proton heat conduction produces different signatures in the electron and proton temperature profiles across the shock layer. Owing to the strong heat conduction, the electron temperature is nearly continuous across the shock, but its gradient has a negative jump across it, thus producing a net heat flux out of the shock layer. The proton temperature exhibits the same characteristics for shocks produced by momentum addition but behaves differently when the shock is formed by the rapid divergence of the flow tube. The adiabatic cooling in a rapidly diverging flow tube reduces the proton temperature so substantially that the proton heat conduction becomes negligible in the vicinity of the shock. As a result, protons experience a positive jump in temperature across the shock. While Alfven waves do not affect the formation of standing shocks, they contribute to the change of the mmomentum and energy balance across them. We also find that for this solar wind model the inclusion of thermal conduction and adiabatic cooling for the elctrons and protons increases significantly the range of parameters characterizing the formation of standing shocks over those previously found for isothermal and polytropic models.

Equations of state and melting curve of boron carbide in the high-pressure range of shock compression

DOE Office of Scientific and Technical Information (OSTI.GOV)

Molodets, A. M., E-mail: molodets@icp.ac.ru; Golyshev, A. A.; Shakhrai, D. V.

We have constructed the equations of state for crystalline boron carbide B{sub 11}C (C–B–C) and its melt under high dynamic and static pressures. A kink on the shock adiabat for boron carbide has been revealed in the pressure range near 100 GPa, and the melting curve with negative curvature in the pressure range 0–120 GPa has been calculated. The results have been used for interpreting the kinks on the shock adiabat for boron carbide in the pressure range of 0–400 GPa.

Shock timing experiments on the National Ignition Facility: Initial results and comparison with simulation

NASA Astrophysics Data System (ADS)

Robey, H. F.; Boehly, T. R.; Celliers, P. M.; Eggert, J. H.; Hicks, D.; Smith, R. F.; Collins, R.; Bowers, M. W.; Krauter, K. G.; Datte, P. S.; Munro, D. H.; Milovich, J. L.; Jones, O. S.; Michel, P. A.; Thomas, C. A.; Olson, R. E.; Pollaine, S.; Town, R. P. J.; Haan, S.; Callahan, D.; Clark, D.; Edwards, J.; Kline, J. L.; Dixit, S.; Schneider, M. B.; Dewald, E. L.; Widmann, K.; Moody, J. D.; Döppner, T.; Radousky, H. B.; Throop, A.; Kalantar, D.; DiNicola, P.; Nikroo, A.; Kroll, J. J.; Hamza, A. V.; Horner, J. B.; Bhandarkar, S. D.; Dzenitis, E.; Alger, E.; Giraldez, E.; Castro, C.; Moreno, K.; Haynam, C.; LaFortune, K. N.; Widmayer, C.; Shaw, M.; Jancaitis, K.; Parham, T.; Holunga, D. M.; Walters, C. F.; Haid, B.; Mapoles, E. R.; Sater, J.; Gibson, C. R.; Malsbury, T.; Fair, J.; Trummer, D.; Coffee, K. R.; Burr, B.; Berzins, L. V.; Choate, C.; Brereton, S. J.; Azevedo, S.; Chandrasekaran, H.; Eder, D. C.; Masters, N. D.; Fisher, A. C.; Sterne, P. A.; Young, B. K.; Landen, O. L.; Van Wonterghem, B. M.; MacGowan, B. J.; Atherton, J.; Lindl, J. D.; Meyerhofer, D. D.; Moses, E.

2012-04-01

Capsule implosions on the National Ignition Facility (NIF) [Lindl et al., Phys. Plasmas 11, 339 (2004)] are underway with the goal of compressing deuterium-tritium (DT) fuel to a sufficiently high areal density (ρR) to sustain a self-propagating burn wave required for fusion power gain greater than unity. These implosions are driven with a carefully tailored sequence of four shock waves that must be timed to very high precision in order to keep the DT fuel on a low adiabat. Initial experiments to measure the strength and relative timing of these shocks have been conducted on NIF in a specially designed surrogate target platform known as the keyhole target. This target geometry and the associated diagnostics are described in detail. The initial data are presented and compared with numerical simulations. As the primary goal of these experiments is to assess and minimize the adiabat in related DT implosions, a methodology is described for quantifying the adiabat from the shock velocity measurements. Results are contrasted between early experiments that exhibited very poor shock timing and subsequent experiments where a modified target geometry demonstrated significant improvement.

PLANETARY EMBRYO BOW SHOCKS AS A MECHANISM FOR CHONDRULE FORMATION

DOE Office of Scientific and Technical Information (OSTI.GOV)

Mann, Christopher R.; Boley, Aaron C.; Morris, Melissa A.

2016-02-20

We use radiation hydrodynamics with direct particle integration to explore the feasibility of chondrule formation in planetary embryo bow shocks. The calculations presented here are used to explore the consequences of a Mars-size planetary embryo traveling on a moderately excited orbit through the dusty, early environment of the solar system. The embryo’s eccentric orbit produces a range of supersonic relative velocities between the embryo and the circularly orbiting gas and dust, prompting the formation of bow shocks. Temporary atmospheres around these embryos, which can be created via volatile outgassing and gas capture from the surrounding nebula, can non-trivially affect thermalmore » profiles of solids entering the shock. We explore the thermal environment of solids that traverse the bow shock at different impact radii, the effects that planetoid atmospheres have on shock morphologies, and the stripping efficiency of planetoidal atmospheres in the presence of high relative winds. Simulations are run using adiabatic and radiative conditions, with multiple treatments for the local opacities. Shock speeds of 5, 6, and 7 km s{sup −1} are explored. We find that a high-mass atmosphere and inefficient radiative conditions can produce peak temperatures and cooling rates that are consistent with the constraints set by chondrule furnace studies. For most conditions, the derived cooling rates are potentially too high to be consistent with chondrule formation.« less

Propagation of exponential shock wave in an axisymmetric rotating non-ideal dusty gas

NASA Astrophysics Data System (ADS)

Nath, G.

2016-09-01

One-dimensional unsteady isothermal and adiabatic flow behind a strong exponential shock wave propagating in a rotational axisymmetric mixture of non-ideal gas and small solid particles, which has variable azimuthal and axial fluid velocities, is analyzed. The shock wave is driven out by a piston moving with time according to exponential law. The azimuthal and axial components of the fluid velocity in the ambient medium are assumed to be varying and obeying exponential laws. In the present work, small solid particles are considered as pseudo-fluid with the assumption that the equilibrium flow-conditions are maintained in the flow-field, and the viscous-stress and heat conduction of the mixture are negligible. Solutions are obtained in both the cases, when the flow between the shock and the piston is isothermal or adiabatic by taking into account the components of vorticity vector and compressibility. It is found that the assumption of zero temperature gradient brings a profound change in the density, axial component of vorticity vector and compressibility distributions as compared to that of the adiabatic case. To investigate the behavior of the flow variables and the influence on the shock wave propagation by the parameter of non-idealness of the gas overline{b} in the mixture as well as by the mass concentration of solid particles in the mixture Kp and by the ratio of the density of solid particles to the initial density of the gas G1 are worked out in detail. It is interesting to note that the shock strength increases with an increase in G1 ; whereas it decreases with an increase in overline{b} . Also, a comparison between the solutions in the cases of isothermal and adiabatic flows is made.

Polarized Balmer line emission from supernova remnant shock waves efficiently accelerating cosmic rays

NASA Astrophysics Data System (ADS)

Shimoda, Jiro; Ohira, Yutaka; Yamazaki, Ryo; Laming, J. Martin; Katsuda, Satoru

2018-01-01

Linearly polarized Balmer line emissions from supernova remnant shocks are studied taking into account the energy loss of the shock owing to the production of non-thermal particles. The polarization degree depends on the downstream temperature and the velocity difference between upstream and downstream regions. The former is derived once the line width of the broad component of the H α emission is observed. Then, the observation of the polarization degree tells us the latter. At the same time, the estimated value of the velocity difference independently predicts adiabatic downstream temperature that is derived from Rankine Hugoniot relations for adiabatic shocks. If the actually observed downstream temperature is lower than the adiabatic temperature, there is a missing thermal energy which is consumed for particle acceleration. It is shown that a larger energy-loss rate leads to more highly polarized H α emission. Furthermore, we find that polarized intensity ratio of H β to H α also depends on the energy-loss rate and that it is independent of uncertain quantities such as electron temperature, the effect of Lyman line trapping and our line of sight.

The adiabatic energy change of plasma electrons and the frame dependence of the cross-shock potential at collisionless magnetosonic shock waves

NASA Technical Reports Server (NTRS)

Goodrich, C. C.; Scudder, J. D.

1984-01-01

The adiabatic energy gain of electrons in the stationary electric and magnetic field structure of collisionless shock waves was examined analytically in reference to conditions of the earth's bow shock. The study was performed to characterize the behavior of electrons interacting with the cross-shock potential. A normal incidence frame (NIF) was adopted in order to calculate the reversible energy change across a time stationary shock, and comparisons were made with predictions made by the de Hoffman-Teller (HT) model (1950). The electron energy gain, about 20-50 eV, is demonstrated to be consistent with a 200-500 eV potential jump in the bow shock quasi-perpendicular geometry. The electrons lose energy working against the solar wind motional electric field. The reversible energy process is close to that modeled by HT, which predicts that the motional electric field vanishes and the electron energy gain from the electric potential is equated to the ion energy loss to the potential.

Bow Shock in Merging Cluster A520: The Edge of the Radio Halo and the Electron–Proton Equilibration Timescale

NASA Astrophysics Data System (ADS)

Wang, Qian H. S.; Giacintucci, Simona; Markevitch, Maxim

2018-04-01

We studied the prominent bow shock in the merging galaxy cluster A520 using a deep Chandra X-ray observation and archival VLA radio data. This shock is a useful diagnostic tool, owing to its clear geometry and relatively high Mach number. At the “nose” of the shock, we measure a Mach number of M={2.4}-0.2+0.4. The shock becomes oblique away from the merger axis, with the Mach number falling to ≃1.6 around 30° from the nose. The electron temperature immediately behind the shock nose is consistent with that from the Rankine–Hugoniot adiabat, and is higher (at a 95% confidence) than expected for adiabatic compression of electrons followed by Coulomb electron–proton equilibration, indicating the presence of equilibration mechanisms faster than Coulomb collisions. This is similar to an earlier finding for the Bullet cluster. We also combined four archival VLA data sets to obtain a better image of the cluster’s giant radio halo at 1.4 GHz. An abrupt edge of the radio halo traces the shock front, and no emission is detected in the pre-shock region. If the radio edge were due only to adiabatic compression of relativistic electrons in pre-shock plasma, we would expect a pre-shock radio emission detectable in this radio data set; however, an interferometric artifact dominates the uncertainty, so we cannot rule this model out. Other interesting features of the radio halo include a peak at the remnant of the cool core, suggesting that the core used to have a radio minihalo, and a peak marking a possible region of high turbulence.

Solar-flare-induced Forbush decreases - Dependence on shock wave geometry

NASA Technical Reports Server (NTRS)

Thomas, B. T.; Gall, R.

1984-01-01

It is argued that the principal mechanism for the association of Forbush decreases with the passage of a solar flare shock wave is prolonged containment of cosmic ray particles behind the flare compression region, which acts as a semipermeable obstacle to particle motion along the field lines, leading to additional adiabatic cooling of the particles. Liouville's theorem is used to calculate the instantaneous distribution function at 1 AU for each particle arriving at the earth. By averaging over a large number of individual estimates, a representative estimate of the omnidirectional phase space density and the corresponding particle intensity is obtained. The energy change of individual particles at the shocks is found to be small in comparison to the energy lost by adiabatic cooling of the cosmic rays between the shock wave and the sun. The effects of particle rigidity, diffusion coefficient, and flare longitude on the magnitude of the Forbush decrease are quantitatively investigated.

Modeling of Particle Acceleration at Multiple Shocks via Diffusive Shock Acceleration: Preliminary Results

NASA Technical Reports Server (NTRS)

Parker, L. Neergaard; Zank, G. P.

2013-01-01

Successful forecasting of energetic particle events in space weather models require algorithms for correctly predicting the spectrum of ions accelerated from a background population of charged particles. We present preliminary results from a model that diffusively accelerates particles at multiple shocks. Our basic approach is related to box models in which a distribution of particles is diffusively accelerated inside the box while simultaneously experiencing decompression through adiabatic expansion and losses from the convection and diffusion of particles outside the box. We adiabatically decompress the accelerated particle distribution between each shock by either the method explored in Melrose and Pope (1993) and Pope and Melrose (1994) or by the approach set forth in Zank et al. (2000) where we solve the transport equation by a method analogous to operator splitting. The second method incorporates the additional loss terms of convection and diffusion and allows for the use of a variable time between shocks. We use a maximum injection energy (E(sub max)) appropriate for quasi-parallel and quasi-perpendicular shocks and provide a preliminary application of the diffusive acceleration of particles by multiple shocks with frequencies appropriate for solar maximum (i.e., a non-Markovian process).

Modeling of Particle Acceleration at Multiple Shocks Via Diffusive Shock Acceleration: Preliminary Results

NASA Technical Reports Server (NTRS)

Parker, Linda Neergaard; Zank, Gary P.

2013-01-01

We present preliminary results from a model that diffusively accelerates particles at multiple shocks. Our basic approach is related to box models (Protheroe and Stanev, 1998; Moraal and Axford, 1983; Ball and Kirk, 1992; Drury et al., 1999) in which a distribution of particles is diffusively accelerated inside the box while simultaneously experiencing decompression through adiabatic expansion and losses from the convection and diffusion of particles outside the box (Melrose and Pope, 1993; Zank et al., 2000). We adiabatically decompress the accelerated particle distribution between each shock by either the method explored in Melrose and Pope (1993) and Pope and Melrose (1994) or by the approach set forth in Zank et al. (2000) where we solve the transport equation by a method analogous to operator splitting. The second method incorporates the additional loss terms of convection and diffusion and allows for the use of a variable time between shocks. We use a maximum injection energy (Emax) appropriate for quasi-parallel and quasi-perpendicular shocks (Zank et al., 2000, 2006; Dosch and Shalchi, 2010) and provide a preliminary application of the diffusive acceleration of particles by multiple shocks with frequencies appropriate for solar maximum (i.e., a non-Markovian process).

Modeling of High-Velocity Flows in ITAM Impulse Facilities

DTIC Science & Technology

2010-04-01

up to 150 ms; Adiabatic compression wind tunnels up to 100 ms; Shock tubes... shock tubes. Basic and applied aerodynamic research has been performed in these wind tunnels in the range of Mach numbers М = 6 20 for many years...passage of a shock wave propagating over a cold rarefied gas filling the wind tunnel . When the gas heated in the shock wave (plug) passes around the

Shock formation in Ne, Ar, Kr, and Xe on deuterium gas puff implosions

DOE Office of Scientific and Technical Information (OSTI.GOV)

Narkis, J.; Rahman, H. U.; Ney, P.

2016-12-29

1- and 2-D simulations of 1-cm radius, gas-puff liners of Ne, Ar, Kr, and Xe imploding onto a deuterium target are conducted using the discharge parameters for the Zebra (1 MA, 130 ns) driver using the resistive MHD code MACH2. This is an implementation of the Staged Z-pinch concept, in which the target is driven to high-energy-density first by shock compression launched by a diffused azimuthal magnetic field (J×B force), and then by the adiabatic compression as the liner converges on axis. During the run-in phase, the initial shock heating preheats the deuterium plasma, with a subsequent stable, adiabatic compressionmore » heating the target to high energy density. Shock compression of the target coincides with the development of a J×B force at the target/liner interface. Stronger B-field transport and earlier shock compression increases with higher-Z liners, which results in an earlier shock arrival on axis. As a result, delayed shock formation in lower-Z liners yields a relative increase in shock heating, however, the 2-D simulations show an increased target isolation from magneto-Rayleigh-Taylor instability penetration, suggesting that an optimal balance between these two effects is reached in an Ar or Kr liner, rather than with Xe.« less

Stabilization of high-compression, indirect-drive inertial confinement fusion implosions using a 4-shock adiabat-shaped drive

DOE Office of Scientific and Technical Information (OSTI.GOV)

MacPhee, A. G.; Peterson, J. L.; Casey, D. T.

Hydrodynamic instabilities and poor fuel compression are major factors for capsule performance degradation in ignition experiments on the National Ignition Facility. Using a recently developed laser drive profile with a decaying first shock to tune the ablative Richtmyer-Meshkov (ARM) instability and subsequent in-flight Rayleigh-Taylor growth, we have demonstrated reduced growth compared to the standard ignition pulse whilst maintaining conditions for a low fuel adiabat needed for increased compression. Using in-flight x-ray radiography of pre-machined modulations, the first growth measurements using this new ARM-tuned drive have demonstrated instability growth reduction of ∼4× compared to the original design at a convergence ratiomore » of ∼2. Corresponding simulations give a fuel adiabat of ∼1.6, similar to the original goal and consistent with ignition requirements.« less

Stabilization of high-compression, indirect-drive inertial confinement fusion implosions using a 4-shock adiabat-shaped drive

DOE Office of Scientific and Technical Information (OSTI.GOV)

MacPhee, A. G.; Peterson, J. L.; Casey, D. T.

Hydrodynamic instabilities and poor fuel compression are major factors for capsule performance degradation in ignition experiments on the National Ignition Facility. Using a recently developed laser drive profile with a decaying first shock to tune the ablative Richtmyer-Meshkov (ARM) instability and subsequent in-flight Rayleigh-Taylor growth, we have demonstrated reduced growth compared to the standard ignition pulse whilst maintaining conditions for a low fuel adiabat needed for increased compression. Here, using in-flight x-ray radiography of pre-machined modulations, the first growth measurements using this new ARM-tuned drive have demonstrated instability growth reduction of ~4× compared to the original design at a convergencemore » ratio of ~2. Corresponding simulations give a fuel adiabat of ~1.6, similar to the original goal and consistent with ignition requirements.« less

Stabilization of high-compression, indirect-drive inertial confinement fusion implosions using a 4-shock adiabat-shaped drive

DOE PAGES

MacPhee, A. G.; Peterson, J. L.; Casey, D. T.; ...

2015-08-01

Hydrodynamic instabilities and poor fuel compression are major factors for capsule performance degradation in ignition experiments on the National Ignition Facility. Using a recently developed laser drive profile with a decaying first shock to tune the ablative Richtmyer-Meshkov (ARM) instability and subsequent in-flight Rayleigh-Taylor growth, we have demonstrated reduced growth compared to the standard ignition pulse whilst maintaining conditions for a low fuel adiabat needed for increased compression. Here, using in-flight x-ray radiography of pre-machined modulations, the first growth measurements using this new ARM-tuned drive have demonstrated instability growth reduction of ~4× compared to the original design at a convergencemore » ratio of ~2. Corresponding simulations give a fuel adiabat of ~1.6, similar to the original goal and consistent with ignition requirements.« less

The Ignition Physics Campaign on NIF: Status and Progress

NASA Astrophysics Data System (ADS)

Edwards, M. J.; Ignition Team

2016-03-01

We have made significant progress in ICF implosion performance on NIF since the 2011 IFSA. Employing a 3-shock, high adiabat CH (“High-Foot”) design, total neutron yields have increased 10-fold to 6.3 x1015 (a yield of ∼ 17 kJ, which is greater than the energy invested in the DT fuel ∼ 12kJ). At that level, the yield from alpha self-heating is essentially equivalent to the compression yield, indicating that we are close to the alpha self-heating regime. Low adiabat, 4-shock High Density Carbon (HDC) capsules have been imploded in conventional gas-filled hohlraums, and employing a 6 ns, 2-shock pulse, HDC capsules were imploded in near-vacuum hohlraums with overall coupling ∼ 98%. Both the 4- and 2-shock HDC capsules had very low mix and high yield over simulated performance. Rugby holraums have demonstrated uniform x-ray drive with minimal Cross Beam Energy Transfer (CBET), and we have made good progress in measuring and modelling growth of ablation front hydro instabilities.

Electron distributions upstream of the Comet Halley bow shock - Evidence for adiabatic heating

NASA Technical Reports Server (NTRS)

Larson, D. E.; Anderson, K. A.; Lin, R. P.; Carlson, C. W.; Reme, H.; Glassmeier, K. H.; Neubauer, F. M.

1992-01-01

Three-dimensional plasma electron (22 eV to 30 keV) observations upstream of Comet Halley bow shock, obtained by the RPA-1 COPERNIC (Reme Plasma Analyzer - Complete Positive Ion, Electron and Ram Negative Ion Measurements near Comet Halley) experiment on the Giotto spacecraft are reported. Besides electron distributions typical of the undisturbed solar wind and backstreaming electrons observed when the magnetic field line intersects the cometary bow shock, a new type of distribution, characterized by enhanced low energy (less than 100 eV) flux which peaks at 90-deg pitch angles is found. These are most prominent when the spacecraft is on field lines which pass close to but are not connected to the bow shock. The 90-deg pitch angle electrons appear to have been adiabatically heated by the increase in the magnetic field strength resulting from the compression of the upstream solar wind plasma by the cometary mass loading. A model calculation of this effect which agrees qualitatively with the observed 90-deg flux enhancements is presented.

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

DOE Office of Scientific and Technical Information (OSTI.GOV)

Kartavykh, Y. Y.; Dröge, W.; Gedalin, M.

2016-03-20

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 formore » 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.« less

Staged Z-pinch for the production of high-flux neutrons and net energy

DOE Office of Scientific and Technical Information (OSTI.GOV)

Wessel, Frank J.; Rahman, Hafiz Ur; Rostoker, Norman

A fusible target is embedded in a high Z liner, ohmically heated and then shock wave heated by implosion of an enveloping high Z liner. The target is adiabatically heated by compression, fusibly ignited and charged-particle heated as it is being ignited. A shock front forms as the liner implodes which shock front detaches from the more slowly moving liner, collides with the outer surface of the target, accelerates inward, rapidly heating the target, adiabatically compressing the target and liner and amplifying the current to converge the liner mass toward a central axis thereby compressing the target to a fusionmore » condition when it begins to ignite and produce charged particles. The charged particles are trapped in a large magnetic field surrounding the target. The energy of the charged particles is deposited into the target to further heat the target to produce an energy gain.« less

Shock-wave flow regimes at entry into the diffuser of a hypersonic ramjet engine: Influence of physical properties of the gas medium

NASA Astrophysics Data System (ADS)

Tarnavskii, G. A.

2006-07-01

The physical aspects of the effective-adiabatic-exponent model making it possible to decompose the total problem on modeling of high-velocity gas flows into individual subproblems (“physicochemical processes” and “ aeromechanics”), which ensures the creation of a universal and efficient computer complex divided into a number of independent units, have been analyzed. Shock-wave structures appearing at entry into the duct of a hypersonic aircraft have been investigated based on this methodology, and the influence of the physical properties of the gas medium in a wide range of variations of the effective adiabatic exponent has been studied.

Modeling of Particle Acceleration at Multiple Shocks Via Diffusive Shock Acceleration: Preliminary Results

NASA Astrophysics Data System (ADS)

Parker, L. N.; Zank, G. P.

2013-12-01

Successful forecasting of energetic particle events in space weather models require algorithms for correctly predicting the spectrum of ions accelerated from a background population of charged particles. We present preliminary results from a model that diffusively accelerates particles at multiple shocks. Our basic approach is related to box models (Protheroe and Stanev, 1998; Moraal and Axford, 1983; Ball and Kirk, 1992; Drury et al., 1999) in which a distribution of particles is diffusively accelerated inside the box while simultaneously experiencing decompression through adiabatic expansion and losses from the convection and diffusion of particles outside the box (Melrose and Pope, 1993; Zank et al., 2000). We adiabatically decompress the accelerated particle distribution between each shock by either the method explored in Melrose and Pope (1993) and Pope and Melrose (1994) or by the approach set forth in Zank et al. (2000) where we solve the transport equation by a method analogous to operator splitting. The second method incorporates the additional loss terms of convection and diffusion and allows for the use of a variable time between shocks. We use a maximum injection energy (Emax) appropriate for quasi-parallel and quasi-perpendicular shocks (Zank et al., 2000, 2006; Dosch and Shalchi, 2010) and provide a preliminary application of the diffusive acceleration of particles by multiple shocks with frequencies appropriate for solar maximum (i.e., a non-Markovian process).

Simulating cosmic ray physics on a moving mesh

NASA Astrophysics Data System (ADS)

Pfrommer, C.; Pakmor, R.; Schaal, K.; Simpson, C. M.; Springel, V.

2017-03-01

We discuss new methods to integrate the cosmic ray (CR) evolution equations coupled to magnetohydrodynamics on an unstructured moving mesh, as realized in the massively parallel AREPO code for cosmological simulations. We account for diffusive shock acceleration of CRs at resolved shocks and at supernova remnants in the interstellar medium (ISM) and follow the advective CR transport within the magnetized plasma, as well as anisotropic diffusive transport of CRs along the local magnetic field. CR losses are included in terms of Coulomb and hadronic interactions with the thermal plasma. We demonstrate the accuracy of our formalism for CR acceleration at shocks through simulations of plane-parallel shock tubes that are compared to newly derived exact solutions of the Riemann shock-tube problem with CR acceleration. We find that the increased compressibility of the post-shock plasma due to the produced CRs decreases the shock speed. However, CR acceleration at spherically expanding blast waves does not significantly break the self-similarity of the Sedov-Taylor solution; the resulting modifications can be approximated by a suitably adjusted, but constant adiabatic index. In first applications of the new CR formalism to simulations of isolated galaxies and cosmic structure formation, we find that CRs add an important pressure component to the ISM that increases the vertical scaleheight of disc galaxies and thus reduces the star formation rate. Strong external structure formation shocks inject CRs into the gas, but the relative pressure of this component decreases towards halo centres as adiabatic compression favours the thermal over the CR pressure.

Shock Propagation In Crustal Rock

DTIC Science & Technology

1991-04-29

liquid produced above the melting point during unloading. Figure 5 displays calculations of release adiabats in the mixed phase regime based on the...muscovite [Bridgman, 1949]. The zero-pressure densities at points 1, 2 and 3 correspond to mixtures of orthoclase + A1203+H20, of wadeite + kyanite ...shocked tantalum and the high pressure melting point , in Shock Waves in Condensed Matter-1983, edited by J. R. Asay, R. A. Graham and G. K. Straub, pp. 91

Estimate of Shock-Hugoniot Adiabat of Liquids from Hydrodynamics

NASA Astrophysics Data System (ADS)

Bouton, E.; Vidal, P.

2007-12-01

Shock states are generally obtained from shock velocity (D) and material velocity (u) measurements. In this paper, we propose a hydrodynamical method for estimating the (D-u) relation of Nitromethane from easily measured properties of the initial state. The method is based upon the differentiation of the Rankine-Hugoniot jump relations with the initial temperature considered as a variable and under the constraint of a unique nondimensional shock-Hugoniot. We then obtain an ordinary differential equation for the shock velocity D in the variable u. Upon integration, this method predicts the shock Hugoniot of liquid Nitromethane with a 5% accuracy for initial temperatures ranging from 250 K to 360 K.

Simulation of the Flow Field Associated with a Rocket Thruster Having an Attached Panel

NASA Technical Reports Server (NTRS)

Davoudzadeh, Farhad; Liu, Nan-Suey

2003-01-01

Two-dimensional inviscid and viscous numerical simulations are performed to predict the flow field induced by a H2-O2 rocket thruster and to provide insight into the heat load on the articles placed in the hot gas exhaust of the thruster under a variety of operating conditions, using the National Combustion Code (NCC). The simulations have captured physical details of the flow field, such as the plume formation and expansion, formation of the shock waves and their effects on the temperature and pressure distributions on the walls of the apparatus and the flat panel. Comparison between the computed results for 2-D and adiabatic walls and the related experimental measurements for 3-D and cooled walls shows that the results of the simulations are consistent with those obtained from the related rig tests.

Quantum effects in the capture of charged particles by dipolar polarizable symmetric top molecules. I. General axially nonadiabatic channel treatment.

PubMed

Auzinsh, M; Dashevskaya, E I; Litvin, I; Nikitin, E E; Troe, J

2013-08-28

The rate coefficients for capture of charged particles by dipolar polarizable symmetric top molecules in the quantum collision regime are calculated within an axially nonadiabatic channel approach. It uses the adiabatic approximation with respect to rotational transitions of the target within first-order charge-dipole interaction and takes into account the gyroscopic effect that decouples the intrinsic angular momentum from the collision axis. The results are valid for a wide range of collision energies (from single-wave capture to the classical limit) and dipole moments (from the Vogt-Wannier and fly-wheel to the adiabatic channel limit).

Analysis of geometric phase effects in the quantum-classical Liouville formalism.

PubMed

Ryabinkin, Ilya G; Hsieh, Chang-Yu; Kapral, Raymond; Izmaylov, Artur F

2014-02-28

We analyze two approaches to the quantum-classical Liouville (QCL) formalism that differ in the order of two operations: Wigner transformation and projection onto adiabatic electronic states. The analysis is carried out on a two-dimensional linear vibronic model where geometric phase (GP) effects arising from a conical intersection profoundly affect nuclear dynamics. We find that the Wigner-then-Adiabatic (WA) QCL approach captures GP effects, whereas the Adiabatic-then-Wigner (AW) QCL approach does not. Moreover, the Wigner transform in AW-QCL leads to an ill-defined Fourier transform of double-valued functions. The double-valued character of these functions stems from the nontrivial GP of adiabatic electronic states in the presence of a conical intersection. In contrast, WA-QCL avoids this issue by starting with the Wigner transform of single-valued quantities of the full problem. As a consequence, GP effects in WA-QCL can be associated with a dynamical term in the corresponding equation of motion. Since the WA-QCL approach uses solely the adiabatic potentials and non-adiabatic derivative couplings as an input, our results indicate that WA-QCL can capture GP effects in two-state crossing problems using first-principles electronic structure calculations without prior diabatization or introduction of explicit phase factors.

Analysis of geometric phase effects in the quantum-classical Liouville formalism

DOE Office of Scientific and Technical Information (OSTI.GOV)

Ryabinkin, Ilya G.; Izmaylov, Artur F.; Chemical Physics Theory Group, Department of Chemistry, University of Toronto, Toronto, Ontario M5S 3H6

2014-02-28

We analyze two approaches to the quantum-classical Liouville (QCL) formalism that differ in the order of two operations: Wigner transformation and projection onto adiabatic electronic states. The analysis is carried out on a two-dimensional linear vibronic model where geometric phase (GP) effects arising from a conical intersection profoundly affect nuclear dynamics. We find that the Wigner-then-Adiabatic (WA) QCL approach captures GP effects, whereas the Adiabatic-then-Wigner (AW) QCL approach does not. Moreover, the Wigner transform in AW-QCL leads to an ill-defined Fourier transform of double-valued functions. The double-valued character of these functions stems from the nontrivial GP of adiabatic electronic statesmore » in the presence of a conical intersection. In contrast, WA-QCL avoids this issue by starting with the Wigner transform of single-valued quantities of the full problem. As a consequence, GP effects in WA-QCL can be associated with a dynamical term in the corresponding equation of motion. Since the WA-QCL approach uses solely the adiabatic potentials and non-adiabatic derivative couplings as an input, our results indicate that WA-QCL can capture GP effects in two-state crossing problems using first-principles electronic structure calculations without prior diabatization or introduction of explicit phase factors.« less

Adiabatic invariants in stellar dynamics, 3: Application to globular cluster evolution

NASA Technical Reports Server (NTRS)

Weinberg, Martin D.

1994-01-01

The previous two companion papers demonstrate that slowly varying perturbations may not result in adiabatic cutoffs and provide a formalism for computing the long-term effects of time-dependent perturbations on stellar systems. Here, the theory is implemented in a Fokker-Planck code and a suite of runs illustrating the effects of shock heating on globular cluster evolution are described. Shock heating alone results in considerable mass loss for clusters with R(sub g) less than or approximately 8 kpc: a concentration c = 1.5 cluster with R(sub g) kpc loses up to 95% of its initial mass in 15 Gyr. Only those with concentration c greater than or approximately 1.3 survive disk shocks inside of this radius. Other effects, such as mass loss by stellar evolution, will decrease this survival bound. Loss of the initial halo together with mass segregation leads to mass spectral indices, x, which may be considerably larger than their initial values.

Electron heating in quasi-perpendicular shocks - A Monte Carlo simulation

NASA Technical Reports Server (NTRS)

Veltri, Pierluigi; Mangeney, Andre; Scudder, Jack D.

1990-01-01

To study the problem of electron heating in quasi-perpendicular shocks, under the combined effects of 'reversible' motion, in the shock electric potential and magnetic field, and wave-particle interactions a diffusion equation was derived, in the drift (adiabatic) approximation and it was solved by using a Monte Carlo method. The results show that most of the observations can be explained within this framework. The simulation has also definitively shown that the electron parallel temperature is determined by the dc electromagnetic field and not by any wave particle induced heating. Wave-particle interactions are effective in smoothing out the large gradients in phase space produced by the 'reversible' motion of the electrons, thus producing a 'cooling' of the electrons. Some constraints on the wave-particle interaction process may be obtained from a detailed comparison between the simulation and observations. In particular, it appears that the adiabatic approximation must be violated in order to explain the observed evolution of the perpendicular temperature.

Precision shock tuning on the national ignition facility.

PubMed

Robey, H F; Celliers, P M; Kline, J L; Mackinnon, A J; Boehly, T R; Landen, O L; Eggert, J H; Hicks, D; Le Pape, S; Farley, D R; Bowers, M W; Krauter, K G; Munro, D H; Jones, O S; Milovich, J L; Clark, D; Spears, B K; Town, R P J; Haan, S W; Dixit, S; Schneider, M B; Dewald, E L; Widmann, K; Moody, J D; Döppner, T D; Radousky, H B; Nikroo, A; Kroll, J J; Hamza, A V; Horner, J B; Bhandarkar, S D; Dzenitis, E; Alger, E; Giraldez, E; Castro, C; Moreno, K; Haynam, C; LaFortune, K N; Widmayer, C; Shaw, M; Jancaitis, K; Parham, T; Holunga, D M; Walters, C F; Haid, B; Malsbury, T; Trummer, D; Coffee, K R; Burr, B; Berzins, L V; Choate, C; Brereton, S J; Azevedo, S; Chandrasekaran, H; Glenzer, S; Caggiano, J A; Knauer, J P; Frenje, J A; Casey, D T; Johnson, M Gatu; Séguin, F H; Young, B K; Edwards, M J; Van Wonterghem, B M; Kilkenny, J; MacGowan, B J; Atherton, J; Lindl, J D; Meyerhofer, D D; Moses, E

2012-05-25

Ignition implosions on the National Ignition Facility [J. D. Lindl et al., Phys. Plasmas 11, 339 (2004)] are underway with the goal of compressing deuterium-tritium fuel to a sufficiently high areal density (ρR) to sustain a self-propagating burn wave required for fusion power gain greater than unity. These implosions are driven with a very carefully tailored sequence of four shock waves that must be timed to very high precision to keep the fuel entropy and adiabat low and ρR high. The first series of precision tuning experiments on the National Ignition Facility, which use optical diagnostics to directly measure the strength and timing of all four shocks inside a hohlraum-driven, cryogenic liquid-deuterium-filled capsule interior have now been performed. The results of these experiments are presented demonstrating a significant decrease in adiabat over previously untuned implosions. The impact of the improved shock timing is confirmed in related deuterium-tritium layered capsule implosions, which show the highest fuel compression (ρR~1.0 g/cm(2)) measured to date, exceeding the previous record [V. Goncharov et al., Phys. Rev. Lett. 104, 165001 (2010)] by more than a factor of 3. The experiments also clearly reveal an issue with the 4th shock velocity, which is observed to be 20% slower than predictions from numerical simulation.

Estimate of shock-Hugoniot adiabat of liquids from hydrodyamics

DOE Office of Scientific and Technical Information (OSTI.GOV)

Bouton, E.; Vidal, P.

2007-12-12

Shock states are generally obtained from shock velocity (D) and material velocity (u) measurements. In this paper, we propose a hydrodynamical method for estimating the (D-u) relation of Nitromethane from easily measured properties of the initial state. The method is based upon the differentiation of the Rankine-Hugoniot jump relations with the initial temperature considered as a variable and under the constraint of a unique nondimensional shock-Hugoniot. We then obtain an ordinary differential equation for the shock velocity D in the variable u. Upon integration, this method predicts the shock Hugoniot of liquid Nitromethane with a 5% accuracy for initial temperaturesmore » ranging from 250 K to 360 K.« less

Electron Heating in Low Mach Number Perpendicular Shocks. II. Dependence on the Pre-shock Conditions

NASA Astrophysics Data System (ADS)

Guo, Xinyi; Sironi, Lorenzo; Narayan, Ramesh

2018-05-01

Recent X-ray observations of merger shocks in galaxy clusters have shown that the post-shock plasma is two-temperature, with the protons being hotter than the electrons. In this work, the second of a series, we investigate the efficiency of irreversible electron heating in perpendicular low Mach number shocks, by means of two-dimensional particle-in-cell simulations. We consider values of plasma beta (the ratio of thermal and magnetic pressures) in the range 4 ≲ β p0 ≲ 32, and sonic Mach number (the ratio of shock speed to pre-shock sound speed) in the range 2 ≲ M s ≲ 5, as appropriate for galaxy cluster shocks. As shown in Paper I, magnetic field amplification—induced by shock compression of the pre-shock field, or by strong proton cyclotron and mirror modes accompanying the relaxation of proton temperature anisotropy—can drive the electron temperature anisotropy beyond the threshold of the electron whistler instability. The growth of whistler waves breaks the electron adiabatic invariance, and allows for efficient entropy production. We find that the post-shock electron temperature T e2 exceeds the adiabatic expectation {T}e2,{ad} by an amount ({T}e2-{T}e2,{ad})/{T}e0≃ 0.044 {M}s({M}s-1) (here, T e0 is the pre-shock temperature), which depends only weakly on the plasma beta over the range 4 ≲ β p0 ≲ 32 that we have explored, as well as on the proton-to-electron mass ratio (the coefficient of ≃0.044 is measured for our fiducial {m}i/{m}e=49, and we estimate that it will decrease to ≃0.03 for the realistic mass ratio). Our results have important implications for current and future observations of galaxy cluster shocks in the radio band (synchrotron emission and Sunyaev–Zel’dovich effect) and at X-ray frequencies.

A velocity probe-based method for continuous detonation and shock measurement in near-field underwater explosion.

PubMed

Li, Kebin; Li, Xiaojie; Yan, Honghao; Wang, Xiaohong; Miao, Yusong

2017-12-01

A new velocity probe which permits recording the time history of detonation and shock waves has been developed by improving the commercial on principle and structure. A method based on the probe is then designed to measure the detonation velocity and near-field shock parameters in a single underwater explosion, by which the oblique shock wave front of cylindrical charges and the peak pressure attenuation curve of spherical explosive are obtained. A further derivation of detonation pressure, adiabatic exponent, and other shock parameters is conducted. The present method offers a novel and reliable parameter determination for near-field underwater explosion.

A velocity probe-based method for continuous detonation and shock measurement in near-field underwater explosion

NASA Astrophysics Data System (ADS)

Li, Kebin; Li, Xiaojie; Yan, Honghao; Wang, Xiaohong; Miao, Yusong

2017-12-01

A new velocity probe which permits recording the time history of detonation and shock waves has been developed by improving the commercial on principle and structure. A method based on the probe is then designed to measure the detonation velocity and near-field shock parameters in a single underwater explosion, by which the oblique shock wave front of cylindrical charges and the peak pressure attenuation curve of spherical explosive are obtained. A further derivation of detonation pressure, adiabatic exponent, and other shock parameters is conducted. The present method offers a novel and reliable parameter determination for near-field underwater explosion.

Cones in Supersonic Flow

NASA Technical Reports Server (NTRS)

Hantzsche, W.; Wendt, H.

1947-01-01

In the case of cones in axially symmetric flow of supersonic velocity, adiabatic compression takes place between shock wave and surface of the cone. Interpolation curves betwen shock polars and the surface are therefore necessary for the complete understanding of this type of flow. They are given in the present report by graphical-numerical integration of the differential equation for all cone angles and airspeeds.

Dissipative advective accretion disc solutions with variable adiabatic index around black holes

NASA Astrophysics Data System (ADS)

Kumar, Rajiv; Chattopadhyay, Indranil

2014-10-01

We investigated accretion on to black holes in presence of viscosity and cooling, by employing an equation of state with variable adiabatic index and multispecies fluid. We obtained the expression of generalized Bernoulli parameter which is a constant of motion for an accretion flow in presence of viscosity and cooling. We obtained all possible transonic solutions for a variety of boundary conditions, viscosity parameters and accretion rates. We identified the solutions with their positions in the parameter space of generalized Bernoulli parameter and the angular momentum on the horizon. We showed that a shocked solution is more luminous than a shock-free one. For particular energies and viscosity parameters, we obtained accretion disc luminosities in the range of 10- 4 - 1.2 times Eddington luminosity, and the radiative efficiency seemed to increase with the mass accretion rate too. We found steady state shock solutions even for high-viscosity parameters, high accretion rates and for wide range of composition of the flow, starting from purely electron-proton to lepton-dominated accretion flow. However, similar to earlier studies of inviscid flow, accretion shock was not obtained for electron-positron pair plasma.

Behavior of dusty real gas on adiabatic propagation of cylindrical imploding strong shock waves

NASA Astrophysics Data System (ADS)

Gangwar, P. K.

2018-05-01

In this paper, CCW method has been used to study the behavior of dusty real gas on adiabatic propagation of cylindrical imploding strong shock waves. The strength of overtaking waves is estimated under the assumption that both C+ and C- disturbances propagate in non-uniform region of same density distribution. It is assumed that the dusty gas is the mixture of a real gas and a large number of small spherical solid particles of uniform size. The solid particles are uniformly distributed in the medium. Maintaining equilibrium flow conditions, the expressions for shock strength has been derived both for freely propagation as well as under the effect of overtaking disturbances. The variation of all flow variables with propagation distance, mass concentration of solid particles in the mixture and the ratio of solid particles to the initial density of gas have been computed and discussed through graphs. It is found that the presence of dust particles in the gases medium has significant effects on the variation of flow variables and the shock is strengthened under the influence of overtaking disturbances. The results accomplished here been compared with those for ideal gas.

The effect of shock dynamics on compressibility of ignition-scale National Ignition Facility implosions

DOE Office of Scientific and Technical Information (OSTI.GOV)

Zylstra, A. B.; Frenje, J. A.; Séguin, F. H.

The effects of shock dynamics on compressibility of indirect-drive ignition-scale surrogate implosions, CH shells filled with D 3He gas, have been studied using charged-particle spectroscopy. Spectral measurements of D 3He protons produced at the shock-bang time probe the shock dynamics and in-flight characteristics of an implosion. The proton shock yield is found to vary by over an order of magnitude. A simple model relates the observed yield to incipient hot-spot adiabat, suggesting that implosions with rapid radiation-power increase during the main drive pulse may have a 2x higher hot-spot adiabat, potentially reducing compressibility. A self-consistent 1-D implosion model was usedmore » to infer the areal density (pR) and the shell center-of-mass radius (R cm) from the downshift of the shock-produced D 3He protons. The observed pR at shock-bang time is substantially higher for implosions, where the laser drive is on until near the compression bang time ('short-coast'), while longer-coasting implosions have lower pR. This corresponds to a much larger temporal difference between the shock- and compression-bang time in the long-coast implosions (~800 ps) than in the short-coast (~400 ps); this will be verified with a future direct bang-time diagnostic. This model-inferred differential bang time contradicts radiation-hydrodynamic simulations, which predict constant 700–800 ps differential independent of coasting time. This result is potentially explained by uncertainties in modeling late-time ablation drive on the capsule. In an ignition experiment, an earlier shock-bang time resulting in an earlier onset of shell deceleration, potentially reducing compression and, thus, fuel pR.« less

The effect of shock dynamics on compressibility of ignition-scale National Ignition Facility implosions

DOE PAGES

Zylstra, A. B.; Frenje, J. A.; Séguin, F. H.; ...

2014-11-03

The effects of shock dynamics on compressibility of indirect-drive ignition-scale surrogate implosions, CH shells filled with D 3He gas, have been studied using charged-particle spectroscopy. Spectral measurements of D 3He protons produced at the shock-bang time probe the shock dynamics and in-flight characteristics of an implosion. The proton shock yield is found to vary by over an order of magnitude. A simple model relates the observed yield to incipient hot-spot adiabat, suggesting that implosions with rapid radiation-power increase during the main drive pulse may have a 2x higher hot-spot adiabat, potentially reducing compressibility. A self-consistent 1-D implosion model was usedmore » to infer the areal density (pR) and the shell center-of-mass radius (R cm) from the downshift of the shock-produced D 3He protons. The observed pR at shock-bang time is substantially higher for implosions, where the laser drive is on until near the compression bang time ('short-coast'), while longer-coasting implosions have lower pR. This corresponds to a much larger temporal difference between the shock- and compression-bang time in the long-coast implosions (~800 ps) than in the short-coast (~400 ps); this will be verified with a future direct bang-time diagnostic. This model-inferred differential bang time contradicts radiation-hydrodynamic simulations, which predict constant 700–800 ps differential independent of coasting time. This result is potentially explained by uncertainties in modeling late-time ablation drive on the capsule. In an ignition experiment, an earlier shock-bang time resulting in an earlier onset of shell deceleration, potentially reducing compression and, thus, fuel pR.« less

On numerical model of one-dimensional time-dependent gas flows through bed of encapsulated phase change material

NASA Astrophysics Data System (ADS)

Lutsenko, N. A.; Fetsov, S. S.

2017-10-01

Mathematical model and numerical method are proposed for investigating the one-dimensional time-dependent gas flows through a packed bed of encapsulated Phase Change Material (PCM). The model is based on the assumption of interacting interpenetrating continua and includes equations of state, continuity, momentum conservation and energy for PCM and gas. The advantage of the method is that it does not require predicting the location of phase transition zone and can define it automatically as in a usual shock-capturing method. One of the applications of the developed numerical model is the simulation of novel Adiabatic Compressed Air Energy Storage system (A-CAES) with Thermal Energy Storage subsystem (TES) based on using the encapsulated PCM in packed bed. Preliminary test calculations give hope that the method can be effectively applied in the future for modelling the charge and discharge processes in such TES with PCM.

Herbig-Haro objects as the heads of radiative jets

NASA Technical Reports Server (NTRS)

Blondin, John M.; Konigl, Arieh; Fryxell, Bruce A.

1989-01-01

The interpretation of certain HH objects as the heads of nonadiabatic supersonic jets is examined using two-dimensional numerical simulations. It is found that radiative jets develop a dense shell between the jet shock and the leading bow shock when the cooling distance behind either one of these shocks is smaller than the jet radius. It is proposed that the radiatively cooling shell may account for the variable emission pattern from objects like HH 1. Also, it is suggested that HH objects with measured space velocities that exceed the spectroscopically inferred shock velocities could correspond to heavy jets in which the bow shock is effectively adiabatic. Low-excitation objects in which these velocities are comparable may represent light jets where the jet shock is nonradiative.

Shock-drift particle acceleration in superluminal shocks - A model for hot spots in extragalactic radio sources

NASA Technical Reports Server (NTRS)

Begelman, Mitchell C.; Kirk, John G.

1990-01-01

Shock-drift acceleration at relativistic shock fronts is investigated using a fully relativistic treatment of both the microphysics of the shock-drift acceleration and the macrophysics of the shock front. By explicitly tracing particle trajectories across shocks, it is shown how the adiabatic invariance of a particle's magnetic moment breaks down as the upstream shock speed becomes relativistic, and is recovered at subrelativistic velocities. These calculations enable the mean increase in energy of a particle which encounters the shock with a given pitch angle to be calculated. The results are used to construct the downstream electron distribution function in terms of the incident distribution function and the bulk properties of the shock. The synchrotron emissivity of the transmitted distribution is calculated, and it is demonstrated that amplification factors are easily obtained which are more than adequate to explain the observed constrasts in surface brightness between jets and hot spots.

Shock-Wave Pulse Compression and Stretching of Dodecane and Mineral Oils

NASA Astrophysics Data System (ADS)

Bannikova, I. A.; Zubareva, A. N.; Utkin, A. V.

2018-04-01

The behavior of dodecane, vacuum, and transformer oils under shock-wave pulse compression and stretching are studied experimentally. The wave profiles are registered using a VISAR laser interferometer. The shock adiabats, the dependence of the sound velocity on the pressure, and the maximum negative pressures developed in the studied liquids are determined. It is shown that the negative pressure value does not depend on the deformation rate in the case of oils and is a strong function of the compression pulse amplitude in the case of dodecane.

Another self-similar blast wave: Early time asymptote with shock heated electrons and high thermal conductivity

NASA Technical Reports Server (NTRS)

Cox, D. P.; Edgar, R. J.

1982-01-01

Accurate approximations are presented for the self-similar structures of nonradiating blast waves with adiabatic ions, isothermal electrons, and equation ion and electron temperatures at the shock. The cases considered evolve in cavities with power law ambient densities (including the uniform density case) and have negligible external pressure. The results provide the early time asymptote for systems with shock heating of electrons and strong thermal conduction. In addition, they provide analytical results against which two fluid numerical hydrodynamic codes can be checked.

An Integrated Development Environment for Adiabatic Quantum Programming

DOE Office of Scientific and Technical Information (OSTI.GOV)

Humble, Travis S; McCaskey, Alex; Bennink, Ryan S

2014-01-01

Adiabatic quantum computing is a promising route to the computational power afforded by quantum information processing. The recent availability of adiabatic hardware raises the question of how well quantum programs perform. Benchmarking behavior is challenging since the multiple steps to synthesize an adiabatic quantum program are highly tunable. We present an adiabatic quantum programming environment called JADE that provides control over all the steps taken during program development. JADE captures the workflow needed to rigorously benchmark performance while also allowing a variety of problem types, programming techniques, and processor configurations. We have also integrated JADE with a quantum simulation enginemore » that enables program profiling using numerical calculation. The computational engine supports plug-ins for simulation methodologies tailored to various metrics and computing resources. We present the design, integration, and deployment of JADE and discuss its use for benchmarking adiabatic quantum programs.« less

Gyroscopic effect in low-energy classical capture of a rotating quadrupolar diatom by an ion.

PubMed

Dashevskaya, Elena; Litvin, Iliya; Nikitin, Evgueni

2006-03-09

The low-energy capture of homonuclear diatoms by ions is due mainly to the long-range part of the interpartner potential with leading terms that correspond to charge-quadrupole interaction and charge-induced dipole interaction. The capture dynamics is described by the perturbed-rotor adiabatic potentials and the Coriolis interaction between manifold of states that belong to a given value of the intrinsic angular momentum. When the latter is large enough, it can noticeably affect the capture cross section calculated in the adiabatic channel approximation due to the gyroscopic property of a rotating diatom. This paper presents the low-energy (low-temperature) state-selected partial and mean capture cross sections (rate coefficients) for the charge-quadrupole interaction that include the gyroscopic effect (decoupling of intrinsic angular momentum from the collision axis), quantum correction for the diatom rotation, and the correction for the charge-induced dipole interaction. These results complement recent studies on the gyroscopic effect in the quantum regime of diatom-ion capture (Dashevskaya, E. I.; Litvin, I.; Nikitin, E. E.; Troe, J. J. Chem. Phys. 2004, 120, 9989-9997).

DOE Office of Scientific and Technical Information (OSTI.GOV)

Arévalo, P.; Churazov, E.; Zhuravleva, I.

X-ray images of galaxy clusters and gas-rich elliptical galaxies show a wealth of small-scale features that reflect fluctuations in density and/or temperature of the intracluster medium. In this paper we study these fluctuations in M87/Virgo to establish whether sound waves/shocks, bubbles, or uplifted cold gas dominate the structure. We exploit the strong dependence of the emissivity on density and temperature in different energy bands to distinguish between these processes. Using simulations we demonstrate that our analysis recovers the leading type of fluctuation even in the presence of projection effects and temperature gradients. We confirm the isobaric nature of cool filamentsmore » of gas entrained by buoyantly rising bubbles, extending to 7′ to the east and southwest, and the adiabatic nature of the weak shocks at 40″ and 3′ from the center. For features of ∼5–10 kpc, we show that the central 4′ × 4′ region is dominated by cool structures in pressure equilibrium with the ambient hotter gas while up to 30% of the variance in this region can be ascribed to adiabatic fluctuations. The remaining part of the central 14′ × 14′ region, excluding the arms and shocks described above, is dominated by apparently isothermal fluctuations (bubbles) with a possible admixture (at the level of ∼30%) of adiabatic (sound waves) and by isobaric structures. Larger features, of about 30 kpc, show a stronger contribution from isobaric fluctuations. The results broadly agree with a model based on feedback from an active galactic nucleus mediated by bubbles of relativistic plasma.« less

Impact induced depolarization of ferroelectric materials

NASA Astrophysics Data System (ADS)

Agrawal, Vinamra; Bhattacharya, Kaushik

2018-06-01

We study the large deformation dynamic behavior and the associated nonlinear electro-thermo-mechanical coupling exhibited by ferroelectric materials in adiabatic environments. This is motivated by a ferroelectric generator which involves pulsed power generation by loading the ferroelectric material with a shock, either by impact or a blast. Upon impact, a shock wave travels through the material inducing a ferroelectric to nonpolar phase transition giving rise to a large voltage difference in an open circuit situation or a large current in a closed circuit situation. In the first part of this paper, we provide a general continuum mechanical treatment of the situation assuming a sharp phase boundary that is possibly charged. We derive the governing laws, as well as the driving force acting on the phase boundary. In the second part, we use the derived equations and a particular constitutive relation that describes the ferroelectric to nonpolar phase transition to study a uniaxial plate impact problem. We develop a numerical method where the phase boundary is tracked but other discontinuities are captured using a finite volume method. We compare our results with experimental observations to find good agreement. Specifically, our model reproduces the observed exponential rise of charge as well as the resistance dependent Hugoniot. We conclude with a parameter study that provides detailed insight into various aspects of the problem.

The Curious Events Leading to the Theory of Shock Waves

NASA Technical Reports Server (NTRS)

Salas, Manuel D.

2006-01-01

We review the history of the development of the modern theory of shock waves. Several attempts at an early-theory quickly collapsed for lack of foundations in mathematics and thermodynamics. It is not until the works of Rankine and later Hugoniot that a full theory is established. Rankine is the first to show that within the shock a non-adiabatic process must occur. Hugoniot showed that in the absence of viscosity and heat conduction conservation of energy implies conservation of entropy in smooth regions and a jump in entropy across a shock. Even after the theory is fully developed, old notions continue to pervade the literature well into the early part of the 20th Century.

Implications of high-velocity interstellar H I absorption features

NASA Technical Reports Server (NTRS)

Cowie, L.; York, D. G.; Laurent, C.; Vidal-Madjar, A.

1979-01-01

Contributions to the interstellar H I column density at high velocities from immediate postshock gas and from the cooling gas behind a shock are compared. The detection of high-velocity H I in L-epsilon and L-delta for Iota Ori is reported and interpreted as cooling gas behind a shock of 100 km/s velocity. The immediate postshock gas should be observable for shock velocities greater than 200 km/s and permits direct determination of the velocities of adiabatic shocks in the interstellar medium. It is pointed out that interstellar L-alpha and L-beta lines may not have purely Lorentzian profiles if high-velocity H I is a widespread phenomenon.

Energetic storm particle events in the outer heliosphere

NASA Technical Reports Server (NTRS)

Mcdonald, F.; Trainor, J.; Mihalov, J.; Wolfe, J.; Webber, W.

1981-01-01

The evolution of energetic particle events with increasing heliocentric distance is studied through events of Pioneers 10 and 11. Beyond 12 AU the events become the dominant type of solar particle event at 1 AU, and the combined effects of adiabatic cooling and volume expansion rule out the possibility that the particles represent the confinement of the original particle population behind the shock. It is not established whether the particles originate from the solar wind by injection via post-shock enhancements or are energetic solar particles further energized by the shock, although their very long lifetime favors the solar wind origin.

Exact Time-Dependent Exchange-Correlation Potential in Electron Scattering Processes

NASA Astrophysics Data System (ADS)

Suzuki, Yasumitsu; Lacombe, Lionel; Watanabe, Kazuyuki; Maitra, Neepa T.

2017-12-01

We identify peak and valley structures in the exact exchange-correlation potential of time-dependent density functional theory that are crucial for time-resolved electron scattering in a model one-dimensional system. These structures are completely missed by adiabatic approximations that, consequently, significantly underestimate the scattering probability. A recently proposed nonadiabatic approximation is shown to correctly capture the approach of the electron to the target when the initial Kohn-Sham state is chosen judiciously, and it is more accurate than standard adiabatic functionals but ultimately fails to accurately capture reflection. These results may explain the underestimation of scattering probabilities in some recent studies on molecules and surfaces.

Adiabatic Compression Sensitivity of Liquid Fuels and Monopropellants

NASA Technical Reports Server (NTRS)

Ismail, Ismail M. K.; Hawkins, Tom W.

2000-01-01

Liquid rocket propellants can be sensitive to rapid compression. Such liquids may undergo decomposition and their handling may be accompanied with risk. Decomposition produces small gas bubbles in the liquid, which upon rapid compression may cause catastrophic explosions. The rapid compression can result from mechanical shocks applied on the tank containing the liquid or from rapid closure of the valves installed on the lines. It is desirable to determine the conditions that may promote explosive reactions. At Air Force Research Laboratory (AFRL), we constructed an apparatus and established a safe procedure for estimating the sensitivity of propellant materials towards mechanical shocks (Adiabatic Compression Tester). A sample is placed on a stainless steel U-tube, held isothermally at a temperature between 20 and 150 C then exposed to an abrupt mechanical shock of nitrogen gas at a pressure between 6.9 and 20.7 MPa (1000 to 3000 psi). The apparatus is computer interfaced and is driven with LABTECH NOTEBOOK-pro (registered) Software. In this presentation, the design of the apparatus is shown, the operating procedure is outlined, and the safety issues are addressed. The results obtained on different energetic materials are presented.

The MSW Effect and Matter Effects in Neutrino Oscillations

NASA Astrophysics Data System (ADS)

Smirnov, A. Yu

2005-01-01

The MSW (Mikheyev-Smirnov-Wolfenstein) effect is the adiabatic or partially adiabatic neutrino flavor conversion in media with varying density. The main notions related to the effect, its dynamics and physical picture are reviewed. The large mixing MSW effect is realized inside the Sun providing a solution of the solar neutrino problem. The small mixing MSW effect driven by the 1 3 mixing can be realized for the supernova (SN) neutrinos. Inside collapsing stars new elements of the MSW dynamics may show up: non-oscillatory transition, non-adiabatic conversion, time dependent adiabaticity violation induced by shock waves. Effects of the resonance enhancement and the parametric enhancement of oscillations can be realized for atmospheric and accelerator neutrinos in the Earth. Precise results for neutrino oscillations in low density media with arbitrary density profile are presented and the attenuation effect is described. The area of applications is the solar and SN neutrinos inside the Earth, and the results are crucial for the neutrino oscillation tomography.

The MSW Effect and Matter Effects in Neutrino Oscillations

NASA Astrophysics Data System (ADS)

Smirnov, A. Yu.

2006-03-01

The MSW (Mikheyev-Smirnov-Wolfenstein) effect is the adiabatic or partially adiabatic neutrino flavor conversion in media with varying density. The main notions related to the effect, its dynamics and physical picture are reviewed. The large mixing MSW effect is realized inside the Sun providing a solution of the solar neutrino problem. The small mixing MSW effect driven by the 1-3 mixing can be realized for the supernova (SN) neutrinos. Inside collapsing stars new elements of the MSW dynamics may show up: non-oscillatory transition, non-adiabatic conversion, time dependent adiabaticity violation induced by shock waves. Effects of the resonance enhancement and the parametric enhancement of oscillations can be realized for atmospheric and accelerator neutrinos in the Earth. Precise results for neutrino oscillations in low density media with arbitrary density profile are presented and the attenuation effect is described. The area of applications is the solar and SN neutrinos inside the Earth, and the results are crucial for the neutrino oscillation tomography.

Floating shock fitting via Lagrangian adaptive meshes

NASA Technical Reports Server (NTRS)

Vanrosendale, John

1995-01-01

In recent work we have formulated a new approach to compressible flow simulation, combining the advantages of shock-fitting and shock-capturing. Using a cell-centered on Roe scheme discretization on unstructured meshes, we warp the mesh while marching to steady state, so that mesh edges align with shocks and other discontinuities. This new algorithm, the Shock-fitting Lagrangian Adaptive Method (SLAM), is, in effect, a reliable shock-capturing algorithm which yields shock-fitted accuracy at convergence.

Magnetohydrodynamic modeling of three Van Allen Probes storms in 2012 and 2013

NASA Astrophysics Data System (ADS)

Paral, J.; Hudson, M. K.; Kress, B. T.; Wiltberger, M. J.; Wygant, J. R.; Singer, H. J.

2015-08-01

Coronal mass ejection (CME)-shock compression of the dayside magnetopause has been observed to cause both prompt enhancement of radiation belt electron flux due to inward radial transport of electrons conserving their first adiabatic invariant and prompt losses which at times entirely eliminate the outer zone. Recent numerical studies suggest that enhanced ultra-low frequency (ULF) wave activity is necessary to explain electron losses deeper inside the magnetosphere than magnetopause incursion following CME-shock arrival. A combination of radial transport and magnetopause shadowing can account for losses observed at radial distances into L = 4.5, well within the computed magnetopause location. We compare ULF wave power from the Electric Field and Waves (EFW) electric field instrument on the Van Allen Probes for the 8 October 2013 storm with ULF wave power simulated using the Lyon-Fedder-Mobarry (LFM) global magnetohydrodynamic (MHD) magnetospheric simulation code coupled to the Rice Convection Model (RCM). Two other storms with strong magnetopause compression, 8-9 October 2012 and 17-18 March 2013, are also examined. We show that the global MHD model captures the azimuthal magnetosonic impulse propagation speed and amplitude observed by the Van Allen Probes which is responsible for prompt acceleration at MeV energies reported for the 8 October 2013 storm. The simulation also captures the ULF wave power in the azimuthal component of the electric field, responsible for acceleration and radial transport of electrons, at frequencies comparable to the electron drift period. This electric field impulse has been shown to explain observations in related studies (Foster et al., 2015) of electron acceleration and drift phase bunching by the Energetic Particle, Composition, and Thermal Plasma Suite (ECT) instrument on the Van Allen Probes.

Evaluation of XHVRB for Capturing Explosive Shock Desensitization

NASA Astrophysics Data System (ADS)

Tuttle, Leah; Schmitt, Robert; Kittell, Dave; Harstad, Eric

2017-06-01

Explosive shock desensitization phenomena have been recognized for some time. It has been demonstrated that pressure-based reactive flow models do not adequately capture the basic nature of the explosive behavior. Historically, replacing the local pressure with a shock captured pressure has dramatically improved the numerical modeling approaches. Models based upon shock pressure or functions of entropy have recently been developed. A pseudo-entropy based formulation using the History Variable Reactive Burn model, as proposed by Starkenberg, was implemented into the Eulerian shock physics code CTH. Improvements in the shock capturing algorithm were made. The model is demonstrated to reproduce single shock behavior consistent with published pop plot data. It is also demonstrated to capture a desensitization effect based on available literature data, and to qualitatively capture dead zones from desensitization in 2D corner turning experiments. This models shows promise for use in modeling and simulation problems that are relevant to the desensitization phenomena. Issues are identified with the current implementation and future work is proposed for improving and expanding model capabilities. Sandia National Laboratories is a multi-mission laboratory managed and operated by Sandia Corporation, a wholly owned subsidiary of Lockheed Martin Corporation, for the U.S. DOE's National Nuclear Security Administration under contract DE-AC04-94AL85000.

Floating shock fitting via Lagrangian adaptive meshes

NASA Technical Reports Server (NTRS)

Vanrosendale, John

1994-01-01

In recent works we have formulated a new approach to compressible flow simulation, combining the advantages of shock-fitting and shock-capturing. Using a cell-centered Roe scheme discretization on unstructured meshes, we warp the mesh while marching to steady state, so that mesh edges align with shocks and other discontinuities. This new algorithm, the Shock-fitting Lagrangian Adaptive Method (SLAM) is, in effect, a reliable shock-capturing algorithm which yields shock-fitted accuracy at convergence. Shock-capturing algorithms like this, which warp the mesh to yield shock-fitted accuracy, are new and relatively untried. However, their potential is clear. In the context of sonic booms, accurate calculation of near-field sonic boom signatures is critical to the design of the High Speed Civil Transport (HSCT). SLAM should allow computation of accurate N-wave pressure signatures on comparatively coarse meshes, significantly enhancing our ability to design low-boom configurations for high-speed aircraft.

Converging shock flows for a Mie-Grüneisen equation of state

NASA Astrophysics Data System (ADS)

Ramsey, Scott D.; Schmidt, Emma M.; Boyd, Zachary M.; Lilieholm, Jennifer F.; Baty, Roy S.

2018-04-01

Previous work has shown that the one-dimensional (1D) inviscid compressible flow (Euler) equations admit a wide variety of scale-invariant solutions (including the famous Noh, Sedov, and Guderley shock solutions) when the included equation of state (EOS) closure model assumes a certain scale-invariant form. However, this scale-invariant EOS class does not include even simple models used for shock compression of crystalline solids, including many broadly applicable representations of Mie-Grüneisen EOS. Intuitively, this incompatibility naturally arises from the presence of multiple dimensional scales in the Mie-Grüneisen EOS, which are otherwise absent from scale-invariant models that feature only dimensionless parameters (such as the adiabatic index in the ideal gas EOS). The current work extends previous efforts intended to rectify this inconsistency, by using a scale-invariant EOS model to approximate a Mie-Grüneisen EOS form. To this end, the adiabatic bulk modulus for the Mie-Grüneisen EOS is constructed, and its key features are used to motivate the selection of a scale-invariant approximation form. The remaining surrogate model parameters are selected through enforcement of the Rankine-Hugoniot jump conditions for an infinitely strong shock in a Mie-Grüneisen material. Finally, the approximate EOS is used in conjunction with the 1D inviscid Euler equations to calculate a semi-analytical Guderley-like imploding shock solution in a metal sphere and to determine if and when the solution may be valid for the underlying Mie-Grüneisen EOS.

Monte Carlo shock-like solutions to the Boltzmann equation with collective scattering

NASA Technical Reports Server (NTRS)

Ellison, D. C.; Eichler, D.

1984-01-01

The results of Monte Carlo simulations of steady state shocks generated by a collision operator that isotropizes the particles by means of elastic scattering in some locally defined frame of reference are presented. The simulations include both the back reaction of accelerated particles on the inflowing plasma and the free escape of high-energy particles from finite shocks. Energetic particles are found to be naturally extracted out of the background plasma by the shock process with an efficiency in good quantitative agreement with an earlier analytic approximation (Eichler, 1983 and 1984) and observations (Gosling et al., 1981) of the entire particle spectrum at a quasi-parallel interplanetary shock. The analytic approximation, which allows a self-consistent determination of the effective adiabatic index of the shocked gas, is used to calculate the overall acceleration efficiency and particle spectrum for cases where ultrarelativistic energies are obtained. It is found that shocks of the strength necessary to produce galactic cosmic rays put approximately 15 percent of the shock energy into relativistic particles.

Onion-shell model of cosmic ray acceleration in supernova remnants

NASA Technical Reports Server (NTRS)

Bogdan, T. J.; Volk, H. J.

1983-01-01

A method is devised to approximate the spatially averaged momentum distribution function for the accelerated particles at the end of the active lifetime of a supernova remnant. The analysis is confined to the test particle approximation and adiabatic losses are oversimplified, but unsteady shock motion, evolving shock strength, and non-uniform gas flow effects on the accelerated particle spectrum are included. Monoenergetic protons are injected at the shock front. It is found that the dominant effect on the resultant accelerated particle spectrum is a changing spectral index with shock strength. High energy particles are produced in early phases, and the resultant distribution function is a slowly varying power law over several orders of magnitude, independent of the specific details of the supernova remnant.

Steady flow on to a conveyor belt - Causal viscosity and shear shocks

NASA Technical Reports Server (NTRS)

Syer, D.; Narayan, Ramesh

1993-01-01

Some hydrodynamical consequences of the adoption of a causal theory of viscosity are explored. Causality is introduced into the theory by letting the coefficient of viscosity go to zero as the flow velocity approaches a designated propagation speed for viscous signals. Consideration is given to a model of viscosity which has a finite propagation speed of shear information, and it is shown that it produces two kinds of shear shock. A 'pure shear shock' corresponds to a transition from a superviscous to a subviscous state with no discontinuity in the velocity. A 'mixed shear shock' has a shear transition occurring at the same location as a normal adiabatic or radiative shock. A generalized version of the Rankine-Hugoniot conditions for mixed shear shocks is derived, and self-consistent numerical solutions to a model 2D problem in which an axisymmetric radially infalling stream encounters a spinning star are presented.

A polar-drive shock-ignition design for the National Ignition Facilitya)

NASA Astrophysics Data System (ADS)

Anderson, K. S.; Betti, R.; McKenty, P. W.; Collins, T. J. B.; Hohenberger, M.; Theobald, W.; Craxton, R. S.; Delettrez, J. A.; Lafon, M.; Marozas, J. A.; Nora, R.; Skupsky, S.; Shvydky, A.

2013-05-01

Shock ignition [R. Betti et al., Phys. Rev. Lett. 98, 155001 (2007)] is being pursued as a viable option to achieve ignition on the National Ignition Facility (NIF). Shock-ignition target designs use a high-intensity laser spike at the end of a low-adiabat assembly pulse to launch a spherically convergent strong shock to ignite the hot spot of an imploding capsule. A shock-ignition target design for the NIF is presented. One-dimensional simulations indicate an ignition threshold factor of 4.1 with a gain of 58. A polar-drive beam-pointing configuration for shock-ignition experiments on the NIF at 750 kJ is proposed. The capsule design is shown to be robust to the various one- and two-dimensional effects and nonuniformities anticipated on the NIF. The target is predicted to ignite with a gain of 38 when including all anticipated levels of nonuniformity and system uncertainty.

Fast Sampling Gas Chromatography (GC) System for Speciation in a Shock Tube

DTIC Science & Technology

2016-10-31

capture similar ethylene decomposition rates for temperature-dependent shock experiments. (a) Papers published in peer-reviewed journals (N/A for none...3 GC Sampling System Validation Experiments ............................................................................... 5 Ethylene ...results for cold shock experiments, and both techniques capture similar ethylene decomposition rates for temperature-dependent shock experiments. Problem

Comment on: Negative ions, molecular electron affinity and orbital structure of cata-condensed polycyclic aromatic hydrocarbons by Rustem V. Khatymov, Mars V. Muftakhov and Pavel V. Shchukin.

PubMed

Chen, Edward S; Chen, Edward C M

2018-02-15

The anion mass spectral lifetimes for several aromatic hydrocarbons reported in the subject article were related to significantly different electron affinities. The different values are rationalized using negative ion mass spectral data. Electron affinities for polycyclic aromatic hydrocarbons are reported from the temperature dependence of unpublished electron capture detector data. These are compared with published values and the largest values are assigned to the ground state. The ground state adiabatic electron affinities: (eV) pentacene, 1.41 (3); tetracene, 1.058 (5); benz(a)pyrene, 0.82 (4); benz(a) anthracene, 0.69 (2) anthracene, 0.68 (2); and pyrene, 0.59 (1) are used to assign excited state adiabatic electron affinities: (eV) tetracene: 0.88 (4); anthracene 0.53 (1); pyrene, 0.41 (1); benz(a)anthracene, 0.39 (10); chrysene, 0.32 (1); and phenanthrene, 0.12 (2) and ground state adiabatic electron affinities: (eV) dibenz(a,j)anthracene, 0.69 (3); dibenz(a,h)anthracene, 0.68 (3); benz(e)pyrene, 0.60 (3); and picene, 0.59 (3) from experimental data. The lifetime of benz(a)pyrene is predicted to be larger than 150 μs and for benzo(c)phenanthrene and picene about 40 μs, from ground state adiabatic electron affinities. The assignments of adiabatic electron affinities of aromatic hydrocarbons determined from electron capture detector and mass spectrometric data to ground and excited states are supported by constant electronegativities. A set of consistent ground state adiabatic electron affinities for 15 polycyclic aromatic hydrocarbons is related to lifetimes from the subject article. Copyright © 2017 John Wiley & Sons, Ltd.

Analysis of temperature measurements at lead(Pb) / transparent window interface under shock compression

NASA Astrophysics Data System (ADS)

Robert, G.; Gillot, F.; Bénier, J.

2014-05-01

In this paper, we show that if the temperature bar TA obtained by a pyrometric measurement on a shock-heated material can be reached with a good precision (~5%), its transformation into a useful temperature bar TT to constrain an equation of state is not straightforward. The effects of interface, in particular the adaptation of impedance, can create a difference between bar TA and bar TT of more than 10%. This impedance correction depends on the shock adiabat of the glue, not known for thin layers of few |am but also of the equation of state of the material and of its lines of phase transition.

The Hugoniot adiabat of crystalline copper based on molecular dynamics simulation and semiempirical equation of state

NASA Astrophysics Data System (ADS)

Gubin, S. A.; Maklashova, I. V.; Mel'nikov, I. N.

2018-01-01

The molecular dynamics (MD) method was used for prediction of properties of copper under shock-wave compression and clarification of the melting region of crystal copper. The embedded atom potential was used for the interatomic interaction. Parameters of Hugonoit adiabats of solid and liquid phases of copper calculated by the semiempirical Grüneisen equation of state are consistent with the results of MD simulations and experimental data. MD simulation allows to visualize the structure of cooper on the atomistic level. The analysis of the radial distribution function and the standard deviation by MD modeling allows to predict the melting area behind the shock wave front. These MD simulation data are required to verify the wide-range equation of state of metals. The melting parameters of copper based on MD simulations and semiempirical equations of state are consistent with experimental and theoretical data, including the region of the melting point of copper.

Particle acceleration at shocks in the inner heliosphere

NASA Astrophysics Data System (ADS)

Parker, Linda Neergaard

This dissertation describes a study of particle acceleration at shocks via the diffusive shock acceleration mechanism. Results for particle acceleration at both quasi-parallel and quasi-perpendicular shocks are presented to address the question of whether there are sufficient particles in the solar wind thermal core, modeled as either a Maxwellian or kappa- distribution, to account for the observed accelerated spectrum. Results of accelerating the theoretical upstream distribution are compared to energetic observations at 1 AU. It is shown that the particle distribution in the solar wind thermal core is sufficient to explain the accelerated particle spectrum downstream of the shock, although the shape of the downstream distribution in some cases does not follow completely the theory of diffusive shock acceleration, indicating possible additional processes at work in the shock for these cases. Results show good to excellent agreement between the theoretical and observed spectral index for one third to one half of both quasi-parallel and quasi-perpendicular shocks studied herein. Coronal mass ejections occurring during periods of high solar activity surrounding solar maximum can produce shocks in excess of 3-8 shocks per day. During solar minimum, diffusive shock acceleration at shocks can generally be understood on the basis of single independent shocks and no other shock necessarily influences the diffusive shock acceleration mechanism. In this sense, diffusive shock acceleration during solar minimum may be regarded as Markovian. By contrast, diffusive shock acceleration of particles at periods of high solar activity (e.g. solar maximum) see frequent, closely spaced shocks that include the effects of particle acceleration at preceding and following shocks. Therefore, diffusive shock acceleration of particles at solar maximum cannot be modeled on the basis of diffusive shock acceleration as a single, independent shock and the process is essentially non-Markovian. A multiple shock model is developed based in part on the box model of (Protheroe and Stanev, 1998; Moraal and Axford, 1983; Ball and Kirk, 1992; Drury et al. 1999) that accelerates particles at multiple shocks and decompresses the particles between shocks via two methods. The first method of decompression is based on the that used by Melrose and Pope (1993), which adiabatically decompresses particles between shocks. The second method solves the cosmic ray transport equation and adiabatically decompresses between shocks and includes the loss of particles through convection and diffusion. The transport method allows for the inclusion of a temporal variability and thus allows for a more representative frequency distribution of shocks. The transport method of decompression and loss is used to accelerate particles at seventy-three shocks in a thirty day time period. Comparisons with observations taken at 1 AU during the same time period are encouraging as the model is able to reproduce the observed amplitude of the accelerated particles and in part the variability. This work provides the basis for developing more sophisticated models that can be applied to a suite of observations

Variable High Order Multiblock Overlapping Grid Methods for Mixed Steady and Unsteady Multiscale Viscous Flows

NASA Technical Reports Server (NTRS)

Sjogreen, Bjoern; Yee, H. C.

2007-01-01

Flows containing steady or nearly steady strong shocks in parts of the flow field, and unsteady turbulence with shocklets on other parts of the flow field are difficult to capture accurately and efficiently employing the same numerical scheme even under the multiblock grid or adaptive grid refinement framework. On one hand, sixth-order or higher shock-capturing methods are appropriate for unsteady turbulence with shocklets. On the other hand, lower order shock-capturing methods are more effective for strong steady shocks in terms of convergence. In order to minimize the shortcomings of low order and high order shock-capturing schemes for the subject flows,a multi- block overlapping grid with different orders of accuracy on different blocks is proposed. Test cases to illustrate the performance of the new solver are included.

Hot spot initiation and chemical reaction in shocked polymeric bonded explosives

NASA Astrophysics Data System (ADS)

An, Qi; Zybin, Sergey; Jaramillo-Botero, Andres; Goddard, William; Materials; Process Simulation Center, Caltech Team

2011-06-01

A polymer bonded explosive (PBX) model based on PBXN-106 is studied via molecular dynamics (MD) simulations using reactive force field (ReaxFF) under shock loading conditions. Hotspot is observed when shock waves pass through the non-planar interface of explosives and elastomers. Adiabatic shear localization is proposed as the main mechanism of hotspot ignition in PBX for high velocity impact. Our simulation also shows that the coupling of shear localization and chemical reactions at hotspot region play important rules at stress relaxtion for explosives. The phenomenon that shock waves are obsorbed by elastomers is also observed in the MD simulations. This research received supports from ARO (W911NF-05-1-0345; W911NF-08-1-0124), ONR (N00014-05-1-0778), and Los Alamos National Laboratory (LANL).

Measurements of Reduced Hydrodynamic Instability Growth in Adiabat Shaped Implosions at the NIF

NASA Astrophysics Data System (ADS)

Casey, Daniel; Macphee, Andrew; Milovich, Jose; Smalyuk, Vladimir; Clark, Dan; Robey, Harry; Peterson, Luc; Baker, Kevin; Weber, Chris

2015-11-01

Hydrodynamic instabilities can cause capsule defects and other perturbations to grow and degrade implosion performance in ignition experiments at the National Ignition Facility (NIF). Radiographic measurements of ablation front perturbation growth were performed using adiabat-shaped drives which are shown to have lower ablation front growth than the low foot drive. This is partly due to faster Richtmyer-Meshkov (RM) oscillations during the shock transit phase of the implosion moving the node in the growth factor spectrum to lower mode numbers reducing the peak growth amplitude. This is demonstrated experimentally by a reversal of the perturbation phase at higher mode numbers (120-160). These results show that the ablation front growth and fuel adiabat can be controlled somewhat-independently and are providing insight into new, more stable, ignition designs. This work was performed under the auspices of the U.S. Department of Energy by LLNL under Contract DE-AC52-07NA27344.

Cryogenic Target-Implosion Experiments on OMEGA

DOE Office of Scientific and Technical Information (OSTI.GOV)

Harding, D.R.; Meyerhofer, D.D.; Sangster, T.C.

The University of Rochester’s Laboratory for Laser Energetics has been imploding thick cryogenic targets for six years. Improvements in the Cryogenic Target Handling System and the ability to accurately design laser pulse shapes that properly time shocks and minimize electron preheat, produced high fuel areal densities in deuterium cryogenic targets (202+/-7 mg/cm^2). The areal density was inferred from the energy loss of secondary protons in the fuel (D2) shell. Targets were driven on a low final adiabat (alpha = 2) employing techniques to radially grade the adiabat (the highest adiabat at the ablation surface). The ice layer meets the target-designmore » toughness specification for DT ice of 1-um rms (all modes), while D2 ice layers average 3.0-um-rms roughness. The implosion experiments and the improvements in the quality and understanding of cryogenic targets are presented.« less

Shocks in the Early Universe.

PubMed

Pen, Ue-Li; Turok, Neil

2016-09-23

We point out a surprising consequence of the usually assumed initial conditions for cosmological perturbations. Namely, a spectrum of Gaussian, linear, adiabatic, scalar, growing mode perturbations not only creates acoustic oscillations of the kind observed on very large scales today, it also leads to the production of shocks in the radiation fluid of the very early Universe. Shocks cause departures from local thermal equilibrium as well as create vorticity and gravitational waves. For a scale-invariant spectrum and standard model physics, shocks form for temperatures 1  GeV

Modeling the Plasma Flow in the Inner Heliosheath with a Spatially Varying Compression Ratio

DOE Office of Scientific and Technical Information (OSTI.GOV)

Nicolaou, G.; Livadiotis, G.

2017-03-20

We examine a semi-analytical non-magnetic model of the termination shock location previously developed by Exarhos and Moussas. In their study, the plasma flow beyond the shock is considered incompressible and irrotational, thus the flow potential is analytically derived from the Laplace equation. Here we examine the characteristics of the downstream flow in the heliosheath in order to resolve several inconsistencies existing in the Exarhos and Moussas model. In particular, the model is modified in order to be consistent with the Rankine–Hugoniot jump conditions and the geometry of the termination shock. It is shown that a shock compression ratio varying alongmore » the latitude can lead to physically correct results. We describe the new model and present several simplified examples for a nearly spherical, strong termination shock. Under those simplifications, the upstream plasma is nearly adiabatic for large (∼100 AU) heliosheath thickness.« less

Ion acoustic shock waves in plasmas with warm ions and kappa distributed electrons and positrons

DOE Office of Scientific and Technical Information (OSTI.GOV)

Hussain, S.; Mahmood, S.; Hafeez Ur-Rehman

2013-06-15

The monotonic and oscillatory ion acoustic shock waves are investigated in electron-positron-ion plasmas (e-p-i) with warm ions (adiabatically heated) and nonthermal kappa distributed electrons and positrons. The dissipation effects are included in the model due to kinematic viscosity of the ions. Using reductive perturbation technique, the Kadomtsev-Petviashvili-Burgers (KPB) equation is derived containing dispersion, dissipation, and diffraction effects (due to perturbation in the transverse direction) in e-p-i plasmas. The analytical solution of KPB equation is obtained by employing tangent hyperbolic (Tanh) method. The analytical condition for the propagation of oscillatory and monotonic shock structures are also discussed in detail. The numericalmore » results of two dimensional monotonic shock structures are obtained for graphical representation. The dependence of shock structures on positron equilibrium density, ion temperature, nonthermal spectral index kappa, and the kinematic viscosity of ions are also discussed.« less

Numerical Solutions for Supersonic Flow of an Ideal Gas Around Blunt Two-Dimensional Bodies

NASA Technical Reports Server (NTRS)

Fuller, Franklyn B.

1961-01-01

The method described is an inverse one; the shock shape is chosen and the solution proceeds downstream to a body. Bodies blunter than circular cylinders are readily accessible, and any adiabatic index can be chosen. The lower limit to the free-stream Mach number available in any case is determined by the extent of the subsonic field, which in turn depends upon the body shape. Some discussion of the stability of the numerical processes is given. A set of solutions for flows about circular cylinders at several Mach numbers and several values of the adiabatic index is included.

Shock formation and the ideal shape of ramp compression waves

NASA Astrophysics Data System (ADS)

Swift, Damian C.; Kraus, Richard G.; Loomis, Eric N.; Hicks, Damien G.; McNaney, James M.; Johnson, Randall P.

2008-12-01

We derive expressions for shock formation based on the local curvature of the flow characteristics during dynamic compression. Given a specific ramp adiabat, calculated for instance from the equation of state for a substance, the ideal nonlinear shape for an applied ramp loading history can be determined. We discuss the region affected by lateral release, which can be presented in compact form for the ideal loading history. Example calculations are given for representative metals and plastic ablators. Continuum dynamics (hydrocode) simulations were in good agreement with the algebraic forms. Example applications are presented for several classes of laser-loading experiment, identifying conditions where shocks are desired but not formed, and where long-duration ramps are desired.

Effects of heat conduction on artificial viscosity methods for shock capturing

DOE PAGES

Cook, Andrew W.

2013-12-01

Here we investigate the efficacy of artificial thermal conductivity for shock capturing. The conductivity model is derived from artificial bulk and shear viscosities, such that stagnation enthalpy remains constant across shocks. By thus fixing the Prandtl number, more physical shock profiles are obtained, only on a larger scale. The conductivity model does not contain any empirical constants. It increases the net dissipation of a computational algorithm but is found to better preserve symmetry and produce more robust solutions for strong-shock problems.

Detonation Propagation in Slabs and Axisymmetric Rate Sticks

NASA Astrophysics Data System (ADS)

Romick, Christopher; Aslam, Tariq

Insensitive high explosives (IHE) have many benefits; however, these IHEs exhibit longer reaction zones than more conventional high explosives (HE). This makes IHEs less ideal explosives and more susceptible to edge effects as well as other performance degradation issues. Thus, there is a resulting reduction in the detonation speed within the explosive. Many HE computational models, e. g. WSD, SURF, CREST, have shock-dependent reaction rates. This dependency places a high value on having an accurate shock speed. In the common practice of shock-capturing, there is ambiguity in the shock-state due to smoothing of the shock-front. Moreover, obtaining an accurate shock speed with shock-capturing becomes prohibitively computationally expensive in multiple dimensions. The use of shock-fitting removes the ambiguity of the shock-state as it is one of the boundaries. As such, the required resolution for a given error in the detonation speed is less than with shock-capturing. This allows for further insight into performance degradation. A two-dimensional shock-fitting scheme has been developed for unconfined slabs and rate sticks of HE. The HE modeling is accomplished by Euler equations utilizing several models with single-step irreversible kinetics in slab and rate stick geometries. Department of Energy - LANL.

Pyroxenes and olivines: Structural implications of shock-wave data for high pressure phases

NASA Technical Reports Server (NTRS)

Jeanloz, R.; Ahrens, T. J.

1975-01-01

The nature of the shock-induced, high-pressure phases of olivine and pyroxene rocks is examined in the light of data for the densities of a new class of perovskite-related silicate structures. Also examined are some new Hugoniot and release adiabat data for bronzite. Reexamining available shock data for magnesian pyroxenes and olivines leads to the conclusion that they define a mixed phase (or disequilibrium) region to about the 100 GPa range, related to the kinetics of phase transformation in these silicates. By recognizing this point, certain discrepancies in previous interpretations of shock data can be explained. A set of theoretical Hugonoits for pyroxene and olivine stoichiometry, perovskite-bearing assemblages was constructed based on their properties deduced from high-pressure work, showing that the shock data is compatible with transformations to perovskites in the 45-7GPa region. Finally, the shock data indicate very similar properties for olivine and pyroxene at high pressures making them both equally likely candidates for the lower mantle.

The analytical solution of the problem of a shock focusing in a gas for one-dimensional case

NASA Astrophysics Data System (ADS)

Shestakovskaya, E. S.; Magazov, F. G.

2018-03-01

The analytical solution of the problem of an imploding shock wave in the vessel with an impermeable wall is constructed for the cases of planar, cylindrical and spherical symmetry. The negative velocity is set at the vessel boundary. The velocity of cold ideal gas is zero. At the initial time the shock spreads from this point into the center of symmetry. The boundary moves under the particular law which conforms to the movement of the shock. In Euler variables it moves but in Lagrangian variables its trajectory is a vertical line. Equations that determine the structure of the gas flow between the shock front and the boundary as a function of time and the Lagrangian coordinate as well as the dependence of the entropy on the shock wave velocity are obtained. Self-similar coefficients and corresponding critical values of self-similar coordinates were found for a wide range of adiabatic index. The problem is solved for Lagrangian coordinates.

Explicit and implicit compact high-resolution shock-capturing methods for multidimensional Euler equations 1: Formulation

NASA Technical Reports Server (NTRS)

Yee, H. C.

1995-01-01

Two classes of explicit compact high-resolution shock-capturing methods for the multidimensional compressible Euler equations for fluid dynamics are constructed. Some of these schemes can be fourth-order accurate away from discontinuities. For the semi-discrete case their shock-capturing properties are of the total variation diminishing (TVD), total variation bounded (TVB), total variation diminishing in the mean (TVDM), essentially nonoscillatory (ENO), or positive type of scheme for 1-D scalar hyperbolic conservation laws and are positive schemes in more than one dimension. These fourth-order schemes require the same grid stencil as their second-order non-compact cousins. One class does not require the standard matrix inversion or a special numerical boundary condition treatment associated with typical compact schemes. Due to the construction, these schemes can be viewed as approximations to genuinely multidimensional schemes in the sense that they might produce less distortion in spherical type shocks and are more accurate in vortex type flows than schemes based purely on one-dimensional extensions. However, one class has a more desirable high-resolution shock-capturing property and a smaller operation count in 3-D than the other class. The extension of these schemes to coupled nonlinear systems can be accomplished using the Roe approximate Riemann solver, the generalized Steger and Warming flux-vector splitting or the van Leer type flux-vector splitting. Modification to existing high-resolution second- or third-order non-compact shock-capturing computer codes is minimal. High-resolution shock-capturing properties can also be achieved via a variant of the second-order Lax-Friedrichs numerical flux without the use of Riemann solvers for coupled nonlinear systems with comparable operations count to their classical shock-capturing counterparts. The simplest extension to viscous flows can be achieved by using the standard fourth-order compact or non-compact formula for the viscous terms.

Gas dynamic and force effects of a solid particle in a shock wave in air

NASA Astrophysics Data System (ADS)

Obruchkova, L. R.; Baldina, E. G.; Efremov, V. P.

2017-03-01

Shock wave interaction with an adiabatic solid microparticle is numerically simulated. In the simulation, the shock wave is initiated by the Riemann problem with instantaneous removal of a diaphragm between the high- and low-pressure chambers. The calculation is performed in the two-dimensional formulation using the ideal gas equation of state. The left end of the tube is impermeable, while outflow from the right end is permitted. The particle is assumed to be motionless, impermeable, and adiabatic, and the simulation is performed for time intervals shorted than the time of velocity and temperature relaxation of the particle. The numerical grid is chosen for each particle size to ensure convergence. For each particle size, the calculated hydraulic resistance coefficient describing the particle force impact on the flow is compared with that obtained from the analytical Stokes formula. It is discovered that the Stokes formula can be used for calculation of hydraulic resistance of a motionless particle in a shock wave flow. The influence of the particle diameter on the flow perturbation behind the shock front is studied. Specific heating of the flow in front of the particle is calculated and a simple estimate is proposed. The whole heated region is divided by the acoustic line into the subsonic and supersonic regions. It is demonstrated that the main heat generated by the particle in the flow is concentrated in the subsonic region. The calculations are performed using two different 2D hydro codes. The energy release in the flow induced by the particle is compared with the maximum possible heating at complete termination of the flow. The results can be used for estimating the possibility of gas ignition in front of the particle by a shock wave whose amplitude is insufficient for initiating detonation in the absence of a particle.

The Accuracy of Shock Capturing in Two Spatial Dimensions

NASA Technical Reports Server (NTRS)

Carpenter, Mark H.; Casper, Jay H.

1997-01-01

An assessment of the accuracy of shock capturing schemes is made for two-dimensional steady flow around a cylindrical projectile. Both a linear fourth-order method and a nonlinear third-order method are used in this study. It is shown, contrary to conventional wisdom, that captured two-dimensional shocks are asymptotically first-order, regardless of the design accuracy of the numerical method. The practical implications of this finding are discussed in the context of the efficacy of high-order numerical methods for discontinuous flows.

Shock ignition of thermonuclear fuel with high areal density.

PubMed

Betti, R; Zhou, C D; Anderson, K S; Perkins, L J; Theobald, W; Solodov, A A

2007-04-13

A novel method by C. Zhou and R. Betti [Bull. Am. Phys. Soc. 50, 140 (2005)] to assemble and ignite thermonuclear fuel is presented. Massive cryogenic shells are first imploded by direct laser light with a low implosion velocity and on a low adiabat leading to fuel assemblies with large areal densities. The assembled fuel is ignited from a central hot spot heated by the collision of a spherically convergent ignitor shock and the return shock. The resulting fuel assembly features a hot-spot pressure greater than the surrounding dense fuel pressure. Such a nonisobaric assembly requires a lower energy threshold for ignition than the conventional isobaric one. The ignitor shock can be launched by a spike in the laser power or by particle beams. The thermonuclear gain can be significantly larger than in conventional isobaric ignition for equal driver energy.

Interplanetary shocks, Plasma waves and turbulence, Kinetic waves and instabilities, STEREO spacecraft

NASA Astrophysics Data System (ADS)

Cohen, Z.; Breneman, A. W.; Cattell, C. A.; Davis, L.; Grul, P.; Kersten, K.; Wilson, L. B., III

2017-12-01

Determining the role of plasma waves in providing energy dissipation at shock waves is of long-standing interest. Interplanetary (IP) shocks serve as a large database of low Mach number shocks. We examine electric field waveforms captured by the Time Domain Sampler (TDS) on the STEREO spacecraft during the ramps of IP shocks, with emphasis on captures lasting 2.1 seconds. Previous work has used captures of shorter duration (66 and 131 ms on STEREO, and 17 ms on WIND), which allowed for observation of waves with maximum (minimum) frequencies of 125 kHz (15 Hz), 62.5 kHz (8 Hz), and 60 kHz (59 Hz), respectively. The maximum frequencies are comparable to 2-8 times the plasma frequency in the solar wind, enabling observation of Langmuir waves, ion acoustic, and some whistler-mode waves. The 2 second captures resolve lower frequencies ( few Hz), which allows us to analyze packet structure of the whistler-mode waves and some ion acoustic waves. The longer capture time also improves the resolvability of simultaneous wave modes and of waves with frequencies on the order of 10s of Hz. Langmuir waves, however, cannot be identified at this sampling rate, since the plasma frequency is usually higher than 3.9 kHz. IP shocks are identified from multiple databases (Helsinki heliospheric shock database at http://ipshocks.fi, and the STEREO level 3 shock database at ftp://stereoftp.nascom.nasa.gov/pub/ins_data/impact/level3/). Our analysis focuses on TDS captures in shock ramp regions, with ramp durations determined from magnetic field data taken at 8 Hz. Software is used to identify multiple wave modes in any given capture and classify waves as Langmuir, ion acoustic, whistler, lower hybrid, electron cyclotron drift instability, or electrostatic solitary waves. Relevant frequencies are determined from density and magnetic field data collected in situ. Preliminary results suggest that large amplitude (∼ 5 mV/m) ion acoustic waves are most prevalent in the ramp, in agreement with Wilson, et al. Other modes are also observed. Statistical results will be presented and compared with previous studies and theoretical predictions.

STEREO Observations of Waves in the Ramp Regions of Interplanetary Shocks

NASA Astrophysics Data System (ADS)

Cohen, Z.; Breneman, A. W.; Cattell, C. A.; Davis, L.; Grul, P.; Kersten, K.; Wilson, L. B., III

2017-12-01

Determining the role of plasma waves in providing energy dissipation at shock waves is of long-standing interest. Interplanetary (IP) shocks serve as a large database of low Mach number shocks. We examine electric field waveforms captured by the Time Domain Sampler (TDS) on the STEREO spacecraft during the ramps of IP shocks, with emphasis on captures lasting 2.1 seconds. Previous work has used captures of shorter duration (66 and 131 ms on STEREO, and 17 ms on WIND), which allowed for observation of waves with maximum (minimum) frequencies of 125 kHz (15 Hz), 62.5 kHz (8 Hz), and 60 kHz (59 Hz), respectively. The maximum frequencies are comparable to 2-8 times the plasma frequency in the solar wind, enabling observation of Langmuir waves, ion acoustic, and some whistler-mode waves. The 2 second captures resolve lower frequencies ( few Hz), which allows us to analyze packet structure of the whistler-mode waves and some ion acoustic waves. The longer capture time also improves the resolvability of simultaneous wave modes and of waves with frequencies on the order of 10s of Hz. Langmuir waves, however, cannot be identified at this sampling rate, since the plasma frequency is usually higher than 3.9 kHz. IP shocks are identified from multiple databases (Helsinki heliospheric shock database at http://ipshocks.fi, and the STEREO level 3 shock database at ftp://stereoftp.nascom.nasa.gov/pub/ins_data/impact/level3/). Our analysis focuses on TDS captures in shock ramp regions, with ramp durations determined from magnetic field data taken at 8 Hz. Software is used to identify multiple wave modes in any given capture and classify waves as Langmuir, ion acoustic, whistler, lower hybrid, electron cyclotron drift instability, or electrostatic solitary waves. Relevant frequencies are determined from density and magnetic field data collected in situ. Preliminary results suggest that large amplitude (≥ 5 mV/m) ion acoustic waves are most prevalent in the ramp, in agreement with Wilson, et al. Other modes are also observed. Statistical results will be presented and compared with previous studies and theoretical predictions.

Oscillations of a standing shock wave generated by the Richtmyer-Meshkov instability

NASA Astrophysics Data System (ADS)

Mikaelian, Karnig O.

2016-07-01

In a typical Richtmyer-Meshkov experiment a fast moving flat shock strikes a stationary perturbed interface between fluids A and B creating a transmitted and a reflected shock, both of which are perturbed. We propose shock tube experiments in which the reflected shock is stationary in the laboratory. Such a standing perturbed shock undergoes well-known damped oscillations. We present the conditions required for producing such a standing shock wave, which greatly facilitates the measurement of the oscillations and their rate of damping. We define a critical density ratio Rcritical, in terms of the adiabatic indices of the two fluids, and a critical Mach number Mscritical of the incident shock wave, which produces a standing reflected wave. If the initial density ratio R of the two fluids is less than Rcritical then a standing shock wave is possible at Ms=Mscritical . Otherwise a standing shock is not possible and the reflected wave always moves in the direction opposite the incident shock. Examples are given for present-day operating shock tubes with sinusoidal or inclined interfaces. We consider the effect of viscosity, which affects the damping rate of the oscillations. We point out that nonlinear bubble and spike amplitudes depend relatively weakly on the viscosity of the fluids and that the interface area is a better diagnostic.

Adiabatic theory in regions of strong field gradients. [in magnetosphere

NASA Technical Reports Server (NTRS)

Whipple, E. C.; Northrop, T. G.; Birmingham, T. J.

1986-01-01

The theory for the generalized first invariant for adiabatic motion of charged particles in regions where there are large gradients in magnetic or electric fields is developed. The general condition for an invariant to exist in such regions is that the potential well in which the particle oscillates change its shape slowly as the particle drifts. It is shown how the Kruskal (1962) procedure can be applied to obtain expressions for the invariant and for drift velocities that are asymptotic in a smallness parameter epsilon. The procedure is illustrated by obtaining the invariant and drift velocities for particles traversing a perpendicular shock, and the generalized invariant is compared with the magnetic moment, and the drift orbits with the actual orbits, for a particular case. In contrast to the magnetic moment, the generalized first invariant is better for large gyroradii (large kinetic energies) than for small gyroradii. Expressions for the invariant when an electrostatic potential jump is imposed across the perpendicular shock, and when the particle traverses a rotational shear layer with a small normal component of the magnetic field are given.

On the Magnetospheric Heating Problem

NASA Astrophysics Data System (ADS)

Nykyri, K.; Moore, T.; Dimmock, A. P.; Ma, X.; Johnson, J.; Delamere, P. A.

2016-12-01

In the Earth's magnetosphere the specific entropy, increases by approximately two orders of magnitude when transitioning from the magnetosheath into the magnetosphere. However, the origin of this non-adiabatic heating is not well understood. In addition, there exists a dawn-dusk temperature asymmetry in the flanks of the plasma sheet - the cold component ions are hotter by 30-40% at the dawnside plasma sheet compared to the duskside plasma sheet. Our recent statistical study of magnetosheath temperatures using 7 years of THEMIS data indicates that ion magnetosheath temperatures downstream of quasi-parallel (dawn-flank for the Parker-Spiral IMF) bow shock are only 15 percent higher than downstream of the quasi-perpendicular shock. This magnetosheath temperature asymmetry is therefore inadequate to cause the observed level of the plasma sheet temperature asymmetry. In this presentation we address the origin of non-adiabatic heating from the magnetosheath into the plasma sheet by utilizing small Cluster spacecraft separations, 9 years of statistical THEMIS data as well as Hall-MHD and hybrid simulations. We present evidence of a new physical mechanism capable of cross-scale energy transport at the flank magnetopause with strong contributions to the non-adiabatic heating observed between the magnetosheath and plasma sheet. This same heating mechanism may occur and drive asymmetries also in the magnetospheres of gas giants: Jupiter and Saturn, as well as play role elsewhere in the universe where significant flow shears are present such as in the solar corona, and other astrophysical and laboratory plasmas.

Investigation of Performance of Axial-Flow Compressor of XT-46 Turbine-Propeller Engine. I - Preliminary Investigation at 50-,70-, and 100-Percent Design Equivalent Speed

NASA Technical Reports Server (NTRS)

Creagh, John W.R.; Sandercrock, Donald M.

1950-01-01

An investigation is being conducted to determine the performance of the 12-stage axial-flow compressor of the XT-46 turbine-propeller engine. This compressor was designed to produce a pressure ratio of 9 at an adiabatic efficiency of 0.86. The design pressure ratios per stage were considerably greater than any employed in current aircraft gas-turbine engines using this type of compressor. The compressor performance was evaluated at two stations. The station near the entrance section of the combustors indicated a peak pressure ratio of 6.3 at an adiabatic efficiency of 0.63 for a corrected weight flow of 23.1 pounds per second. The other, located one blade-chord downstream of the last stator row, indicated a peak pressure ratio of 6.97 at an adiabatic efficiency of 0.81 for a corrected weight flow of 30.4 pounds per second. The difference in performance obtained at the two stations is attributed to shock waves in the vicinity of the last stator row. These shock waves and the accompanying flow choking, together with interstage circulatory flows, shift the compressor operating curves into the region where surge would normally occur. The inability of the compressor to meet design pressure ratio is probably due to boundary-layer buildup in the last stages, which cause axial velocities greater than design values that, in turn, adversely affect the angles of attack and turning angles in these blade rows.

A bulk viscosity approach for shock capturing on unstructured grids

NASA Astrophysics Data System (ADS)

Shoeybi, Mohammad; Larsson, Nils Johan; Ham, Frank; Moin, Parviz

2008-11-01

The bulk viscosity approach for shock capturing (Cook and Cabot, JCP, 2005) augments the bulk part of the viscous stress tensor. The intention is to capture shock waves without dissipating turbulent structures. The present work extends and modifies this method for unstructured grids. We propose a method that properly scales the bulk viscosity with the grid spacing normal to the shock for unstructured grid for which the shock is not necessarily aligned with the grid. The magnitude of the strain rate tensor used in the original formulation is replaced with the dilatation, which appears to be more appropriate in the vortical turbulent flow regions (Mani et al., 2008). The original form of the model is found to have an impact on dilatational motions away form the shock wave, which is eliminated by a proposed localization of the bulk viscosity. Finally, to allow for grid adaptation around shock waves, an explicit/implicit time advancement scheme has been developed that adaptively identifies the stiff regions. The full method has been verified with several test cases, including 2D shock-vorticity entropy interaction, homogenous isotropic turbulence, and turbulent flow over a cylinder.

Aspherical Supernovae and Oblique Shock Breakout

NASA Astrophysics Data System (ADS)

Afsariardchi, Niloufar; Matzner, Christopher D.

2017-02-01

In an aspherical supernova explosion, shock emergence is not simultaneous and non-radial flows develop near the stellar surface. Oblique shock breakouts tend to be easily developed in compact progenitors like stripped-envelop core collapse supernovae. According to Matzner et al. (2013), non-spherical explosions develop non-radial flows that alters the observable emission and radiation of a supernova explosion. These flows can limit ejecta speed, change the distribution of matter and heat of the ejecta, suppress the breakout flash, and most importantly engender collisions outside the star. We construct a global numerical FLASH hydrodynamic simulation in a two dimensional spherical coordinate, focusing on the non-relativistic, adiabatic limit in a polytropic envelope to see how these fundamental differences affect the early light curve of core-collapse SNe.

Effects of molecular dissociation on the hydrogen equation of state

NASA Astrophysics Data System (ADS)

Bonev, Stanimir; Schwegler, Eric; Galli, Giulia; Gygi, Francois

2002-03-01

It has been suggested recently(François Gygi and G. Galli, submitted to Phys. Rev. Lett.) that the physical mechanism behind the larger compressibility of liquid deuterium observed in laser shock experiments as compared to ab initio simulations may be related to shock-induced electronic excitations. A possible result of such non-adiabatic processes is hindering of the molecular dissociation. This has motivated us to study the importance of molecular dissociation on the hydrogen equation of state. To this end, we have carried out ab initio molecular dynamics simulations of liquid deuterium where intramolecular dissociation is prevented by the use of bond length contraints. Simulations at both fixed thermodynamic conditions and dynamical simulations of shocked deuterium will be discussed.

The Strength and Radial Profile of the Coronal Magnetic Field from the Standoff Distance of a Coronal Mass Ejection-Driven Shock

NASA Technical Reports Server (NTRS)

Gopalswamy, Nat; Yashiro, Seiji

2011-01-01

We determine the coronal magnetic field strength in the heliocentric distance range 6-23 solar radii (Rs) by measuring the shock standoff distance and the radius of curvature of the flux rope during the 2008 March 25 coronal mass ejection imaged by white-light coronagraphs. Assuming the adiabatic index, we determine the Alfven Mach number, and hence the Alfven speed in the ambient medium using the measured shock speed. By measuring the upstream plasma density using polarization brightness images, we finally get the magnetic field strength upstream of the shock. The estimated magnetic field decreases from approximately 48 mG around 6 Rs to 8 mG at 23 Rs. The radial profile of the magnetic field can be described by a power law in agreement with other estimates at similar heliocentric distances.

Driving mechanism of unsteady separation shock motion in hypersonic interactive flow

NASA Technical Reports Server (NTRS)

Dolling, D. S.; Narlo, J. C., II

1987-01-01

Wall pressure fluctuations were measured under the steady separation shock waves in Mach 5 turbulent interactions induced by unswept circular cylinders on a flat plate. The wall temperature was adiabatic. A conditional sampling algorithm was developed to examine the statistics of the shock wave motion. The same algorithm was used to examine data taken in earlier studies in the Princeton University Mach 3 blowdown tunnel. In these earlier studies, hemicylindrically blunted fins of different leading-edge diameters were tested in boundary layers which developed on the tunnel floor and on a flat plate. A description of the algorithm, the reasons why it was developed and the sensitivity of the results to the threshold settings, are discussed. The results from the algorithm, together with cross correlations and power spectral density estimates suggests that the shock motion is driven by the low-frequency unsteadiness of the downstream separated, vortical flow.

3D Simulations of the ``Keyhole'' Hohlraum for Shock Timing on NIF

NASA Astrophysics Data System (ADS)

Robey, H. F.; Marinak, M. M.; Munro, D. H.; Jones, O. S.

2007-11-01

Ignition implosions planned for the National Ignition Facility (NIF) require a pulse shape with a carefully designed series of steps, which launch a series of shocks through the ablator and DT fuel. The relative timing of these shocks must be tuned to better than +/- 100ps to maintain the DT fuel on a sufficiently low adiabat. To meet these requirements, pre-ignition tuning experiments using a modified hohlraum geometry are being planned. This modified geometry, known as the ``keyhole'' hohlraum, adds a re-entrant gold cone, which passes through the hohlraum and capsule walls, to provide an optical line-of-sight to directly measure the shocks as they break out of the ablator. In order to assess the surrogacy of this modified geometry, 3D simulations using HYDRA [1] have been performed. The drive conditions and the resulting effect on shock timing in the keyhole hohlraum will be compared with the corresponding results for the standard ignition hohlraum. [1] M.M. Marinak, et al., Phys. Plasmas 8, 2275 (2001).

The unifying role of dissipative action in the dynamic failure of solids

DOE PAGES

Grady, Dennis

2015-05-19

Dissipative action, the product of dissipation energy and transport time, is fundamental to the dynamic failure of solids. Invariance of the dissipative action underlies the fourth-power nature of structured shock waves observed in selected solid metals and compounds. Dynamic failure through shock compaction, tensile spall and adiabatic shear are also governed by a constancy of the dissipative action. This commonality underlying the various modes of dynamic failure is described and leads to deeper insights into failure of solids in the intense shock wave event. These insights are in turn leading to a better understanding of the shock deformation processes underlyingmore » the fourth-power law. Experimental result and material models encompassing the dynamic failure of solids are explored for the purpose of demonstrating commonalities leading to invariance of the dissipation action. As a result, calculations are extended to aluminum and uranium metals with the intent of predicting micro-scale energetics and spatial scales in the structured shock wave.« less

Oscillations of a standing shock wave generated by the Richtmyer-Meshkov instability

DOE PAGES

Mikaelian, Karnig O.

2016-07-13

In a typical Richtmyer-Meshkov experiment a fast moving flat shock strikes a stationary perturbed interface between fluids A and B creating a transmitted and a reflected shock, both of which are perturbed. We propose shock tube experiments in which the reflected shock is stationary in the laboratory. Such a standing perturbed shock undergoes well-known damped oscillations. We present the conditions required for producing such a standing shock wave, which greatly facilitates the measurement of the oscillations and their rate of damping. We define a critical density ratio R critical, in terms of the adiabatic indices of the two fluids, andmore » a critical Mach number M critical s of the incident shock wave, which produces a standing reflected wave. If the initial density ratio R of the two fluids is less than R critical then a standing shock wave is possible at M s=M critical s. Otherwise a standing shock is not possible and the reflected wave always moves in the direction opposite the incident shock. Examples are given for present-day operating shock tubes with sinusoidal or inclined interfaces. We consider the effect of viscosity, which affects the damping rate of the oscillations. Furthermore, we point out that nonlinear bubble and spike amplitudes depend relatively weakly on the viscosity of the fluids and that the interface area is a better diagnostic.« less

Electron capture cross sections by O+ from atomic He

NASA Astrophysics Data System (ADS)

Joseph, Dwayne C.; Saha, Bidhan C.

2009-11-01

The adiabatic representation is used in both the quantal and semi classical molecular orbital close coupling methods (MOCC) to evaluate charge exchange cross sections. Our results show good agreement with experimental cross sections

Investigation of the crater-like microdefects induced by laser shock processing with aluminum foil as absorbent layer

NASA Astrophysics Data System (ADS)

Ye, Y. X.; Xuan, T.; Lian, Z. C.; Feng, Y. Y.; Hua, X. J.

2015-06-01

This paper reports that 3D crater-like microdefects form on the metal surface when laser shock processing (LSP) is applied. LSP was conducted on pure copper block using the aluminum foil as the absorbent material and water as the confining layer. There existed the bonding material to attach the aluminum foil on the metal target closely. The surface morphologies and metallographs of copper surfaces were characterized with 3D profiler, the optical microscopy (OM) or the scanning electron microscopy (SEM). Temperature increases of metal surface due to LSP were evaluated theoretically. It was found that, when aluminum foil was used as the absorbent material, and if there existed air bubbles in the bonding material, the air temperatures within the bubbles rose rapidly because of the adiabatic compression. So at the locations of the air bubbles, the metal materials melted and micromelting pool formed. Then under the subsequent expanding of the air bubbles, a secondary shock wave was launched against the micromelting pool and produced the crater-like microdefects on the metal surface. The temperature increases due to shock heat and high-speed deformation were not enough to melt the metal target. The temperature increase induced by the adiabatic compression of the air bubbles may also cause the gasification of the metal target. This will also help form the crater-like microdefects. The results of this paper can help to improve the surface quality of a metal target during the application of LSP. In addition, the results provide another method to fabricate 3D crater-like dents on metal surface. This has a potential application in mechanical engineering.

Heat Pipe Vapor Dynamics. Ph.D. Thesis

NASA Technical Reports Server (NTRS)

Issacci, Farrokh

1990-01-01

The dynamic behavior of the vapor flow in heat pipes is investigated at startup and during operational transients. The vapor is modeled as two-dimensional, compressible viscous flow in an enclosure with inflow and outflow boundary conditions. For steady-state and operating transients, the SIMPLER method is used. In this method a control volume approach is employed on a staggered grid which makes the scheme very stable. It is shown that for relatively low input heat fluxes the compressibility of the vapor flow is low and the SIMPLER scheme is suitable for the study of transient vapor dynamics. When the input heat flux is high or the process under a startup operation starts at very low pressures and temperatures, the vapor is highly compressible and a shock wave is created in the evaporator. It is shown that for a wide range of input heat fluxes, the standard methods, including the SIMPLER scheme, are not suitable. A nonlinear filtering technique, along with the centered difference scheme, are then used for shock capturing as well as for the solution of the cell Reynolds-number problem. For high heat flux, the startup transient phase involves multiple shock reflections in the evaporator region. Each shock reflection causes a significant increase in the local pressure and a large pressure drop along the heat pipe. Furthermore, shock reflections cause flow reversal in the evaporation region and flow circulations in the adiabatic region. The maximum and maximum-averaged pressure drops in different sections of the heat pipe oscillate periodically with time because of multiple shock reflections. The pressure drop converges to a constant value at steady state. However, it is significantly higher than its steady-state value at the initiation of the startup transient. The time for the vapor core to reach steady-state condition depends on the input heat flux, the heat pipe geometry, the working fluid, and the condenser conditions. However, the vapor transient time, for an Na-filled heat pipe is on the order of seconds. Depending on the time constant for the overall system, the vapor transient time may be very short. Therefore, the vapor core may be assumed to be quasi-steady in the transient analysis of a heat pipe operation.

Interaction of a shock with a longitudinal vortex

NASA Technical Reports Server (NTRS)

Erlebacher, Gordon; Hussaini, M. Y.; Shu, Chi-Wang

1996-01-01

In this paper we study the shock/longitudinal vortex interaction problem in axisymmetric geometry. Linearized analysis for small vortex strength is performed, and compared with results from a high order axisymmetric shock-fitted Euler solution obtained for this purpose. It is confirmed that for weak vortices, predictions from linear theory agree well with results from nonlinear numerical simulations at the shock location. To handle very strong longitudinal vortices, which may ultimately break the shock, we use an axisymmetric high order essentially non-oscillatory (ENO) shock capturing scheme. Comparison of shock-captured and shock-fitted results are performed in their regions of common validity. We also study the vortex breakdown as a function of Mach number ranging from 1.3 to 10, thus extending the range of existing results. For vortex strengths above a critical value. a triple point forms on the shock and a secondary shock forms to provide the necessary deceleration so that the fluid velocity can adjust to downstream conditions at the shock.

DOUBLE POWER LAWS IN THE EVENT-INTEGRATED SOLAR ENERGETIC PARTICLE SPECTRUM

DOE Office of Scientific and Technical Information (OSTI.GOV)

Zhao, Lulu; Zhang, Ming; Rassoul, Hamid K., E-mail: lzhao@fit.edu

2016-04-10

A double power law or a power law with exponential rollover at a few to tens of MeV nucleon{sup −1} of the event-integrated differential spectra has been reported in many solar energetic particle (SEP) events. The rollover energies per nucleon of different elements correlate with a particle's charge-to-mass ratio (Q/A). The probable causes are suggested as residing in shock finite lifetimes, shock finite sizes, shock geometry, and an adiabatic cooling effect. In this work, we conduct a numerical simulation to investigate a particle's transport process in the inner heliosphere. We solve the focused transport equation using a time-backward Markov stochasticmore » approach. The convection, magnetic focusing, adiabatic cooling effect, and pitch-angle scattering are included. The effects that the interplanetary turbulence imposes on the shape of the resulting SEP spectra are examined. By assuming a pure power-law differential spectrum at the Sun, a perfect double-power-law feature with a break energy ranging from 10 to 120 MeV nucleon{sup −1} is obtained at 1 au. We found that the double power law of the differential energy spectrum is a robust result of SEP interplanetary propagation. It works for many assumptions of interplanetary turbulence spectra that give various forms of momentum dependence of a particle's mean free path. The different spectral shapes in low-energy and high-energy ends are not just a transition from the convection-dominated propagation to diffusion-dominated propagation.« less

Vapor Flow Patterns During a Start-Up Transient in Heat Pipes

NASA Technical Reports Server (NTRS)

Issacci, F.; Ghoniem, N, M.; Catton, I.

1996-01-01

The vapor flow patterns in heat pipes are examined during the start-up transient phase. The vapor core is modelled as a channel flow using a two dimensional compressible flow model. A nonlinear filtering technique is used as a post process to eliminate the non-physical oscillations of the flow variables. For high-input heat flux, multiple shock reflections are observed in the evaporation region. The reflections cause a reverse flow in the evaporation and circulations in the adiabatic region. Furthermore, each shock reflection causes a significant increase in the local pressure and a large pressure drop along the heat pipe.

Improved Performance of High Areal Density Indirect Drive Implosions at the National Ignition Facility using a Four-Shock Adiabat Shaped Drive

DOE PAGES

Casey, D. T.; Milovich, J. L.; Smalyuk, V. A.; ...

2015-09-01

Hydrodynamic instabilities can cause capsule defects and other perturbations to grow and degrade implosion performance in ignition experiments at the National Ignition Facility (NIF). Here, we show the first experimental demonstration that a strong unsupported first shock in indirect drive implosions at the NIF reduces ablation front instability growth leading to a 3 to 10 times higher yield with fuel ρR > 1 g=cm 2. This work shows the importance of ablation front instability growth during the National Ignition Campaign and may provide a path to improved performance at the high compression necessary for ignition.

Addition of Improved Shock-Capturing Schemes to OVERFLOW 2.1

NASA Technical Reports Server (NTRS)

Burning, Pieter G.; Nichols, Robert H.; Tramel, Robert W.

2009-01-01

Existing approximate Riemann solvers do not perform well when the grid is not aligned with strong shocks in the flow field. Three new approximate Riemann algorithms are investigated to improve solution accuracy and stability in the vicinity of strong shocks. The new algorithms are compared to the existing upwind algorithms in OVERFLOW 2.1. The new algorithms use a multidimensional pressure gradient based switch to transition to a more numerically dissipative algorithm in the vicinity of strong shocks. One new algorithm also attempts to artificially thicken captured shocks in order to alleviate the errors in the solution introduced by "stair-stepping" of the shock resulting from the approximate Riemann solver. This algorithm performed well for all the example cases and produced results that were almost insensitive to the alignment of the grid and the shock.

In-Flight Boundary-Layer Transition on a Large Flat Plate at Supersonic Speeds

NASA Technical Reports Server (NTRS)

Banks, Daniel W.; Fredericks, Michael Alan; Tracy, Richard R.; Matisheck, Jason R.; Vanecek, Neal D.

2012-01-01

A flight experiment was conducted to investigate the pressure distribution, local flow conditions, and boundary-layer transition characteristics on a large flat plate in flight at supersonic speeds up to Mach 2.0. The primary objective of the test was to characterize the local flow field in preparation for future tests of a high Reynolds number natural laminar flow test article. The tests used a F-15B testbed aircraft with a bottom centerline mounted test fixture. A second objective was to determine the boundary-layer transition characteristics on the flat plate and the effectiveness of using a simplified surface coating for future laminar flow flight tests employing infrared thermography. Boundary-layer transition was captured using an onboard infrared imaging system. The infrared imagery was captured in both analog and digital formats. Surface pressures were measured with electronically scanned pressure modules connected to 60 surface-mounted pressure orifices. The local flow field was measured with five 5-hole conical probes mounted near the leading edge of the test fixture. Flow field measurements revealed the local flow characteristics including downwash, sidewash, and local Mach number. Results also indicated that the simplified surface coating did not provide sufficient insulation from the metallic structure, which likely had a substantial effect on boundary-layer transition compared with that of an adiabatic surface. Cold wall conditions were predominant during the acceleration to maximum Mach number, and warm wall conditions were evident during the subsequent deceleration. The infrared imaging system was able to capture shock wave impingement on the surface of the flat plate in addition to indicating laminar-to-turbulent boundary-layer transition.

Modeling and Laboratory Investigations of Radiative Shocks

NASA Astrophysics Data System (ADS)

Grun, Jacob; Laming, J. Martin; Manka, Charles; Moore, Christopher; Jones, Ted; Tam, Daniel

2001-10-01

Supernova remnants are often inhomogeneous, with knots or clumps of material expanding in ambient plasma. This structure may be initiated by hydrodynamic instabilities occurring during the explosion, but it may plausibly be amplified by instabilities of the expanding shocks such as, for example, corrugation instabilities described by D’yakov in 1954, Vishniac in 1983, and observed in the laboratory by Grun et al. in 1991. Shock instability can occur when radiation lowers the effective adiabatic index of the gas. In view of the difficulty of modeling radiation in non-equilibrium plasmas, and the dependence of shock instabilities on such radiation, we are performing a laboratory experiment to study radiative shocks. The shocks are generated in a miniature, laser-driven shock tube. The gas density inside the tube at any instant in time is measured using time and space-resolved interferometry, and the emission spectrum of the gas is measured with time-resolved spectroscopy. We simulate the experiment with a 1D code that models time dependent post-shock ionization and non-equilibrium radiative cooling. S. P. D’yakov, Zhurnal Eksperimentalnoi Teoreticheskoi Fiziki 27, 288 (1954); see also section 90 in L.D. Landau and E.M. Lifshitz, Fluid Mechanics (Butterworth-Heinemann 1987); E.T. Vishniac, Astrophys. J. 236, 880 (1983); J. Grun, et al., Phys. Rev. Lett., 66, 2738 (1991)

Shock compression of Fe-FeS mixture up to 204 GPa

NASA Astrophysics Data System (ADS)

Huang, Haijun; Wu, Shijie; Hu, Xiaojun; Wang, Qingsong; Wang, Xiang; Fei, Yingwei

2013-02-01