Equation of State in Relativistic Magnetohydrodynamics: variable versus constant adiabatic index

The role of the equation of state for a perfectly conducting, relativistic magnetized fluid is the main subject of this work. The ideal constant $\Gamma$-law equation of state, commonly adopted in a wide range of astrophysical applications, is compared with a more realistic equation of state that better approximates the single-specie relativistic gas. The paper focus on three different topics. First, the influence of a more realistic equation of state on the propagation of fast magneto-sonic shocks is investigated. This calls into question the validity of the constant $\Gamma$-law equation of state in problems where the temperature of the gas substantially changes across hydromagnetic waves. Second, we present a new inversion scheme to recover primitive variables (such as rest-mass density and pressure) from conservative ones that allows for a general equation of state and avoids catastrophic numerical cancellations in the non-relativistic and ultrarelativistic limits. Finally, selected numerical tests of astrophysical relevance (including magnetized accretion flows around Kerr black holes) are compared using different equations of state. Our main conclusion is that the choice of a realistic equation of state can considerably bear upon the solution when transitions from cold to hot gas (or viceversa) are present. Under these circumstances, a polytropic equation of state can significantly endanger the solution.Comment: 14 pages, 14 figure

Two-component jet simulations: II. Combining analytical disk and stellar MHD outflow solutions

Theoretical arguments along with observational data of YSO jets suggest the presence of two steady components: a disk wind type outflow needed to explain the observed high mass loss rates and a stellar wind type outflow probably accounting for the observed stellar spin down. Each component's contribution depends on the intrinsic physical properties of the YSO-disk system and its evolutionary stage. The main goal of this paper is to understand some of the basic features of the evolution, interaction and co-existence of the two jet components over a parameter space and when time variability is enforced. Having studied separately the numerical evolution of each type of the complementary disk and stellar analytical wind solutions in Paper I of this series, we proceed here to mix together the two models inside the computational box. The evolution in time is performed with the PLUTO code, investigating the dynamics of the two-component jets, the modifications each solution undergoes and the potential steady state reached.Comment: accepted for publication in A&

TPCI: The PLUTO-CLOUDY Interface

We present an interface between the (magneto-) hydrodynamics code PLUTO and the plasma simulation and spectral synthesis code CLOUDY. By combining these codes, we constructed a new photoionization hydrodynamics solver: The PLUTO-CLOUDY Interface (TPCI), which is well suited to simulate photoevaporative flows under strong irradiation. The code includes the electromagnetic spectrum from X-rays to the radio range and solves the photoionization and chemical network of the 30 lightest elements. TPCI follows an iterative numerical scheme: First, the equilibrium state of the medium is solved for a given radiation field by CLOUDY, resulting in a net radiative heating or cooling. In the second step, the latter influences the (magneto-) hydrodynamic evolution calculated by PLUTO. Here, we validated the one-dimensional version of the code on the basis of four test problems: Photoevaporation of a cool hydrogen cloud, cooling of coronal plasma, formation of a Stroemgren sphere, and the evaporating atmosphere of a hot Jupiter. This combination of an equilibrium photoionization solver with a general MHD code provides an advanced simulation tool applicable to a variety of astrophysical problems.Comment: 13 pages, 10 figures, accepted for publication in A&

The Athena Astrophysical MHD Code in Cylindrical Geometry

A method for implementing cylindrical coordinates in the Athena magnetohydrodynamics (MHD) code is described. The extension follows the approach of Athena's original developers and has been designed to alter the existing Cartesian-coordinates code as minimally and transparently as possible. The numerical equations in cylindrical coordinates are formulated to maintain consistency with constrained transport, a central feature of the Athena algorithm, while making use of previously implemented code modules such as the Riemann solvers. Angular-momentum transport, which is critical in astrophysical disk systems dominated by rotation, is treated carefully. We describe modifications for cylindrical coordinates of the higher-order spatial reconstruction and characteristic evolution steps as well as the finite-volume and constrained transport updates. Finally, we present a test suite of standard and novel problems in one-, two-, and three-dimensions designed to validate our algorithms and implementation and to be of use to other code developers. The code is suitable for use in a wide variety of astrophysical applications and is freely available for download on the web

Young stellar object jet models: From theory to synthetic observations

Astronomical observations, analytical solutions and numerical simulations have provided the building blocks to formulate the current theory of young stellar object jets. Although each approach has made great progress independently, it is only during the last decade that significant efforts are being made to bring the separate pieces together. Building on previous work that combined analytical solutions and numerical simulations, we apply a sophisticated cooling function to incorporate optically thin energy losses in the dynamics. On the one hand, this allows a self-consistent treatment of the jet evolution and on the other, it provides the necessary data to generate synthetic emission maps. Firstly, analytical disk and stellar outflow solutions are properly combined to initialize numerical two-component jet models inside the computational box. Secondly, magneto-hydrodynamical simulations are performed in 2.5D, following properly the ionization and recombination of a maximum of $29$ ions. Finally, the outputs are post-processed to produce artificial observational data. The first two-component jet simulations, based on analytical models, that include ionization and optically thin radiation losses demonstrate promising results for modeling specific young stellar object outflows. The generation of synthetic emission maps provides the link to observations, as well as the necessary feedback for the further improvement of the available models.Comment: accepted for publication A&A, 20 pages, 11 figure

