The Dynamics of Radiative Shock Waves: Linear and Nonlinear Evolution

33 pages, 12 figures, accepted for publication on the Astrophysical Journa

The Piecewise Parabolic Method for Multidimensional Relativistic Fluid Dynamics

We present an extension of the Piecewise Parabolic Method to special relativistic fluid dynamics in multidimensions. The scheme is conservative, dimensionally unsplit, and suitable for a general equation of state. Temporal evolution is second-order accurate and employs characteristic projection operators; spatial interpolation is piece-wise parabolic making the scheme third-order accurate in smooth regions of the flow away from discontinuities. The algorithm is written for a general system of orthogonal curvilinear coordinates and can be used for computations in non-cartesian geometries. A non-linear iterative Riemann solver based on the two-shock approximation is used in flux calculation. In this approximation, an initial discontinuity decays into a set of discontinuous waves only implying that, in particular, rarefaction waves are treated as flow discontinuities. We also present a new and simple equation of state which approximates the exact result for the relativistic perfect gas with high accuracy. The strength of the new method is demonstrated in a series of numerical tests and more complex simulations in one, two and three dimensions

Linear stability analysis of magnetized relativistic rotating jets

We carry out a linear stability analysis of a magnetized relativistic rotating cylindrical jet flow using the approximation of zero thermal pressure. We identify several modes of instability in the jet: Kelvin-Helmholtz, current driven and two kinds of centrifugal-buoyancy modes -- toroidal and poloidal. The Kelvin-Helmholtz mode is found at low magnetization and its growth rate depends very weakly on the pitch parameter of the background magnetic field and on rotation. The current driven mode is found at high magnetization, the values of its growth rate and the wavenumber, corresponding to the maximum growth, increase as we decrease the pitch parameter of the background magnetic field. This mode is stabilized by rotation, especially, at high magnetization. The centrifugal-buoyancy modes, arising due to rotation, tend also to be more stable when magnetization is increased. Overall, relativistic jet flows appear to be more stable with respect to their non-relativistic counterpart.Comment: 15 pages, 15 figures, accepted for pubblication in MNRA

Radiation hydrodynamics integrated in the code PLUTO

The transport of energy through radiation is very important in many astrophysical phenomena. In dynamical problems the time-dependent equations of radiation hydrodynamics have to be solved. We present a newly developed radiation-hydrodynamics module specifically designed for the versatile MHD code PLUTO. The solver is based on the flux-limited diffusion approximation in the two-temperature approach. All equations are solved in the co-moving frame in the frequency independent (grey) approximation. The hydrodynamics is solved by the different Godunov schemes implemented in PLUTO, and for the radiation transport we use a fully implicit scheme. The resulting system of linear equations is solved either using the successive over-relaxation (SOR) method (for testing purposes), or matrix solvers that are available in the PETSc library. We state in detail the methodology and describe several test cases in order to verify the correctness of our implementation. The solver works in standard coordinate systems, such as Cartesian, cylindrical and spherical, and also for non-equidistant grids. We have presented a new radiation-hydrodynamics solver coupled to the MHD-code \PLUTO that is a modern, versatile and efficient new module for treating complex radiation hydrodynamical problems in astrophysics. As test cases, either purely radiative situations, or full radiation-hydrodynamical setups (including radiative shocks and convection in accretion discs) have been studied successfully. The new module scales very well on parallel computers using MPI. For problems in star or planet formation, we have added the possibility of irradiation by a central source.Comment: 13 pages, 11 figures, accepted by Astronomy & Astrophysic

Scalable explicit implementation of anisotropic diffusion with Runge-Kutta-Legendre super-time-stepping

