103 research outputs found
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
A fluid-particle hybrid framework for the PLUTO code: applications to non-thermal emission in jets
We present an implementation of a fully parallel hybrid framework for the evolution of Lagrangian particles coupled to a MHD fluid for the PLUTO code. For the applications of interest, particles represent ensembles of electrons whose spectral energy distribution is governed by a kinetic transport equation that takes into account different physical processes such as diffusive shock acceleration, synchrotron emission and adiabatic expansion. An application to model non-thermal emission from extragalactic jets shows the effectiveness and strength of the approach in describing not only the dynamics but also the radiation properties of jets and, in general, of high-energy astrophysical plasma environments
Making Faranoff-Riley I radio sources. I. Numerical hydrodynamic 3D simulations of low-power jets
Context. Extragalactic radio sources have been classified into two classes, Fanaroff-Riley I and II, which differ in morphology and radio power. Strongly emitting sources belong to the edge-brightened FR II class, and weakly emitting sources to the edge-darkened FR I class. The origin of this dichotomy is not yet fully understood. Numerical simulations are successful in generating FR II morphologies, but they fail to reproduce the diffuse structure of FR Is. Aims: By means of hydro-dynamical 3D simulations of supersonic jets, we investigate how the displayed morphologies depend on the jet parameters. Bow shocks and Mach disks at the jet head, which are probably responsible for the hot spots in the FR II sources, disappear for a jet kinetic power ℒkin ≲ 1043 erg s-1. This threshold compares favorably with the luminosity at which the FR I/FR II transition is observed. Methods: The problem is addressed by numerical means carrying out 3D HD simulations of supersonic jets that propagate in a non-homogeneous medium with the ambient temperature that increases with distance from the jet origin, which maintains constant pressure. Results: The jet energy in the lower power sources, instead of being deposited at the terminal shock, is gradually dissipated by the turbulence. The jets spread out while propagating, and they smoothly decelerate while mixing with the ambient medium and produce the plumes characteristic of FR I objects. Conclusions: Three-dimensionality is an essential ingredient to explore the FR I evolution becausethe properties of turbulence in two and three dimensions are very different, since there is no energy cascade to small scales in two dimensions, and two-dimensional simulations with the same parameters lead to FRII-like behavior
Relativistic MHD Simulations of Jets with Toroidal Magnetic Fields
This paper presents an application of the recent relativistic HLLC
approximate Riemann solver by Mignone & Bodo to magnetized flows with vanishing
normal component of the magnetic field.
The numerical scheme is validated in two dimensions by investigating the
propagation of axisymmetric jets with toroidal magnetic fields.
The selected jet models show that the HLLC solver yields sharper resolution
of contact and shear waves and better convergence properties over the
traditional HLL approach.Comment: 12 pages, 5 figure
- …