TIME-DEPENDENT SIMULATION OF OBLIQUE MHD COSMIC-RAY SHOCKS USING THE 2-FLUID MODEL

Abstract

Using a new, second-order accurate numerical method we present dynamical simulations of oblique MHD cosmic-ray (CR)-modified plane shock evolution. Most of the calculations are done with a two-fluid model for diffusive shock acceleration, but we provide also comparisons between a typical shock computed that way against calculations carried out using the more complete, momentum-dependent, diffusion-advection equation. We also illustrate a test showing that these simulations evolve to dynamical equilibra consistent with previously published steady state analytic calculations for such shocks. In order to improve understanding of the dynamical role of magnetic fields in shocks modified by CR pressure we have explored for time asymptotic states the parameter space of upstream fast mode Mach number, M(f), and plasma beta. We compile the results into maps of dynamical steady state CR acceleration efficiency, epsilon(c). Since the models are simplifications, such maps should not be used to predict quantitatively epsilon(c) in real shocks; however, they are internally consistent, so that they can enable us to compare various competing dynamical effects. The maps, along with additional numerical experiments, show that epsilon(c) is reduced through the action of compressive work on tangential magnetic fields in CR-MHD shocks. Thus epsilon(c) in low beta, moderate M(f) shocks tends to be smaller in quasi-perpendicular shocks than it would be for high beta shocks of the same M(f). This result supports earlier conclusions that strong, oblique magnetic fields inhibit diffusive shock acceleration. For quasi-parallel shocks with beta < 1, on the other hand, epsilon(c) seems to be increased at a given M(f) when compared to high beta shocks. The apparent contradiction to the first conclusion results, however, from the fact that for small beta quasi-parallel shocks, the fast mode Mach number is not a good measure of compression through the shock. That is better reflected in the sonic Mach number, which is greater in these instances, Acceleration efficiencies for high and low beta having comparable sonic Mach numbers are more similar. Time evolution of CR-MHD shocks is qualitatively similar to CR-gasdynamical shocks. However, several potentially interesting differences are apparent. We have run simulations using constant, and nonisotropic, obliquity (and hence spatially) dependent forms of the diffusion coefficient kappa. Comparison of the results shows that while the final steady states achieved are the same in each case, the history of CR-MHD shocks can be strongly modified by variations in kappa and, therefore, in the acceleration timescale. Also, the coupling of CR and MHD in low beta, oblique shocks substantially influences the transient density spike that forms in strongly CR-modified shocks. We find that inside the density spike a MHD slow mode wave can be generated that eventually steepens into a shock. A strong shear layer develops within the density spike, driven by MHD stresses. We conjecture that currents in the shear layer could, in nonplanar flows, result in enhanced particle accretion through drift acceleration