An important ingredient in numerical modelling of high temperature magnetised astrophysical plasmas is the anisotropic transport of heat along magnetic field lines from higher to lower temperatures.Magnetohydrodynamics (MHD) typically involves solving the hyperbolic set of conservation equations along with the induction equation. Incorporating anisotropic thermal conduction requires to also treat parabolic terms arising from the diffusion operator. An explicit treatment of parabolic terms will considerably reduce the simulation time step due to its dependence on the square of the grid resolution ($\Delta x$) for stability. Although an implicit scheme relaxes the constraint on stability, it is difficult to distribute efficiently on a parallel architecture. Treating parabolic terms with accelerated super-time stepping (STS) methods has been discussed in literature but these methods suffer from poor accuracy (first order in time) and also have difficult-to-choose tuneable stability parameters. In this work we highlight a second order (in time) Runge Kutta Legendre (RKL) scheme (first described by Meyer et. al. 2012) that is robust, fast and accurate in treating parabolic terms alongside the hyperbolic conversation laws. We demonstrate its superiority over the first order super time stepping schemes with standard tests and astrophysical applications. We also show that explicit conduction is particularly robust in handling saturated thermal conduction. Parallel scaling of explicit conduction using RKL scheme is demonstrated up to more than $10^4$ processors.Comment: 15 pages, 9 figures, incorporated comments from the referee. This version is now accepted for publication in MNRA

A Particle Module for the PLUTO code: II - Hybrid Framework for Modeling Non-thermal emission from Relativistic Magnetized flows

We describe a new hybrid framework to model non-thermal spectral signatures from highly energetic particles embedded in a large-scale classical or relativistic MHD flow. Our method makes use of \textit{Lagrangian} particles moving through an Eulerian grid where the (relativistic) MHD equations are solved concurrently. Lagrangian particles follow fluid streamlines and represent ensembles of (real) relativistic particles with a finite energy distribution. The spectral distribution of each particle is updated in time by solving the relativistic cosmic ray transport equation based on local fluid conditions. This enables us to account for a number of physical processes, such as adiabatic expansion, synchrotron and inverse Compton emission. An accurate semi-analytically numerical scheme that combines the method of characteristics with a Lagrangian discretization in the energy coordinate is described. In presence of (relativistic) magnetized shocks, a novel approach to consistently model particle energization due to diffusive shock acceleration has been presented. Our approach relies on a refined shock-detection algorithm and updates the particle energy distribution based on the shock compression ratio, magnetic field orientation and amount of (parameterized) turbulence. The evolved distribution from each \textit{Lagrangian} particle is further used to produce observational signatures like emission maps and polarization signals accounting for proper relativistic corrections. We further demonstrate the validity of this hybrid framework using standard numerical benchmarks and evaluate the applicability of such a tool to study high energy emission from extra-galactic jets.Comment: 23 pages, 14 figures, Accepted for publication in The Astrophysical Journa

MHD simulations of three-dimensional Resistive Reconnection in a cylindrical plasma column

Magnetic reconnection is a plasma phenomenon where a topological rearrangement of magnetic field lines with opposite polarity results in dissipation of magnetic energy into heat, kinetic energy and particle acceleration. Such a phenomenon is considered as an efficient mechanism for energy release in laboratory and astrophysical plasmas. An important question is how to make the process fast enough to account for observed explosive energy releases. The classical model for steady state magnetic reconnection predicts reconnection times scaling as $S^{1/2}$ (where $S$ is the Lundquist number) and yields times scales several order of magnitude larger than the observed ones. Earlier two-dimensional MHD simulations showed that for large Lundquist number the reconnection time becomes independent of $S$ ("fast reconnection" regime) due to the presence of the secondary tearing instability that takes place for $S \gtrsim 1 \times 10^4$. We report on our 3D MHD simulations of magnetic reconnection in a magnetically confined cylindrical plasma column under either a pressure balanced or a force-free equilibrium and compare the results with 2D simulations of a circular current sheet. We find that the 3D instabilities acting on these configurations result in a fragmentation of the initial current sheet in small filaments, leading to enhanced dissipation rate that becomes independent of the Lundquist number already at $S \simeq 1\times 10^3$.Comment: 11 pages, 11 figures, accepted for publication in MNRA