A numerical simulation of the evolution and fate of a Fanaroff–Riley type I jet. The case of 3C 31

• Published in: Monthly notices of the Royal Astronomical Society Vol. 382; no. 2; pp. 526 - 542
• Main Authors: Perucho, M., Martí, J. M.
• Format: Journal Article
• Language: English
• Published: Oxford, UK Blackwell Publishing Ltd 01.12.2007; Blackwell Science; Oxford University Press
• Subjects: Astronomy; Astrophysics; Earth, ocean, space; Exact sciences and technology; galaxies: active; galaxies: individual: 3C 31; galaxies: jets; galaxies: nuclei; Luminosity; Radio astronomy; radio continuum: galaxies; Simulation; Stars & galaxies; X ray observation; Jet flow; Radio galaxies; galaxies: jets; Deceleration; Relativistic gas; Digital simulation; galaxies: individual: 3C 31; Pressure gradients; Continuum; Active galaxies; Lifetime; Equations of state; Theoretical model; galaxies: active; Galaxy nuclei; Cosmic radio sources; Density gradient; Subsonic jet; galaxies: nuclei; Age; radio continuum: galaxies
• ISSN: 0035-8711; 1365-2966
• DOI: 10.1111/j.1365-2966.2007.12454.x

The evolution of Fanaroff–Riley type I (FR I) jets has long been studied in the framework of the FR I–FR II dichotomy. The present paradigm consists of the expansion of overpressured jets in the ambient medium and the generation of standing recollimation shocks, followed by mass entrainment from the external medium that decelerates the jets to subsonic speeds. In this paper, we test the present theoretical and observational models via a relativistic numerical simulation of the jets in the radio galaxy 3C 31. We use the parameters derived from the modelling presented by Laing & Bridle as input parameters for the simulation of the evolution of the source, thus assuming that they have not varied over the lifetime of the source. We simulate about 10 per cent of the total lifetime of the jets in 3C 31. Realistic density and pressure gradients for the atmosphere are used. The simulation includes an equation of state for a two-component relativistic gas that allows a separate treatment of leptonic and baryonic matter. We compare our results with the modelling of the observational data of the source. Our results show that the bow shock evolves self-similarly at a quasi-constant speed, with slight deceleration by the end of the simulation, in agreement with recent X-ray observations that show the presence of bow shocks in FR I sources. The jet expands until it becomes underpressured with respect to the ambient medium, and then recollimates. Subsequent oscillations around pressure equilibrium and generation of standing shocks lead to the mass-loading and disruption of the jet flow. We derive an estimate for the minimum age of the source of t > 1 × 108 yr, which may imply continuous activity of 3C 31 since the triggering of its activity. The simulation shows that weak compact steep spectrum sources may be the young counterparts of FR I sources. We conclude that the observed properties of the jets in 3C 31 are basically recovered by the standing shock scenario.