TESS: A Relativistic Hydrodynamics Code on a Moving Voronoi Mesh

We have generalized a method for the numerical solution of hyperbolic systems of equations using a dynamic Voronoi tessellation of the computational domain. The Voronoi tessellation is used to generate moving computational meshes for the solution of multi-dimensional systems of conservation laws in finite-volume form. The mesh generating points are free to move with arbitrary velocity, with the choice of zero velocity resulting in an Eulerian formulation. Moving the points at the local fluid velocity makes the formulation effectively Lagrangian. We have written the TESS code to solve the equations of compressible hydrodynamics and magnetohydrodynamics for both relativistic and non-relativistic fluids on a dynamic Voronoi mesh. When run in Lagrangian mode, TESS is significantly less diffusive than fixed mesh codes and thus preserves contact discontinuities to high precision while also accurately capturing strong shock waves. TESS is written for Cartesian, spherical and cylindrical coordinates and is modular so that auxilliary physics solvers are readily integrated into the TESS framework and so that the TESS framework can be readily adapted to solve general systems of equations. We present results from a series of test problems to demonstrate the performance of TESS and to highlight some of the advantages of the dynamic tessellation method for solving challenging problems in astrophysical fluid dynamics.Comment: ApJS, 197, 1

The Dynamics of Radiative Shock Waves: Linear and Nonlinear Evolution

The stability properties of one-dimensional radiative shocks with a power-law cooling function of the form $\Lambda \propto \rho^2T^\alpha$ are the main subject of this work. The linear analysis originally presented by Chevalier & Imamura, is thoroughfully reviewed for several values of the cooling index $\alpha$ and higher overtone modes. Consistently with previous results, it is shown that the spectrum of the linear operator consists in a series of modes with increasing oscillation frequency. For each mode a critical value of the cooling index, $\alpha_\textrm{c}$, can be defined so that modes with $\alpha < \alpha_\textrm{c}$ are unstable, while modes with $\alpha > \alpha_\textrm{c}$ are stable. The perturbative analysis is complemented by several numerical simulations to follow the time-dependent evolution of the system for different values of $\alpha$. Particular attention is given to the comparison between numerical and analytical results (during the early phases of the evolution) and to the role played by different boundary conditions. It is shown that an appropriate treatment of the lower boundary yields results that closely follow the predicted linear behavior. During the nonlinear regime, the shock oscillations saturate at a finite amplitude and tend to a quasi-periodic cycle. The modes of oscillations during this phase do not necessarily coincide with those predicted by linear theory, but may be accounted for by mode-mode coupling.Comment: 33 pages, 12 figures, accepted for publication on the Astrophysical Journa

Intestinal helminths of red foxes (Vulpes vulpes) in north-west Italy

SummaryA total of 180 foxes (Vulpes vulpes) from an area scarcely investigated of north-west Italy, were examined for intestinal helminths using sedimentation and counting technique (SCT). Faecal samples were submitted to centrifugation with 50 % zinc sulphate used as flotation solution.No fox was found completely negative for intestinal helminths. The most frequently identified nematodes were Uncinaria stenocephala (70.0 %), Molineus legerae (27.2 %), Toxocara canis (26.7 %), Toxascaris leonina (25.6 %), Trichuris vulpis (21.1 %), Aonchotheca putorii (8.9 %), Pterygodermatites affinis (5.6 %). Genus Mesocestoides (81.7 %), family Dilepididae (29.4 %) and Taenia spp. (8.3 %) were the most prevalent cestodes. All foxes were negative for E. multilocularis and E. granulosus. In two foxes trematodes belonging to the family Plagiorchidae were found.The study highlighted that foxes are hosts of intestinal helminths of veterinary and medical importance which may be transmitted to dogs and humans

The PLUTO Code for Adaptive Mesh Computations in Astrophysical Fluid Dynamics

We present a description of the adaptive mesh refinement (AMR) implementation of the PLUTO code for solving the equations of classical and special relativistic magnetohydrodynamics (MHD and RMHD). The current release exploits, in addition to the static grid version of the code, the distributed infrastructure of the CHOMBO library for multidimensional parallel computations over block-structured, adaptively refined grids. We employ a conservative finite-volume approach where primary flow quantities are discretized at the cell-center in a dimensionally unsplit fashion using the Corner Transport Upwind (CTU) method. Time stepping relies on a characteristic tracing step where piecewise parabolic method (PPM), weighted essentially non-oscillatory (WENO) or slope-limited linear interpolation schemes can be handily adopted. A characteristic decomposition-free version of the scheme is also illustrated. The solenoidal condition of the magnetic field is enforced by augmenting the equations with a generalized Lagrange multiplier (GLM) providing propagation and damping of divergence errors through a mixed hyperbolic/parabolic explicit cleaning step. Among the novel features, we describe an extension of the scheme to include non-ideal dissipative processes such as viscosity, resistivity and anisotropic thermal conduction without operator splitting. Finally, we illustrate an efficient treatment of point-local, potentially stiff source terms over hierarchical nested grids by taking advantage of the adaptivity in time. Several multidimensional benchmarks and applications to problems of astrophysical relevance assess the potentiality of the AMR version of PLUTO in resolving flow features separated by large spatial and temporal disparities.Comment: 34 pages, 34 figures, accepted for publication in ApJ