Showing 1–50 of 51 results for author: Martí, J M

Magnetic and thermal acceleration in extragalactic jets: An application to NGC 315

Authors: L. Ricci, M. Perucho, J. López-Miralles, J. M. Martí, B. Boccardi

Abstract: Aims. Relativistic jets launched from active galactic nuclei accelerate up to highly relativistic velocities within a few parsecs to tens of parsecs. The precise way in which this process takes place is still under study. While magnetic acceleration is known to be able to accelerate relativistic outflows, little attention has been paid to the role of thermal acceleration. The latter has been assum… ▽ More Aims. Relativistic jets launched from active galactic nuclei accelerate up to highly relativistic velocities within a few parsecs to tens of parsecs. The precise way in which this process takes place is still under study. While magnetic acceleration is known to be able to accelerate relativistic outflows, little attention has been paid to the role of thermal acceleration. The latter has been assumed to act only on compact regions, very close to the central engine, and to become negligible on parsec scales. However, this holds under the assumption of small internal energies as compared to the magnetic ones, and whether this is true or what happens when we drop this assumption is currently uncertain. Methods. We use a 2D relativistic magnetohydrodynamical code to explore jet acceleration from sub-parsec to parsec scales. As initial conditions for our models, we use observational constraints on jet properties derived by means of very long baseline interferometry observations for a Fanaroff Riley I radio galaxy, NGC\,315. We investigate the parameter space established for this source and perform a number of simulations of magnetically, thermally or kinetically dominated jets at injection, and compare our results with the observed ones. Results. Our simulated jets show that when thermal energy is comparable to or exceeds magnetic energy, thermal acceleration becomes significant at parsec scales. This result has important consequences, potentially extending the acceleration region far beyond the collimation scales, as thermal acceleration can effectively operate within a conically expanding jet. In all the models, we observe acceleration to be driven by expansion, as expected. A number of our models allow us to reproduce the acceleration and opening angles observed in NGC\,315. Finally, our results indicate that disk-launched winds might play an important role in the jet propagation. △ Less

Submitted 11 December, 2023; originally announced December 2023.

Comments: 25 pages, 16 figures. Accepted for publication in A&A

Filamentary structures as the origin of blazar jet radio variability

Authors: Antonio Fuentes, José L. Gómez, José M. Martí, Manel Perucho, Guang-Yao Zhao, Rocco Lico, Andrei P. Lobanov, Gabriele Bruni, Yuri Y. Kovalev, Andrew Chael, Kazunori Akiyama, Katherine L. Bouman, He Sun, Ilje Cho, Efthalia Traianou, Teresa Toscano, Rohan Dahale, Marianna Foschi, Leonid I. Gurvits, Svetlana Jorstad, Jae-Young Kim, Alan P. Marscher, Yosuke Mizuno, Eduardo Ros, Tuomas Savolainen

Abstract: Supermassive black holes at the centre of active galactic nuclei power some of the most luminous objects in the Universe. Typically, very long baseline interferometric (VLBI) observations of blazars have revealed only funnel-like morphologies with little information of the ejected plasma internal structure, or lacked the sufficient dynamic range to reconstruct the extended jet emission. Here we sh… ▽ More Supermassive black holes at the centre of active galactic nuclei power some of the most luminous objects in the Universe. Typically, very long baseline interferometric (VLBI) observations of blazars have revealed only funnel-like morphologies with little information of the ejected plasma internal structure, or lacked the sufficient dynamic range to reconstruct the extended jet emission. Here we show microarcsecond-scale angular resolution images of the blazar 3C 279 obtained at 22 GHz with the space VLBI mission RadioAstron, which allowed us to resolve the jet transversely and reveal several filaments produced by plasma instabilities in a kinetically dominated flow. Our high angular resolution and dynamic range image suggests that emission features traveling down the jet may manifest as a result of differential Doppler-boosting within the filaments, as opposed to the standard shock-in-jet model invoked to explain blazar jet radio variability. Moreover, we infer that the filaments in 3C 279 are possibly threaded by a helical magnetic field rotating clockwise, as seen in the direction of the flow motion, with an intrinsic helix pitch angle of ~45 degrees in a jet with a Lorentz factor of ~13 at the time of observation. △ Less

Submitted 3 November, 2023; originally announced November 2023.

Comments: 15 pages, 7 figures. Initial version of an article published in Nature Astronomy

Journal ref: Nature Astronomy 7 (2023) 1359

On the application of Jacobian-free Riemann solvers for relativistic radiation magnetohydrodynamics under M1 closure

Authors: Jose López-Miralles, Jose María Martí, Manel Perucho

Abstract: Radiative transfer plays a major role in high-energy astrophysics. In multiple scenarios and in a broad range of energy scales, the coupling between matter and radiation is essential to understand the interplay between theory, observations and numerical simulations. In this paper, we present a novel scheme for solving the equations of radiation relativistic magnetohydrodynamics within the parallel… ▽ More Radiative transfer plays a major role in high-energy astrophysics. In multiple scenarios and in a broad range of energy scales, the coupling between matter and radiation is essential to understand the interplay between theory, observations and numerical simulations. In this paper, we present a novel scheme for solving the equations of radiation relativistic magnetohydrodynamics within the parallel code Lóstrego. These equations, which are formulated taking successive moments of the Boltzmann radiative transfer equation, are solved under the gray-body approximation and the M1 closure using an IMEX time integration scheme. The main novelty of our scheme is that we introduce for the first time in the context of radiation magnetohydrodynamics a family of Jacobian-free Riemann solvers based on internal approximations to the Polynomial Viscosity Matrix, which were demonstrated to be robust and accurate for non-radiative applications. The robustness and the limitations of the new algorithms are tested by solving a collection of one-dimensional and multi-dimensional test problems, both in the free-streaming and in the diffusion radiation transport limits. Due to its stable performance, the applicability of the scheme presented in this paper to real astrophysical scenarios in high-energy astrophysics is promising. In future simulations, we expect to be able to explore the dynamical relevance of photon-matter interactions in the context of relativistic jets and accretion discs, from microquasars and AGN to gamma-ray bursts. △ Less

Submitted 1 December, 2022; originally announced December 2022.

Comments: 21 pages, 13 figures. Accepted for publication in Computer Physics Communications

3D RMHD simulations of jet-wind interactions in High Mass X-ray Binaries

Authors: J. López-Miralles, M. Perucho, J. M. Martí, S. Migliari, V. Bosch-Ramon

Abstract: The interaction of jets in High-Mass X-ray Binaries (HMXBs) with the strong winds driven by the hot companion star in the vicinity of the compact object is fundamental to understand the jet dynamics, non-thermal emission and long-term stability. However, the role of the jet magnetic field in this process is unclear. We study the dynamical role of weak and moderate-to-strong toroidal magnetic field… ▽ More The interaction of jets in High-Mass X-ray Binaries (HMXBs) with the strong winds driven by the hot companion star in the vicinity of the compact object is fundamental to understand the jet dynamics, non-thermal emission and long-term stability. However, the role of the jet magnetic field in this process is unclear. We study the dynamical role of weak and moderate-to-strong toroidal magnetic fields during the first hundreds of seconds of jet propagation, focusing on the magnetized flow dynamics and the mechanisms of energy conversion. We have developed the code Lóstrego v1.0, a new 3D RMHD code to simulate astrophysical plasmas in Cartesian coordinates. Using this tool, we performed the first 3D RMHD numerical simulations of relativistic magnetized jets propagating through the clumpy stellar wind in a HMXB. The overall morphology and dynamics of weakly magnetized jet models is similar to previous hydrodynamical simulations, where the jet head generates a strong shock in the ambient medium and the initial over-pressure with respect to the stellar wind drives one or more recollimation shocks. In the time scales of our simulations, these jets are ballistic and seem to be more stable against internal instabilities than jets with the same power in the absence of fields. However, moderate-to-strong toroidal magnetic fields favour the development of current-driven instabilities and the disruption of the jet within the binary. A detailed analysis of the energy distribution in the relativistic outflow and the ambient medium reveals that both magnetic and internal energies can contribute to the effective acceleration of the jet. We certify that the jet feedback into the ambient medium is highly dependent on the jet energy distribution at injection, where hotter, more dilute and/or more magnetized jets are more efficient, as anticipated by feedback studies in the case of jets in active galaxies. △ Less

Submitted 22 February, 2022; originally announced February 2022.

Comments: 28 pages, 23 figures. Abstract shortened for arXiv. Accepted for publication in A&A

Journal ref: A&A 661, A117 (2022)

Reading M87's DNA: A Double Helix revealing a large scale Helical Magnetic Field

Authors: Alice Pasetto, Carlos Carrasco-Gonzalez, Jose L. Gomez, Jose M. Marti, Manel Perucho, Shane P. O'Sullivan, Craig Anderson, Daniel Jacobo Diaz-Gonzalez, Antonio Fuentes, John Wardle

Abstract: We present unprecedented high fidelity radio images of the M87 jet. We analyzed Jansky Very Large Array (VLA) broadband, full polarization, radio data from 4 to 18 GHz. The observations were taken with the most extended configuration (A configuration), which allow the study of the emission of the jet up to kpc scales with a linear resolution $\sim$10 pc. The high sensitivity and resolution of our… ▽ More We present unprecedented high fidelity radio images of the M87 jet. We analyzed Jansky Very Large Array (VLA) broadband, full polarization, radio data from 4 to 18 GHz. The observations were taken with the most extended configuration (A configuration), which allow the study of the emission of the jet up to kpc scales with a linear resolution $\sim$10 pc. The high sensitivity and resolution of our data allow to resolve the jet width. We confirm a double-helix morphology of the jet material between $\sim$300 pc and $\sim$1 kpc. We found a gradient of the polarization degree with a minimum at the projected axis and maxima at the jet edges, and a gradient in the Faraday depth with opposite signs at the jet edges. We also found that the behavior of the polarization properties along the wide range of frequencies is consistent with internal Faraday depolarization. All these characteristics strongly support the presence of a helical magnetic field in the M87 jet up to 1 kpc from the central black hole although the jet is most likely particle dominated at these large scales. Therefore, we propose a plausible scenario in which the helical configuration of the magnetic field has been maintained to large scales thanks to the presence of Kelvin-Helmholtz instabilities. △ Less

Submitted 13 December, 2021; originally announced December 2021.

Comments: Published on ApJL

Long-term FRII jet evolution in dense environments

Authors: Manel Perucho, José María Martí, Vicent Quilis

Abstract: We present long-term numerical three-dimensional simulations of a relativistic outflow propagating through a galactic ambient medium and environment, up to distances $\sim 100$~kpc. Our aim is to study the role of dense media in the global dynamics of the radio source. We use a relativistic gas equation of state, and a basic description of thermal cooling terms. In previous work, we showed that a… ▽ More We present long-term numerical three-dimensional simulations of a relativistic outflow propagating through a galactic ambient medium and environment, up to distances $\sim 100$~kpc. Our aim is to study the role of dense media in the global dynamics of the radio source. We use a relativistic gas equation of state, and a basic description of thermal cooling terms. In previous work, we showed that a linear perturbation could enhance the jet propagation during the early phases of evolution, by introducing obliquity to the jet reverse shock. Here, we show that this effect is reduced in denser media. We find that the dentist-drill effect acts earlier, due to slower jet propagation and an increased growth of the helical instability. The global morphology of the jet is less elongated, with more prominent lobes. The fundamental physical parameters of the jet generated structure derived from our simulations fall within the estimated values derived for FRII jets in the 3C sample. In agreement with previous axisymmetric and three dimensional simulations in lower density media, we conclude that shock heating of the interstellar and intergalactic media is very efficient in the case of powerful, relativistic jets. △ Less

Submitted 7 December, 2021; v1 submitted 6 December, 2021; originally announced December 2021.

Comments: Accepted for publication in Monthly Notices of the Royal Astronomical Society

On the deceleration of Fanaroff-Riley Class I jets: mass loading of magnetized jets by stellar winds

Authors: Andreu Anglés-Castillo, Manel Perucho, José María Martí, Robert A. Laing

Abstract: In this paper we present steady-state RMHD simulations that include a mass-load term to study the process of jet deceleration. The mass-load mimics the injection of a proton-electron plasma from stellar winds within the host galaxy into initially pair plasma jets, with mean stellar mass-losses ranging from $10^{-14}$ to $10^{-9}\,{M_\odot\,yr^{-1}}$. The spatial jet evolution covers… ▽ More In this paper we present steady-state RMHD simulations that include a mass-load term to study the process of jet deceleration. The mass-load mimics the injection of a proton-electron plasma from stellar winds within the host galaxy into initially pair plasma jets, with mean stellar mass-losses ranging from $10^{-14}$ to $10^{-9}\,{M_\odot\,yr^{-1}}$. The spatial jet evolution covers $\sim 500\,{\rm pc}$ from jet injection in the grid at 10~pc from the jet nozzle. Our simulations use a relativistic gas equation of state and a pressure profile for the ambient medium. We compare these simulations with previous dynamical simulations of relativistic, non-magnetised jets. Our results show that toroidal magnetic fields can prevent fast jet expansion and the subsequent embedding of further stars via magnetic tension. In this sense, magnetic fields avoid a runaway deceleration process. Furthermore, when the mass-load is large enough to increase the jet density and produce fast, differential jet expansion, the conversion of magnetic energy flux into kinetic energy flux (i.e., magnetic acceleration), helps to delay the deceleration process with respect to non-magnetised jets. We conclude that the typical stellar population in elliptical galaxies cannot explain jet deceleration in classical FRI radio galaxies. However, we observe a significant change in the jet composition, thermodynamical parameters and energy dissipation along its evolution, even for moderate values of the mass-load. △ Less

Submitted 20 October, 2020; originally announced October 2020.

Comments: Accepted for publication in Monthly Notices of the Royal Astronomical Society

Cosmic magnetic fields with MASCLET: an application to galaxy clusters

Authors: Vicent Quilis, José Ma. Martí, Susana Planelles

Abstract: We describe and test a new version of the adaptive mesh refinement (AMR) cosmological code MASCLET. The new version of the code includes all the ingredients of its previous version plus a description of the evolution of the magnetic field under the approximation of the ideal magneto-hydrodynamics (MHD). To preserve the divergence-free condition of MHD, the original divergence cleaning algorithm of… ▽ More We describe and test a new version of the adaptive mesh refinement (AMR) cosmological code MASCLET. The new version of the code includes all the ingredients of its previous version plus a description of the evolution of the magnetic field under the approximation of the ideal magneto-hydrodynamics (MHD). To preserve the divergence-free condition of MHD, the original divergence cleaning algorithm of Dedner et al. (2002) is implemented. We present a set of well-known 1D and 2D tests, such as several shock-tube problems, the fast rotor and the Orszag-Tang vortex. The performance of the code in all the tests is excellent with estimated median relative errors of $\nabla \cdot {\bf B}$ in the 2D tests smaller than $5 \times 10^{-5}$ for the fast rotor test, and $5 \times 10^{-3}$ for the Orszag-Tang vortex. As an astrophysical application of the code, we present a simulation of a cosmological box of 40 comoving Mpc side length in which a primordial uniform comoving magnetic field of strength 0.1 nG is seeded. The simulation shows how the magnetic field is channelled along the filaments of gas and is concentrated and amplified within galaxy clusters. Comparison with the values expected from pure compression reveals an additional amplification of the magnetic field caused by turbulence in the central region of the cluster. Values of the order of $\sim 1 μ$G are obtained in clusters at $z\sim 0$ with median relative errors of $\nabla \cdot {\bf B}$ below 0.4\%. The implications of a proper description of the dynamics of the magnetic field and their possible observational counterparts in future facilities are discussed. △ Less

Submitted 27 March, 2020; originally announced March 2020.

Comments: 14 pages, 13 figures, accepted for publication in MNRAS

Total and linearly polarized synchrotron emission from overpressured magnetized relativistic jets

Authors: Antonio Fuentes, José L. Gómez, José M. Martí, Manel Perucho

Abstract: We present relativistic magnetohydrodynamic (RMHD) simulations of stationary overpressured magnetized relativistic jets which are characterized by their dominant type of energy, namely internal, kinetic, or magnetic. Each model is threaded by a helical magnetic field with a pitch angle of $45^\circ$ and features a series of recollimation shocks produced by the initial pressure mismatch, whose stre… ▽ More We present relativistic magnetohydrodynamic (RMHD) simulations of stationary overpressured magnetized relativistic jets which are characterized by their dominant type of energy, namely internal, kinetic, or magnetic. Each model is threaded by a helical magnetic field with a pitch angle of $45^\circ$ and features a series of recollimation shocks produced by the initial pressure mismatch, whose strength and number varies as a function of the dominant type of energy. We perform a study of the polarization signatures from these models by integrating the radiative transfer equations for synchrotron radiation using as inputs the RMHD solutions. These simulations show a top-down emission asymmetry produced by the helical magnetic field and a progressive confinement of the emission into a jet spine as the magnetization increases and the internal energy of the non-thermal population is considered to be a constant fraction of the thermal one. Bright stationary components associated with the recollimation shocks appear presenting a relative intensity modulated by the Doppler boosting ratio between the pre-shock and post-shock states. Small viewing angles show a roughly bimodal distribution in the polarization angle due to the helical structure of the magnetic field, which is also responsible for the highly stratified degree of linear polarization across the jet width. In addition, small variations of the order of $26^\circ$ are observed in the polarization angle of the stationary components, which can be used to identify recollimation shocks in astrophysical jets. △ Less

Submitted 5 April, 2018; originally announced April 2018.

Comments: Accepted for publication in ApJ

Magnetized, Relativistic Jets

Authors: Manel Perucho, José María Martí, José Luís Gómez, Antonio Fuentes

Abstract: Extragalactic relativistic jets are composed by charged particles and magnetic fields, as inferred from the synchrotron emission that we receive from them. The Larmor radii of the particles propagating along the magnetic field are much smaller than the scales of the problem, providing the necessary coherence to the system to treat is as a flow. We can thus study them using relativistic magnetohydr… ▽ More Extragalactic relativistic jets are composed by charged particles and magnetic fields, as inferred from the synchrotron emission that we receive from them. The Larmor radii of the particles propagating along the magnetic field are much smaller than the scales of the problem, providing the necessary coherence to the system to treat is as a flow. We can thus study them using relativistic magnetohydrodynamics. As a first step, we have studied the structure of steady-state configurations of jets by using numerical simulations. We have used a helical field configuration and have changed different relevant parameters that control the way in which the energy flux is distributed in jets (namely, the proportion of the energy flux carried by internal, kinetic or magnetic energy). Our results show significant differences among the different kinds of jets. Finally, we also report on results based on synthetic maps of our simulated jets. △ Less

Submitted 4 August, 2017; originally announced August 2017.

Comments: Submitted for publication in Proceedings of Science, as contribution to the proceedings of the XII Multifrequency Behaviour of High Energy Cosmic Sources Workshop, held in Palermo, 12-17 June 2017

Radio mode feedback: Does relativity matter?

Authors: Manel Perucho, José María Martí, Vicent Quilis, Marina Borja-Lloret

Abstract: Radio mode feedback, associated with the propagation of powerful outflows in active galaxies, is a crucial ingredient in galaxy evolution. Extragalactic jets are well collimated and relativistic, both in terms of thermodynamics and kinematics. They generate strong shocks in the ambient medium, associated with observed hotspots, and carve cavities that are filled with the shocked jet flow. In this… ▽ More Radio mode feedback, associated with the propagation of powerful outflows in active galaxies, is a crucial ingredient in galaxy evolution. Extragalactic jets are well collimated and relativistic, both in terms of thermodynamics and kinematics. They generate strong shocks in the ambient medium, associated with observed hotspots, and carve cavities that are filled with the shocked jet flow. In this Letter, we compare the pressure evolution in the hotspot and the cavity generated by relativistic and classical jets. Our results show that the classical approach underestimates the cavity pressure by a factor larger or equal to 2 for a given shocked volume during the whole active phase. The tension between both approaches can only be alleviated by unrealistic jet flow densities or gigantic jet areas in the classical case. As a consequence, the efficiency of a relativistic jet heating the ambient is typically around 20% larger compared with a classical jet, and the heated volume is 2 to 10 times larger during the time evolution. This conflict translates into two substantially disparate manners, both spatially and temporal, of heating the ambient medium. These differences are expected to have relevant implications on the star formation rates of the host galaxies and their evolution. △ Less

Submitted 18 July, 2017; originally announced July 2017.

Comments: Accepted for publication in MNRAS Letters

The internal structure of overpressured, magnetized, relativistic jets

Authors: José M. Martí, Manel Perucho, José L. Gómez

Abstract: This work presents the first characterization of the internal structure of overpressured steady superfast magnetosonic relativistic jets in connection with their dominant type of energy. To this aim, relativistic magnetohydrodynamic simulations of different jet models threaded by a helical magnetic field have been analyzed covering a wide region in the magnetosonic Mach number - specific internal… ▽ More This work presents the first characterization of the internal structure of overpressured steady superfast magnetosonic relativistic jets in connection with their dominant type of energy. To this aim, relativistic magnetohydrodynamic simulations of different jet models threaded by a helical magnetic field have been analyzed covering a wide region in the magnetosonic Mach number - specific internal energy plane. The merit of this plane is that models dominated by different types of energy (internal energy: hot jets; rest-mass energy: kinetically dominated jets; magnetic energy: Poynting-flux dominated jets) occupy well separated regions. The analyzed models also cover a wide range of magnetizations. Models dominated by the internal energy (i.e., hot models, or Poynting-flux dominated jets with magnetizations larger than but close to 1) have a rich internal structure characterized by a series of recollimation shocks and present the largest variations in the flow Lorentz factor (and internal energy density). Conversely, in kinetically dominated models there is not much internal nor magnetic energy to be converted into kinetic one and the jets are featureless, with small variations in the flow Lorentz factor. The presence of a significant toroidal magnetic field threading the jet produces large gradients in the transversal profile of the internal energy density. Poynting-flux dominated models with high magnetization (>10) are prone to be unstable against magnetic pinch modes, which sets limits to the expected magnetization in parsec-scale AGN jets {and/or constrains their magnetic field configuration}. △ Less

Submitted 12 September, 2016; v1 submitted 2 September, 2016; originally announced September 2016.

Comments: 22 pages, 16 figures, changed title, accepted for publication in The Astrophysical Journal

The Spanish Square Kilometre Array White Book

Authors: J. A. Acosta-Pulido, I. Agudo, A. Alberdi, J. Alcolea, E. J. Alfaro, A. Alonso-Herrero, G. Anglada, P. Arnalte-Mur, Y. Ascasibar, B. Ascaso, R. Azulay, R. Bachiller, A. Baez-Rubio, E. Battaner, J. Blasco, C. B. Brook, V. Bujarrabal, G. Busquet, M. D. Caballero-Garcia, C. Carrasco-Gonzalez, J. Casares, A. J. Castro-Tirado, L. Colina, F. Colomer, I. de Gregorio-Monsalvo , et al. (94 additional authors not shown)

Abstract: The Square Kilometre Array (SKA) is called to revolutionise essentially all areas of Astrophysics. With a collecting area of about a square kilometre, the SKA will be a transformational instrument, and its scientific potential will go beyond the interests of astronomers. Its technological challenges and huge cost requires a multinational effort, and Europe has recognised this by putting the SKA on… ▽ More The Square Kilometre Array (SKA) is called to revolutionise essentially all areas of Astrophysics. With a collecting area of about a square kilometre, the SKA will be a transformational instrument, and its scientific potential will go beyond the interests of astronomers. Its technological challenges and huge cost requires a multinational effort, and Europe has recognised this by putting the SKA on the roadmap of the European Strategy Forum for Research Infrastructures (ESFRI). The Spanish SKA White Book is the result of the coordinated effort of 120 astronomers from 40 different research centers. The book shows the enormous scientific interest of the Spanish astronomical community in the SKA and warrants an optimum scientific exploitation of the SKA by Spanish researchers, if Spain enters the SKA project. △ Less

Submitted 17 June, 2015; v1 submitted 10 June, 2015; originally announced June 2015.

Comments: Editors: Miguel A. Perez-Torres (Editor-in-chief), L. Verdes-Montenegro, J. C. Guirado, A. Alberdi, J. Martin-Pintado, R. Bachiller, D. Herranz, J. M. Girart, J. Gorgas, C. Hernandez-Monteagudo, S. Migliari and J. M. Rodriguez Espinosa; 289 pages, published in June 2015 by Sociedad Española de Astronomía

On the correction of conserved variables for numerical RMHD with staggered constrained transport

Authors: José M. Martí

Abstract: Despite the success of the combination of conservative schemes and staggered constrained transport algorithms in the last fifteen years, the accurate description of highly magnetized, relativistic flows with strong shocks represents still a challenge in numerical RMHD. The present paper focusses in the accuracy and robustness of several correction algorithms for the conserved variables, which has… ▽ More Despite the success of the combination of conservative schemes and staggered constrained transport algorithms in the last fifteen years, the accurate description of highly magnetized, relativistic flows with strong shocks represents still a challenge in numerical RMHD. The present paper focusses in the accuracy and robustness of several correction algorithms for the conserved variables, which has become a crucial ingredient in the numerical simulation of problems where the magnetic pressure dominates over the thermal pressure by more than two orders of magnitude. Two versions of non-relativistic and fully relativistic corrections have been tested and compared using a magnetized cylindrical explosion with high magnetization ($ \ge 10^4$) as test. In the non-relativistic corrections, the total energy is corrected for the difference in the classical magnetic energy term between the average of the staggered fields and the conservative ones, before (CA1) and after (CA1') recovering the primitive variables. These corrections are unable to pass the test at any numerical resolution. The two relativistic approaches (CA2 and CA2'), correcting also the magnetic terms depending on the flow speed in both the momentum and the total energy, reveal as much more robust. These algorithms pass the test succesfully and with very small deviations of the energy conservation ($\le 10^{-4}$), and very low values of the total momentum ($\le 10^{-8}$). In particular, the algorithm CA2' (that corrects the conserved variables after recovering the primitive variables) passes the test at all resolutions. The numerical code used to run all the test cases is briefly described. △ Less

Submitted 12 January, 2015; originally announced January 2015.

Comments: 17 pages, 2 figures, accepted for publication in Computer Physics Communications

Large-scale jets from active galactic nuclei as a source of ICM heating: cavities and shocks

Authors: Manel Perucho, José María Martí, Vicent Quilis, Elena Ricciardelli

Abstract: The evolution of powerful extragalactic jets is not only interesting by itself, but also for its impact on the evolution of the host galaxy and its surroundings. We have performed long-term axisymmetric numerical simulations of relativistic jets with different powers to study their evolution through an environment with a pressure and density gradient. Our results show key differences in the evolut… ▽ More The evolution of powerful extragalactic jets is not only interesting by itself, but also for its impact on the evolution of the host galaxy and its surroundings. We have performed long-term axisymmetric numerical simulations of relativistic jets with different powers to study their evolution through an environment with a pressure and density gradient. Our results show key differences in the evolution of jets with different powers in terms of the spatial and temporal scales of energy deposition. According to our results, the observed morphology in X-ray cavities requires that an important fraction of the jet's energetic budget is in the form of internal energy. Thus, light, lepton-dominated jets are favoured. In all cases, heating is mainly produced by shocks. Cavity overpressure is sustained by an important population of thermal particles. Our simulations reproduce the cool-core structure in projected, luminosity-weighted temperature. We have performed an additional simulation of a slow, massive jet and discuss the differences with its relativistic counterparts. Important qualitative and quantitative differences are found between the non-relativistic and the relativistic jets. Our conclusions point towards a dual-mode of AGN kinetic feedback, depending on the jet power. △ Less

Submitted 11 September, 2014; originally announced September 2014.

Comments: Accepted for publication in MNRAS. The quality of some figures has been reduced to adjust their size to arxiv requirements

On the deceleration of FRI jets: mass loading by stellar-winds

Authors: Manel Perucho, José María Martí, Robert A. Laing, Philip E. Hardee

Abstract: Jets in low-luminosity radio galaxies are known to decelerate from relativistic speeds on parsec scales to mildly or sub-relativistic speeds on kiloparsec scales. Several mechanisms have been proposed to explain this effect, including strong reconfinement shocks and the growth of instabilities (both leading to boundary-layer entrainment) and mass loading from stellar winds or molecular clouds. We… ▽ More Jets in low-luminosity radio galaxies are known to decelerate from relativistic speeds on parsec scales to mildly or sub-relativistic speeds on kiloparsec scales. Several mechanisms have been proposed to explain this effect, including strong reconfinement shocks and the growth of instabilities (both leading to boundary-layer entrainment) and mass loading from stellar winds or molecular clouds. We have performed a series of axisymmetric simulations of the early evolution of jets in a realistic ambient medium to probe the effects of mass loading from stellar winds using the code Ratpenat. We study the evolution of Fanaroff-Riley Class I (FRI) jets, with kinetic powers L_j \sim 1.e41-1.e44 erg/s, within the first 2 kpc of their evolution, where deceleration by stellar mass loading should be most effective. Mass entrainment rates consistent with present models of stellar mass loss in elliptical galaxies produce deceleration and effective decollimation of weak FRI jets within the first kiloparsec. However, powerful FRI jets are not decelerated significantly. In those cases where the mass loading is important, the jets show larger opening angles and decollimate at smaller distances, but the overall structure and dynamics of the bow-shock are similar to those of unloaded jets with the same power and thrust. According to our results, the flaring observed on kpc scales is initiated by mass loading in the weaker FRI jets and by reconfinement shocks or the growth of instabilities in the more powerful jets. The final mechanism of decollimation and deceleration is always the development of disruptive pinching modes. △ Less

Submitted 4 April, 2014; originally announced April 2014.

Comments: Accepted for publication in Monthly Notices of the Royal Astronomical Society

Intracluster Medium reheating by relativistic jets

Authors: Manel Perucho, Vicent Quilis, José María Martí

Abstract: Galactic jets are powerful energy sources reheating the intra-cluster medium in galaxy clusters. Their crucial role in the cosmic puzzle, motivated by observations, has been established by a great number of numerical simulations missing the relativistic nature of these jets. We present the first relativistic simulations of the very long term evolution of realistic galactic jets. Unexpectedly, our… ▽ More Galactic jets are powerful energy sources reheating the intra-cluster medium in galaxy clusters. Their crucial role in the cosmic puzzle, motivated by observations, has been established by a great number of numerical simulations missing the relativistic nature of these jets. We present the first relativistic simulations of the very long term evolution of realistic galactic jets. Unexpectedly, our results show no buoyant bubbles, but large cocoon regions compatible with the observed X-ray cavities. The reheating is more efficient and faster than in previous scenarios, and it is produced by the shock wave driven by the jet, that survives for several hundreds of Myrs. Therefore, the X-ray cavities in clusters produced by powerful relativistic jets would remain confined by weak shocks for extremely long periods, whose detection could be an observational challenge. △ Less

Submitted 13 October, 2011; v1 submitted 27 July, 2011; originally announced July 2011.

Comments: Accepted for publication in ApJ

A Roe-type Riemann solver based on the spectral decomposition of the equations of Relativistic Magnetohydrodynamics

Authors: J. M. Ibáñez, M. A. Aloy, P. Mimica, L. Antón, J. A. Miralles, J. M. Martí

Abstract: In a recent paper (Antón et al. 2010) we have derived sets of right and left eigenvectors of the Jacobians of the relativistic MHD equations, which are regular and span a complete basis in any physical state including degenerate ones. We present a summary of the main steps followed in the above derivation and the numerical experiments carried out with the linearized (Roe-type) Riemann solver we ha… ▽ More In a recent paper (Antón et al. 2010) we have derived sets of right and left eigenvectors of the Jacobians of the relativistic MHD equations, which are regular and span a complete basis in any physical state including degenerate ones. We present a summary of the main steps followed in the above derivation and the numerical experiments carried out with the linearized (Roe-type) Riemann solver we have developed, and some note on the (non-)convex character of the relativistic MHD equations. △ Less

Submitted 11 December, 2010; originally announced December 2010.

Comments: 6 pages, Proceedings of the ASTRONUM-2010 meeting held at San Diego (USA). To by published in ASP Conf. Ser

Stability of three-dimensional relativistic jets: implications for jet collimation

Authors: M. Perucho, J. M. Marti, J. M. Cela, M. Hanasz, R. de la Cruz, F. Rubio

Abstract: The stable propagation of jets in FRII sources is remarkable if one takes into account that large-scale jets are subjected to potentially highly disruptive three-dimensional (3D) Kelvin-Helmholtz instabilities. Numerical simulations can address this problem and help clarify the causes of this remarkable stability. Following previous studies of the stability of relativistic flows in two dimensions… ▽ More The stable propagation of jets in FRII sources is remarkable if one takes into account that large-scale jets are subjected to potentially highly disruptive three-dimensional (3D) Kelvin-Helmholtz instabilities. Numerical simulations can address this problem and help clarify the causes of this remarkable stability. Following previous studies of the stability of relativistic flows in two dimensions (2D), it is our aim to test and extend the conclusions of such works to three dimensions. We present numerical simulations for the study of the stability properties of 3D, sheared, relativistic flows. This work uses a fully parallelized code Ratpenat that solves equations of relativistic hydrodynamics in 3D. The results of the present simulations confirm those in 2D. We conclude that the growth of resonant modes in sheared relativistic flows could be important in explaining the long-term collimation of extragalactic jets. △ Less

Submitted 24 May, 2010; originally announced May 2010.

Comments: Accepted for publication in A&A

Relativistic Magnetohydrodynamics: Renormalized eigenvectors and full wave decomposition Riemann solver

Authors: L. Anton, J. A. Miralles, J. M. Marti, J. M. Ibanez, M. A. Aloy, P. Mimica

Abstract: We obtain renormalized sets of right and left eigenvectors of the flux vector Jacobians of the relativistic MHD equations, which are regular and span a complete basis in any physical state including degenerate ones. The renormalization procedure relies on the characterization of the degeneracy types in terms of the normal and tangential components of the magnetic field to the wavefront in the fl… ▽ More We obtain renormalized sets of right and left eigenvectors of the flux vector Jacobians of the relativistic MHD equations, which are regular and span a complete basis in any physical state including degenerate ones. The renormalization procedure relies on the characterization of the degeneracy types in terms of the normal and tangential components of the magnetic field to the wavefront in the fluid rest frame. Proper expressions of the renormalized eigenvectors in conserved variables are obtained through the corresponding matrix transformations. Our work completes previous analysis that present different sets of right eigenvectors for non-degenerate and degenerate states, and can be seen as a relativistic generalization of earlier work performed in classical MHD. Based on the full wave decomposition (FWD) provided by the the renormalized set of eigenvectors in conserved variables, we have also developed a linearized (Roe-type) Riemann solver. Extensive testing against one- and two-dimensional standard numerical problems allows us to conclude that our solver is very robust. When compared with a family of simpler solvers that avoid the knowledge of the full characteristic structure of the equations in the computation of the numerical fluxes, our solver turns out to be less diffusive than HLL and HLLC, and comparable in accuracy to the HLLD solver. The amount of operations needed by the FWD solver makes it less efficient computationally than those of the HLL family in one-dimensional problems. However its relative efficiency increases in multidimensional simulations. △ Less

Submitted 23 December, 2009; originally announced December 2009.

Comments: 50 pages, 17 figures (2 in color). Submitted to ApJ Suppl. Ser

Journal ref: Astrophys.J.Suppl.188:1-31,2010

Spectral evolution of superluminal components in parsec-scale jets

Authors: P. Mimica, M. A. Aloy, I. Agudo, J. M. Marti, J. L. Gomez, J. A. Miralles

Abstract: (Abridged) We present numerical simulations of the spectral evolution and emission of radio components in relativistic jets. We have developed an algorithm (SPEV) for the transport of a population of non-thermal electrons including radiative losses. For large values of the ratio of gas pressure to magnetic field energy density, \ab \sim 6\times 10^4, quiescent jet models show substantial spectra… ▽ More (Abridged) We present numerical simulations of the spectral evolution and emission of radio components in relativistic jets. We have developed an algorithm (SPEV) for the transport of a population of non-thermal electrons including radiative losses. For large values of the ratio of gas pressure to magnetic field energy density, \ab \sim 6\times 10^4, quiescent jet models show substantial spectral evolution, with observational consequences only above radio frequencies. Larger values of the magnetic field (\ab \sim 6\times 10^2), such that synchrotron losses are moderately important at radio frequencies, present a larger ratio of shocked-to-unshocked regions brightness than the models without radiative losses, despite the fact that they correspond to the same underlying hydrodynamic structure. We also show that jets with a positive photon spectral index result if the lower limit γ_min of the non-thermal particle energy distribution is large enough. A temporary increase of the Lorentz factor at the jet inlet produces a traveling perturbation that appears in the synthetic maps as a superluminal component. We show that trailing components can be originated not only in pressure matched jets, but also in over-pressured ones, where the existence of recollimation shocks does not allow for a direct identification of such features as Kelvin-Helmholtz modes, and its observational imprint depends on the observing frequency. If the magnetic field is large (\ab \sim 6\times 10^2), the spectral index in the rarefaction trailing the traveling perturbation does not change much with respect to the same model without any hydrodynamic perturbation. If the synchrotron losses are considered the spectral index displays a smaller value than in the corresponding region of the quiescent jet model. △ Less

Submitted 20 February, 2009; v1 submitted 7 November, 2008; originally announced November 2008.

Comments: Accepted for publication in ApJ. A PDF with high-resolution figures and animations referenced in the article can be found at http://www.uv.es/aloy/SPEV_MOVIES

Journal ref: Astrophys.J.696:1142-1163,2009

The Influence of Helical Magnetic Fields in the Dynamics and Emission of Relativistic Jets

Authors: M. Roca-Sogorb, M. Perucho, J. L. Gomez, J. M. Marti, L. Anton, M. A. Aloy, I. Agudo

Abstract: We present numerical relativistic magnetohydrodynamic and emission simulations aimed to study the role played by the magnetic field in the dynamics and emission of relativistic jets in Active Galactic Nuclei. We focus our analysis on the study of the emission from recollimation shocks since they may provide an interpretation for the stationary components seen at parsec-scales in multiple sources… ▽ More We present numerical relativistic magnetohydrodynamic and emission simulations aimed to study the role played by the magnetic field in the dynamics and emission of relativistic jets in Active Galactic Nuclei. We focus our analysis on the study of the emission from recollimation shocks since they may provide an interpretation for the stationary components seen at parsec-scales in multiple sources. We show that the relative brightness of the knots associated with the recollimation shocks decreases with increasing jet magnetization, suggesting that jets presenting stationary components may have a relatively weak magnetization, with magnetic fields of the order of equipartition or below. △ Less

Submitted 26 February, 2008; originally announced February 2008.

Comments: 4 pages, 3 figures. Proceedings from Astronomical Society of the Pacific Conference Series (ASPCS)

A numerical simulation of the evolution and fate of a FRI jet. The case of 3C 31

Authors: M. Perucho, J. M. Marti

Abstract: The evolution of FRI jets has been long studied in the framework of the FRI-FRII dichotomy. 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 \cite{lb02a,lb02b} as input parameters for the simulation of the evolution of the source, thu… ▽ More The evolution of FRI jets has been long studied in the framework of the FRI-FRII dichotomy. 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 \cite{lb02a,lb02b} 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 % 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 FRI 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. 10^8 \rm{yrs}$, which may imply continuous activity of 3C 31 since the triggering of its activity. The simulation shows that weak CSS sources may be the young counterparts of FRIs. We conclude that the observed properties of the jets in 3C 31 are basically recovered by the standing shock scenario. △ Less

Submitted 12 September, 2007; originally announced September 2007.

Comments: Accepted for publication in MNRAS. For better quality figures, please check http://www.mpifr-bonn.mpg.de/staff/mperucho/Research.html

Analytical and numerical studies of fluid instabilities in relativistic jets

Authors: M. Perucho, A. P. Lobanov, J. M. Marti

Abstract: Relativistic outflows represent one of the best-suited tools to probe the physics of AGN. Numerical modelling of internal structure of the relativistic outflows on parsec scales provides important clues about the conditions and dynamics of the material in the immediate vicinity of the central black holes in AGN. We investigate possible causes of the structural patterns and regularities observed… ▽ More Relativistic outflows represent one of the best-suited tools to probe the physics of AGN. Numerical modelling of internal structure of the relativistic outflows on parsec scales provides important clues about the conditions and dynamics of the material in the immediate vicinity of the central black holes in AGN. We investigate possible causes of the structural patterns and regularities observed in the parsec-scale jet of the well known quasar 3C 273. We compare the model with the radio structure observed in 3C 273 on parsec scales using very long baseline interferometry (VLBI) and constrain the basic properties of the flow. Our results show that Kelvin-Helmholtz instabilities are the most plausible mechanism to generate the observed structures. △ Less

Submitted 8 February, 2006; originally announced February 2006.

Comments: Proceedings of Relativistic Astrophysics & Cosmology (Munchen, November 2005). 3 pages

Nonlinear stability of relativistic sheared planar jets

Authors: M. Perucho, J. M. Marti, M. Hanasz

Abstract: The linear and non-linear stability of sheared, relativistic planar jets is studied by means of linear stability analysis and numerical hydrodynamical simulations. Our results extend the previous Kelvin-Hemlholtz stability studies for relativistic, planar jets in the vortex sheet approximation performed by Perucho et al. (2004a,b) by including a shear layer between the jet and the external mediu… ▽ More The linear and non-linear stability of sheared, relativistic planar jets is studied by means of linear stability analysis and numerical hydrodynamical simulations. Our results extend the previous Kelvin-Hemlholtz stability studies for relativistic, planar jets in the vortex sheet approximation performed by Perucho et al. (2004a,b) by including a shear layer between the jet and the external medium and more general perturbations. The models considered span a wide range of Lorentz factors ($2.5-20$) and internal energies ($0.08 c^2-60 c^2$) and are classified into three classes according to the main characteristics of their long-term, non-linear evolution. We observe a clear separation of these three groups in a relativistic Mach-number Lorentz-factor plane. Jets with a low Lorentz factor and small relativistic Mach number are disrupted after saturation. Those with a large Lorentz factor and large relativistic Mach number are the stablest, due to the appearance of short wavelength resonant modes which generate local mixing and heating in the shear layer around a fast, unmixed core, giving a plausible solution for the problem of the long-term stability of relativistic jets. A third group is present between them, including jets with intermediate values of Lorentz factor and relativistic Mach number, which are disrupted by a slow process of mixing favored by an efficient and continuous conversion of kinetic into internal energy. In the long term, all the models develop a distinct transversal structure (shear/transition layers) as a consequence of KH perturbation growth, depending on the class they belong to. The properties of these shear layers are analyzed in connection with the parameters of the original jet models. △ Less

Submitted 14 October, 2005; originally announced October 2005.

Comments: accepted for publication in A&A (in press). High resolution plots, figures and Appendices of the paper will be found in the online version of the paper in A&A, and on request to perucho@mpifr-bonn.mpg.de

The exact solution of the Riemann problem in relativistic MHD with tangential magnetic fields

Authors: Roberto Romero, Jose M. Marti, Jose A. Pons, Jose M. Ibanez, Juan A. Miralles

Abstract: We have extended the procedure to find the exact solution of the Riemann problem in relativistic hydrodynamics to a particular case of relativistic magnetohydrodynamics in which the magnetic field of the initial states is tangential to the discontinuity and orthogonal to the flow velocity. The wave pattern produced after the break up of the initial discontinuity is analogous to the non--magnetic… ▽ More We have extended the procedure to find the exact solution of the Riemann problem in relativistic hydrodynamics to a particular case of relativistic magnetohydrodynamics in which the magnetic field of the initial states is tangential to the discontinuity and orthogonal to the flow velocity. The wave pattern produced after the break up of the initial discontinuity is analogous to the non--magnetic case and we show that the problem can be understood as a purely relativistic hydrodynamical problem with a modified equation of state. The new degree of freedom introduced by the non-zero component of the magnetic field results in interesting effects consisting in the change of the wave patterns for given initial thermodynamical states, in a similar way to the effects arising from the introduction of tangential velocities. Secondly, when the magnetic field dominates the thermodynamical pressure and energy, the wave speeds approach the speed of light leading to fast shocks and fast and arbitrarily thin rarefaction waves. Our approach is the first non-trivial exact solution of a Riemann problem in relativistic magnetohydrodynamics and it can also be of great interest to test numerical codes against known analytical or exact solutions. △ Less

Submitted 22 June, 2005; originally announced June 2005.

Comments: 17 pages, 5 figures, accepted for publication in the Journal of Fluid Mechanics

Journal ref: J.Fluid Mech.544:323,2005

Numerical 3+1 general relativistic magnetohydrodynamics: a local characteristic approach

Authors: L. Anton, O. Zanotti, J. A. Miralles, J. M. Marti, J. M. Ibanez, J. A. Font, J. A. Pons

Abstract: We present a general procedure to solve numerically the general relativistic magnetohydrodynamics (GRMHD) equations within the framework of the 3+1 formalism. The work reported here extends our previous investigation in general relativistic hydrodynamics (Banyuls et al. 1997) where magnetic fields were not considered. The GRMHD equations are written in conservative form to exploit their hyperbol… ▽ More We present a general procedure to solve numerically the general relativistic magnetohydrodynamics (GRMHD) equations within the framework of the 3+1 formalism. The work reported here extends our previous investigation in general relativistic hydrodynamics (Banyuls et al. 1997) where magnetic fields were not considered. The GRMHD equations are written in conservative form to exploit their hyperbolic character in the solution procedure. All theoretical ingredients necessary to build up high-resolution shock-capturing schemes based on the solution of local Riemann problems (i.e. Godunov-type schemes) are described. In particular, we use a renormalized set of regular eigenvectors of the flux Jacobians of the relativistic magnetohydrodynamics equations. In addition, the paper describes a procedure based on the equivalence principle of general relativity that allows the use of Riemann solvers designed for special relativistic magnetohydrodynamics in GRMHD. Our formulation and numerical methodology are assessed by performing various test simulations recently considered by different authors. These include magnetized shock tubes, spherical accretion onto a Schwarzschild black hole, equatorial accretion onto a Kerr black hole, and magnetized thick accretion disks around a black hole prone to the magnetorotational instability. △ Less

Submitted 2 June, 2005; originally announced June 2005.

Comments: 18 pages, 8 figures, submitted to ApJ

Journal ref: Astrophys.J.637:296-312,2006

Real vs. simulated relativistic jets

Authors: J. L. Gomez, J. M. Marti, I. Agudo, A. P. Marscher, S. G. Jorstad, M. A. Aloy

Abstract: Intensive VLBI monitoring programs of jets in AGN are showing the existence of intricate emission patterns, such as upstream motions or slow moving and quasi-stationary componentes trailing superluminal features. Relativistic hydrodynamic and emission simulations of jets are in very good agreement with these observations, proving as a powerful tool for the understanding of the physical processes… ▽ More Intensive VLBI monitoring programs of jets in AGN are showing the existence of intricate emission patterns, such as upstream motions or slow moving and quasi-stationary componentes trailing superluminal features. Relativistic hydrodynamic and emission simulations of jets are in very good agreement with these observations, proving as a powerful tool for the understanding of the physical processes taking place in the jets of AGN, microquasars and GRBs. These simulations show that the variability of the jet emission is the result of a complex combination of phase motions, viewing angle selection effects, and non-linear interactions between perturbations and the underlying jet and/or ambient medium. Both observations and simulations suggest that shock-in-jet models may be an overly simplistic idealization when interpreting the emission structure observed in actual jets. △ Less

Submitted 3 January, 2005; originally announced January 2005.

Comments: 6 pages. 4 figures. Proceedings of the 7th European VLBI Network Symposium held in Toledo, Spain on October 12-15, 2004. Editors: R. Bachiller, F. Colomer, J.-F. Desmurs, P. de Vicente (Observatorio Astronomico Nacional), p. 93-98. Needs evn2004.cls

Simulations of the relativistic parsec-scale jet in 3C273

Authors: M. Perucho, A. P. Lobanov, J. M. Marti

Abstract: We present a hydrodynamical 3D simulation of the relativistic jet in 3C273, in comparison to previous linear perturbation analysis of Kelvin-Helmholtz instability developing in the jet. Our aim is to assess advantages and limitations of both analytical and numerical approaches and to identify spatial and temporal scales on which the linear regime of Kelvin-Helmholtz instability can be applied in… ▽ More We present a hydrodynamical 3D simulation of the relativistic jet in 3C273, in comparison to previous linear perturbation analysis of Kelvin-Helmholtz instability developing in the jet. Our aim is to assess advantages and limitations of both analytical and numerical approaches and to identify spatial and temporal scales on which the linear regime of Kelvin-Helmholtz instability can be applied in studies of morphology and kinematics of parsec-scale jets. △ Less

Submitted 17 December, 2004; originally announced December 2004.

Comments: 4 pages, 3 figures; to be published in Proceedings of the workshop "Multiband Approach to AGN", held on Sep.30-Oct.2 in Bonn. Publication: Memorie della Societa Astronomica Italiana, v. 26, No.1 (2005). Reduced figure resolution! Version with original figures is availavble at http://www.mpifr-bonn.mpg.de/bonn04/proceedings/perucho.pdf

A dynamical model for the CSO-MSO-FRII evolution: hints from hot spot properties

Authors: M. Perucho, J. M. Marti

Abstract: Compact Symmetric Objects (CSOs) are considered the young counterparts of large doubles according to advance speeds measured or inferred from spectral ageing. Here we present a simple power law model for the CSO/FRII evolution based on the study of sources with well defined hot-spots. The luminosity of the hot spots is estimated under minimum energy conditions. The advance of the source is consi… ▽ More Compact Symmetric Objects (CSOs) are considered the young counterparts of large doubles according to advance speeds measured or inferred from spectral ageing. Here we present a simple power law model for the CSO/FRII evolution based on the study of sources with well defined hot-spots. The luminosity of the hot spots is estimated under minimum energy conditions. The advance of the source is considered to proceed in ram pressure equilibrium with the ambient medium. Finally, we also assume that the jets feeding the hot spots are relativistic and have a time dependent power. Comparison with observational data allows to interpret the CSO-FRII evolution in terms of decreasing jet power with time. △ Less

Submitted 28 July, 2004; originally announced July 2004.

Comments: 7 pages

Journal ref: Third GPS/CSS Workshop, 2003 PASA, 20 Number 1 pp.94-97

Stability of hydrodynamical relativistic planar jets. II. Long-term nonlinear evolution

Authors: M. Perucho, J. M. Marti, M. Hanasz

Abstract: In this paper we continue our study of the Kelvin-Helmholtz (KH) instability in relativistic planar jets following the long-term evolution of the numerical simulations which were introduced in Paper I. The models have been classified into four classes (I to IV) with regard to their evolution in the nonlinear phase, characterized by the process of jet/ambient mixing and momentum transfer. Models… ▽ More In this paper we continue our study of the Kelvin-Helmholtz (KH) instability in relativistic planar jets following the long-term evolution of the numerical simulations which were introduced in Paper I. The models have been classified into four classes (I to IV) with regard to their evolution in the nonlinear phase, characterized by the process of jet/ambient mixing and momentum transfer. Models undergoing qualitatively different non-linear evolution are clearly grouped in well-separated regions in a jet Lorentz factor/jet-to-ambient enthalpy diagram. Jets with a low Lorentz factor and small enthalpy ratio are disrupted by a strong shock after saturation. Those with a large Lorentz factor and enthalpy ratio are unstable although the process of mixing and momentum exchange proceeds to a longer time scale due to a steady conversion of kinetic to internal energy in the jet. In these cases, the high value of the initial Lorentz seems to prevent transversal velocity from growing far enough to generate the strong shock that breaks the slower jets. Finally, jets with either high Lorentz factors and small enthalpy ratios or low Lorentz factors and large enthalpy ratios appear as the most stable. △ Less

Submitted 28 July, 2004; originally announced July 2004.

Comments: Paper II, 16 pages. Accepted for publication in Astronomy & Astrophysics. Due to arXiv limits, figures have low quality. Better quality ones will be available in the published paper or on request to manuel.perucho@uv.es

Journal ref: Astron.Astrophys. 427 (2004) 431-444

Stability of hydrodynamical relativistic planar jets. I. Linear evolution and saturation of Kelvin-Helmholtz modes

Authors: M. Perucho, M. Hanasz, J. M. Marti, H. Sol

Abstract: The effects of relativistic dynamics and thermodynamics in the development of Kelvin-Helmholtz instabilities in planar, relativistic jets along the early phases (namely linear and saturation phases) of evolution has been studied by a combination of linear stability analysis and high-resolution numerical simulations for the most unstable first reflection modes in the temporal approach. Three diff… ▽ More The effects of relativistic dynamics and thermodynamics in the development of Kelvin-Helmholtz instabilities in planar, relativistic jets along the early phases (namely linear and saturation phases) of evolution has been studied by a combination of linear stability analysis and high-resolution numerical simulations for the most unstable first reflection modes in the temporal approach. Three different values of the jet Lorentz factor (5, 10 and 20) and a few different values of specific internal energy of the jet matter (from 0.08 to $60.0 c^2$) have been considered. Figures illustrating the evolution of the perturbations are also shown. △ Less

Submitted 27 July, 2004; originally announced July 2004.

Comments: Paper I, 17 pages. Accepted for publication in Astronomy & Astrophysics. Image size and quality are reduced due to arxiv limits. Better quality images will be available in the published paper or on request to manuel.perucho@uv.es

Journal ref: Astron.Astrophys. 427 (2004) 415-429

Assessment of a high-resolution central scheme for the solution of the relativistic hydrodynamics equations

Authors: Arturo Lucas-Serrano, Jose A. Font, Jose M. Ibanez, Jose M. Marti

Abstract: We assess the suitability of a recent high-resolution central scheme developed by Kurganov & Tadmor (2000) for the solution of the relativistic hydrodynamics equations. The novelty of this approach relies on the absence of Riemann solvers in the solution procedure. The computations we present are performed in one and two spatial dimensions in Minkowski spacetime. Standard numerical experiments s… ▽ More We assess the suitability of a recent high-resolution central scheme developed by Kurganov & Tadmor (2000) for the solution of the relativistic hydrodynamics equations. The novelty of this approach relies on the absence of Riemann solvers in the solution procedure. The computations we present are performed in one and two spatial dimensions in Minkowski spacetime. Standard numerical experiments such as shock tubes and the relativistic flat-faced step test are performed. As an astrophysical application the article includes two-dimensional simulations of the propagation of relativistic jets using both Cartesian and cylindrical coordinates. The simulations reported clearly show the capabilities of the numerical scheme to yield satisfactory results, with an accuracy comparable to that obtained by the so-called high-resolution shock-capturing schemes based upon Riemann solvers (Godunov-type schemes), even well inside the ultrarelativistic regime. Such central scheme can be straightforwardly applied to hyperbolic systems of conservation laws for which the characteristic structure is not explicitly known, or in cases where the exact solution of the Riemann problem is prohibitively expensive to compute numerically. Finally, we present comparisons with results obtained using various Godunov-type schemes as well as with those obtained using other high-resolution central schemes which have recently been reported in the literature. △ Less

Submitted 27 July, 2004; originally announced July 2004.

Comments: 14 pages, 12 figures, to appear in A&A

Journal ref: Astron.Astrophys. 428 (2004) 703-715

3D Simulations of Relativistic Precessing Jets Probing the Structure of Superluminal Sources

Authors: Miguel Angel Aloy, Jose Maria Marti, Jose Luis Gomez, Ivan Agudo, Ewald Mueller, Jose Maria Ibanyez

Abstract: We present the results of a three-dimensional, relativistic, hydrodynamic simulation of a precessing jet into which a compact blob of matter is injected. A comparison of synthetic radio maps computed from the hydrodynamic model, taking into account the appropriate light travel time delays, with those obtained from observations of actual superluminal sources shows that the variability of the jet… ▽ More We present the results of a three-dimensional, relativistic, hydrodynamic simulation of a precessing jet into which a compact blob of matter is injected. A comparison of synthetic radio maps computed from the hydrodynamic model, taking into account the appropriate light travel time delays, with those obtained from observations of actual superluminal sources shows that the variability of the jet emission is the result of a complex combination of phase motions, viewing angle selection effects, and non-linear interactions between perturbations and the underlying jet and/or the external medium. These results question the hydrodynamic properties inferred from observed apparent motions and radio structures, and reveal that shock-in-jet models may be overly simplistic. △ Less

Submitted 6 February, 2003; originally announced February 2003.

Comments: Accepted for publication in ApJL. 4 pages, 5 figures (4 in color)

Journal ref: Astrophys.J. 585 (2003) L109-L112

Influence of internal energy on the stability of relativistic flows

Authors: Manuel Perucho, Michal Hanasz, Jose Maria Marti, Helene Sol

Abstract: A set of simulations concerning the influence of internal energy on the stability of relativistic jets is presented. Results show that perturbations saturate when the amplitude of the velocity perturbation approaches the speed of light limit. Also, contrary to what predicted by linear stability theory, jets with higher specific internal energy appear to be more stable. A set of simulations concerning the influence of internal energy on the stability of relativistic jets is presented. Results show that perturbations saturate when the amplitude of the velocity perturbation approaches the speed of light limit. Also, contrary to what predicted by linear stability theory, jets with higher specific internal energy appear to be more stable. △ Less

Submitted 5 December, 2002; originally announced December 2002.

Comments: 4 pages, 2 figures. To be published in proceedings of the V Scientific Meeting of the SEA (Highlights of Spanish Astrophysics III), Toledo, Spain, September 9-13, 2002. Ed. J. Gallego, J. Zamorano and N. Cardiel, ASSL, Kluwer. Uploaded version of figure 3 is low resolution due to size problems: available high resolution version on request to Manuel Perucho (Manuel.Perucho@uv.es)

Physical parameters in the hot spots and jets of Compact Symmetric Objects

Authors: Manuel Perucho, Jose Maria Marti

Abstract: We present a model to determine the physical parameters of jets and hot spots of a sample of CSOs under very basic assumptions like synchrotron emission and minimum energy conditions. Based on this model we propose a simple evolutionary scenario for these sources assuming that they evolve in ram pressure equilibrium with the external medium and constant jet power. The parameters of our model are… ▽ More We present a model to determine the physical parameters of jets and hot spots of a sample of CSOs under very basic assumptions like synchrotron emission and minimum energy conditions. Based on this model we propose a simple evolutionary scenario for these sources assuming that they evolve in ram pressure equilibrium with the external medium and constant jet power. The parameters of our model are constrained from fits of observational data (radio luminosity, hot spot radius and hot spot advance speed) versus projected linear size. From these plots we conclude that CSOs evolve self-similarly and that their radio luminosity increases with linear size along the first kiloparsec. Assuming that the jets feeding CSOs are relativistic from both kinematical and thermodynamical points of view, we use the values of the pressure and particle number density within the hot spots to estimate the fluxes of momentum (thrust), energy, and particles of these relativistic jets. The mean jet power obtained in this way is within an order of magnitude that inferred for FRII sources, which is consistent with CSOs being the possible precursors of large doubles. The inferred flux of particles corresponds to, for a barionic jet, about a 10% of the mass accreted by a black hole of $10^8 {\rm M_{\odot}}$ at the Eddington limit, pointing towards a very efficient conversion of accretion flow into ejection, or to a leptonic composition of jets. △ Less

Submitted 30 November, 2001; originally announced November 2001.

Comments: 11 pages, 2 figures. Accepted for publication in Astrophysical Journal

Does the plasma composition affect the long term evolution of relativistic jets?

Authors: L. Scheck, M. A. Aloy, J. M. Marti, J. L. Gomez, E. Mueller

Abstract: We study the influence of the matter content of extragalactic jets on their morphology, dynamics and emission properties. For this purpose we consider jets of extremely different compositions including pure leptonic and baryonic plasmas. Our work is based on two-dimensional relativistic hydrodynamic simulations of the long-term evolution of powerful extragalactic jets propagating into a homogene… ▽ More We study the influence of the matter content of extragalactic jets on their morphology, dynamics and emission properties. For this purpose we consider jets of extremely different compositions including pure leptonic and baryonic plasmas. Our work is based on two-dimensional relativistic hydrodynamic simulations of the long-term evolution of powerful extragalactic jets propagating into a homogeneous environment. The equation of state used in the simulations accounts for an arbitrary mixture of electrons, protons and electron-positron pairs. Using the hydrodynamic models we have also computed synthetic radio maps and the thermal Bremsstrahlung X-ray emission from their cavities. Although there is a difference of about three orders of magnitude in the temperatures of the cavities inflated by the simulated jets, we find that both the morphology and the dynamic behaviour are almost independent on the assumed composition of the jets. Their evolution proceeds in two distinct epochs. During the first one multidimensional effects are unimportant and the jets propagate ballistically. The second epoch starts when the first larger vortices are produced near the jet head causing the beam cross section to increase and the jet to decelerate. The evolution of the cocoon and cavity is in agreement with a simple theoretical model. The beam velocities are relativistic ($Γ\simeq 4$) at kiloparsec scales supporting the idea that the X-ray emission of several extragalactic jets may be due to relativistically boosted CMB photons. The radio emission of all models is dominated by the contribution of the hot spots. All models exhibit a depression in the X-rays surface brightness of the cavity interior in agreement with recent observations. △ Less

Submitted 7 November, 2002; v1 submitted 19 November, 2001; originally announced November 2001.

Comments: 41 pages, 12 png figures. See also high quality versions of the figures and several movies under http://www.mpa-garching.mpg.de/~maa/JETS/LongTermEv/index.html

Journal ref: Mon.Not.Roy.Astron.Soc. 331 (2002) 615

Polarimetric VLBI observations of 0735+178

Authors: I. Agudo, J. L. Gomez, D. C. Gabuzda, J. C. Guirado, A. Alberdi, A. P. Marscher, M. A. Aloy, J. M. Marti

Abstract: We present a new centimeter polarimetric VLBI image of the BL Lac object 0735+178. This source exhibits one of the most pronounced curvatures observed in jets of AGNs, with two sharp apparent bends of 90 degrees within the inner 2 milliarcseconds from the core. Through the analysis of the data gathered over the past decades we study whether this curvature is produced by precession of the jet wit… ▽ More We present a new centimeter polarimetric VLBI image of the BL Lac object 0735+178. This source exhibits one of the most pronounced curvatures observed in jets of AGNs, with two sharp apparent bends of 90 degrees within the inner 2 milliarcseconds from the core. Through the analysis of the data gathered over the past decades we study whether this curvature is produced by precession of the jet with ejection of ballistic components, or a precession in such a way that components' velocity vector are always parallel to the jet axis. These possibilities are also studied by comparison with 3D hydrodynamic relativistic simulations of precessing jets. △ Less

Submitted 26 January, 2001; originally announced January 2001.

Comments: LaTeX, 4 pages, 4 PostScript figures. To be published in the Proceedings of the 5th EVN Symposium, eds. Conway, J., Polatidis, A. & Booth., R

Jet stability and the generation of superluminal and stationary components

Authors: I. Agudo, J. L. Gomez, J. M. Marti, J. M Ibanez, A. P. Marscher, A. Alberdi, M. A. Aloy, P. E. Hardee

Abstract: We present a numerical simulation of the response of an expanding relativistic jet to the ejection of a superluminal component. The simulation has been performed with a relativistic time-dependent hydrodynamical code from which simulated radio maps are computed by integrating the transfer equations for synchrotron radiation. The interaction of the superluminal component with the underlying jet r… ▽ More We present a numerical simulation of the response of an expanding relativistic jet to the ejection of a superluminal component. The simulation has been performed with a relativistic time-dependent hydrodynamical code from which simulated radio maps are computed by integrating the transfer equations for synchrotron radiation. The interaction of the superluminal component with the underlying jet results in the formation of multiple conical shocks behind the main perturbation. These trailing components can be easily distinguished because they appear to be released from the primary superluminal component, instead of being ejected from the core. Their oblique nature should also result in distinct polarization properties. Those appearing closer to the core show small apparent motions and a very slow secular decrease in brightness, and could be identified as stationary components. Those appearing farther downstream are weaker and can reach superluminal apparent motions. The existence of these trailing components indicates that not all observed components necessarily represent major perturbations at the jet inlet; rather, multiple emission components can be generated by a single disturbance in the jet. While the superluminal component associated with the primary perturbation exhibits a rather stable pattern speed, trailing components have velocities that increase with distance from the core but move at less than the jet speed. The trailing components exhibit motion and structure consistent with the triggering of pinch modes by the superluminal component. △ Less

Submitted 11 January, 2001; originally announced January 2001.

Comments: Accepted by ApJ Letters. LaTeX, 19 pages, 4 PostScript figures

Physical parameters in relativistic jets from Compact Symmetric Objects

Authors: Manuel Perucho, Jose Ma. Marti

Abstract: Compact symmetric objects conform a class of sources characterized by high luminosity radio emission located symmetrically on both sides of the active galactic nucleus on linear scales of less than 1 kpc. Given their small size, the hot spots of the jets in CSOs provide a unique laboratory for the study of the physics of relativistic jets and their environment close to the central engine. We p… ▽ More Compact symmetric objects conform a class of sources characterized by high luminosity radio emission located symmetrically on both sides of the active galactic nucleus on linear scales of less than 1 kpc. Given their small size, the hot spots of the jets in CSOs provide a unique laboratory for the study of the physics of relativistic jets and their environment close to the central engine. We present a simple model for the hot spots in CSOs assuming synchrotron emission, minimum energy and ram pressure equilibrium with the external medium. Further comparison of our model with observational data allows us to constrain the physical parameters in the hot spots and the jets feeding them, and the density profile of the external medium. △ Less

Submitted 15 November, 2000; originally announced November 2000.

Comments: 4 pages, 1 figure, proceedings of the IV scientific meeting of the Sociedad Espanola de Astronomia (SEA)

The exact solution of the Riemann problem with non-zero tangential velocities in relativistic hydrodynamics

Authors: Jose A. Pons, Jose Ma. Marti, Ewald Mueller

Abstract: We have generalised the exact solution of the Riemann problem in special relativistic hydrodynamics for arbitrary tangential flow velocities. The solution is obtained by solving the jump conditions across shocks plus an ordinary differential equation arising from the self-similarity condition along rarefaction waves, in a similar way as in purely normal flow. The dependence of the solution on th… ▽ More We have generalised the exact solution of the Riemann problem in special relativistic hydrodynamics for arbitrary tangential flow velocities. The solution is obtained by solving the jump conditions across shocks plus an ordinary differential equation arising from the self-similarity condition along rarefaction waves, in a similar way as in purely normal flow. The dependence of the solution on the tangential velocities is analysed, and the impact of this result on the development of multidimensional relativistic hydrodynamic codes (of Godunov type) is discussed. △ Less

Submitted 2 May, 2000; originally announced May 2000.

Comments: 26 pages, 4 figures. Accepted for publication in Journal of Fluid Mechanics

Journal ref: J.Fluid Mech. 422 (2000) 125

Relativistic Jets from Collapsars

Authors: Miguel A. Aloy, Ewald Mueller, Jose M. Ibanez, Jose M. Marti, Andrew MacFadyen

Abstract: Using a collapsar progenitor model of MacFadyen & Woosley we have simulated the propagation of an axisymmetric jet through a collapsing rotating massive star with the GENESIS multi-dimensional relativistic hydrodynamic code. The jet forms as a consequence of an assumed (constant or variable) energy deposition in the range $10^{50}$ erg s$^{-1}$ to $10^{51}$ erg s$^{-1}$ within a $30^{\circ}$ con… ▽ More Using a collapsar progenitor model of MacFadyen & Woosley we have simulated the propagation of an axisymmetric jet through a collapsing rotating massive star with the GENESIS multi-dimensional relativistic hydrodynamic code. The jet forms as a consequence of an assumed (constant or variable) energy deposition in the range $10^{50}$ erg s$^{-1}$ to $10^{51}$ erg s$^{-1}$ within a $30^{\circ}$ cone around the rotation axis. The jet flow is strongly beamed ($\la$ few degrees), spatially inhomogeneous, and time dependent. The jet reaches the surface of the stellar progenitor ($R_{\ast} = 2.98 \times 10^{10} $cm) intact. At breakout the maximum Lorentz factor of the jet flow is 33. After breakout the jet accelerates into the circumstellar medium, whose density is assumed to decrease exponentially and then being constant $ρ_{\rm ext} = 10^{-5}$ gcm$^{-3}$. Outside the star the flow begins to expand also laterally ($v \sim c$), but the beam remains very well collimated. At a distance of $2.54 R_{\ast}$, where the simulation ends, the Lorentz factor has increased to 44. △ Less

Submitted 17 January, 2000; v1 submitted 5 November, 1999; originally announced November 1999.

Comments: 7 pages, 4 figures, ApJL accepted

Riemann Solvers in General Relativistic Hydrodynamics

Authors: J. M. Ibanez, M. A. Aloy, J. A. Font, J. M. Marti, J. A. Miralles, J. A. Pons

Abstract: Our contribution concerns with the numerical solution of the 3D general relativistic hydrodynamical system of equations within the framework of the 3+1 formalism. We summarize the theoretical ingredients which are necessary in order to build up a numerical scheme based on the solution of local Riemann problems. Hence, the full spectral decomposition of the Jacobian matrices of the system, i.e.,… ▽ More Our contribution concerns with the numerical solution of the 3D general relativistic hydrodynamical system of equations within the framework of the 3+1 formalism. We summarize the theoretical ingredients which are necessary in order to build up a numerical scheme based on the solution of local Riemann problems. Hence, the full spectral decomposition of the Jacobian matrices of the system, i.e., the eigenvalues and the right and left eigenvectors, is explicitly shown. An alternative approach consists in using any of the special relativistic Riemann solvers recently developed for describing the evolution of special relativistic flows. Our proposal relies on a local change of coordinates in terms of which the spacetime metric is locally Minkowskian and permits an accurate description of numerical general relativistic hydrodynamics. △ Less

Submitted 3 November, 1999; originally announced November 1999.

Comments: 10 pages, to appear in the proceedings of the conference "Godunov methods: theory and applications", Oxford, October 1999

Relativistic Jets from Collapsars

Authors: Miguel A. Aloy, Ewald Mueller, Jose M. Ibanez, Jose M. Marti, Andrew MacFadyen

Abstract: We have studied the relativistic beamed outflow proposed to occur in the collapsar model of gamma-ray bursts. A jet forms as a consequence of an assumed energy deposition of $\sim 10^{50}- 10^{51}$ erg/s within a $30^{\circ}$ cone around the rotation axis of the progenitor star. The generated jet flow is strongly beamed ($\la$ few degrees) and reaches the surface of the stellar progenitor (r… ▽ More We have studied the relativistic beamed outflow proposed to occur in the collapsar model of gamma-ray bursts. A jet forms as a consequence of an assumed energy deposition of $\sim 10^{50}- 10^{51}$ erg/s within a $30^{\circ}$ cone around the rotation axis of the progenitor star. The generated jet flow is strongly beamed ($\la$ few degrees) and reaches the surface of the stellar progenitor (r $\approx 3 10^{10} $cm) intact. At break-out the maximum Lorentz factor of the jet flow is about 33. Simulations have been performed with the GENESIS multi-dimensional relativistic hydrodynamic code. △ Less

Submitted 26 October, 1999; originally announced October 1999.

Comments: 6 pages, 2 figures, to appear in the proceedings of the conference "Godunov methods: theory and applications", Oxford, October 1999

Simulations of Relativistic Jets with GENESIS

Authors: Miguel A. Aloy, Jose M. Ibanez, Jose M. Marti, Jose L. Gomez, Ewald Mueller

Abstract: The multidimensional relativistic hydrodynamical code GENESIS has been used to obtain first results of {\it 3D} simulations of relativistic jets. We have studied the influence of a slight perturbation of the injection velocity field on the morphodynamics of otherwise axisymmetric relativistic jets. The multidimensional relativistic hydrodynamical code GENESIS has been used to obtain first results of {\it 3D} simulations of relativistic jets. We have studied the influence of a slight perturbation of the injection velocity field on the morphodynamics of otherwise axisymmetric relativistic jets. △ Less

Submitted 26 October, 1999; originally announced October 1999.

Comments: 6 pages, 2 figures, to appear in the proceedings of the conference "Godunov methods: theory and applications", Oxford, October 1999

An exact Riemann Solver for multidimensional special relativistic hydrodynamics

Authors: J. Pons, J. Ma. Marti, E. Muller

Abstract: We have generalised the exact solution of the Riemann problem in special relativistic hydrodynamics for arbitrary tangential flow velocities. The solution is obtained by solving the jump conditions across shocks plus an ordinary differential equation arising from the self-similarity condition along rarefaction waves, in a similar way as in purely normal flow. The dependence of the solution on th… ▽ More We have generalised the exact solution of the Riemann problem in special relativistic hydrodynamics for arbitrary tangential flow velocities. The solution is obtained by solving the jump conditions across shocks plus an ordinary differential equation arising from the self-similarity condition along rarefaction waves, in a similar way as in purely normal flow. The dependence of the solution on the tangential velocities is analysed. This solution has been used to build up an exact Riemann solver implemented in a multidimensional relativistic (Godunov-type) hydro-code. △ Less

Submitted 26 October, 1999; originally announced October 1999.

Comments: 8 pages, 2 figures. To appear in the proceedings of the conference "Godunov methods: theory and applications". Oxford, October-1999

Numerical hydrodynamics in special relativity

Authors: Jose M. Marti, Ewald Mueller

Abstract: This review is concerned with a discussion of numerical methods for the solution of the equations of special relativistic hydrodynamics (SRHD). Particular emphasis is put on a comprehensive review of the application of high-resolution shock-capturing methods in SRHD. Results obtained with different numerical SRHD methods are compared, and two astrophysical applications of SRHD flows are discusse… ▽ More This review is concerned with a discussion of numerical methods for the solution of the equations of special relativistic hydrodynamics (SRHD). Particular emphasis is put on a comprehensive review of the application of high-resolution shock-capturing methods in SRHD. Results obtained with different numerical SRHD methods are compared, and two astrophysical applications of SRHD flows are discussed. An evaluation of the various numerical methods is given and future developments are analyzed. △ Less

Submitted 21 June, 1999; originally announced June 1999.

Comments: submitted to the electronic journal Living Reviews in Relativity, at http://www.livingreviews.org 72 pages, 14 postscript figures included

Journal ref: LivingRev.Rel.2:3,1999

GENESIS: A high-resolution code for 3D relativistic hydrodynamics

Authors: M. A. Aloy, J. M. Ibanez, J. M. Marti, E. Muller

Abstract: The main features of a three dimensional, high-resolution special relativistic hydro code based on relativistic Riemann solvers are described. The capabilities and performance of the code are discussed. In particular, we present the results of extensive test calculations which demonstrate that the code can accurately and efficiently handle strong shocks in three spatial dimensions. Results of th… ▽ More The main features of a three dimensional, high-resolution special relativistic hydro code based on relativistic Riemann solvers are described. The capabilities and performance of the code are discussed. In particular, we present the results of extensive test calculations which demonstrate that the code can accurately and efficiently handle strong shocks in three spatial dimensions. Results of the performance of the code on single and multi-processor machines are given. Simulations (in double precision) with 7x10^6 computational cells require less than 1 Gb of RAM memory and ~7x10^-5 CPU seconds per zone and time step (on a SCI Cray-Origin 2000 with a R10000 processor). Currently, a version of the numerical code is under development, which is suited for massively parallel computers with distributed memory architecture (like, e.g., Cray T3E). △ Less

Submitted 23 March, 1999; originally announced March 1999.

Comments: 39 pages, 8 figures, 8 tables. Accepted to be publish in the ApJS. Tar+gzip document

Journal ref: Astrophys.J.Suppl. 122 (1999) 151

General Relativistic Hydrodynamics with Special Relativistic Riemann Solvers

Authors: Jose A. Pons, Jose A. Font, Jose M. Ibanez, Jose M. Marti, Juan A. Miralles

Abstract: We present a general and practical procedure to solve the general relativistic hydrodynamic equations by using any of the special relativistic Riemann solvers recently developed for describing the evolution of special relativistic flows. Our proposal relies on a local change of coordinates in terms of which the spacetime metric is locally Minkowskian and permits accurate numerical calculations o… ▽ More We present a general and practical procedure to solve the general relativistic hydrodynamic equations by using any of the special relativistic Riemann solvers recently developed for describing the evolution of special relativistic flows. Our proposal relies on a local change of coordinates in terms of which the spacetime metric is locally Minkowskian and permits accurate numerical calculations of general relativistic hydrodynamics problems using the numerical tools developed for the special relativistic case with negligible computational cost. The feasibility of the method has been confirmed by a number of numerical experiments. △ Less

Submitted 21 July, 1998; originally announced July 1998.

Comments: 12 pages, 2 figures, accepted for publication in Astronomy and Astrophysics

Journal ref: Astron.Astrophys. 339 (1998) 638-642

Bending relativistic jets in AGNs

Authors: T. Plewa, J. Ma. Marti, E. Mueller, M. Rozyczka, M. Sikora

Abstract: We present simulations of relativistic jets propagating in a nonuniform medium. Specifically, we study the bending of jets propagating obliquely to the vector of the density-gradient. Our results are applied to the NGC 4258, where such a medium is assumed to be provided by the atmosphere of the sub-parsec accretion disk tilted with respect to the original direction of the jet propagation. As a r… ▽ More We present simulations of relativistic jets propagating in a nonuniform medium. Specifically, we study the bending of jets propagating obliquely to the vector of the density-gradient. Our results are applied to the NGC 4258, where such a medium is assumed to be provided by the atmosphere of the sub-parsec accretion disk tilted with respect to the original direction of the jet propagation. As a result, the jet is bent on a scale comparable to the density scaleheight of the disk atmosphere. The magnitude of the bending effect is found to be largest for light jets with low Lorentz factors. The predicted direction of bending is consistent with the observations. △ Less

Submitted 23 June, 1997; v1 submitted 12 June, 1997; originally announced June 1997.

Comments: 6 pages, tex, uses lecproc.cmm, 3 figures (included); also available at http://www.camk.edu.pl/~tomek/papers/rjet_s.uu (reduced resolution) or http://www.camk.edu.pl/~tomek/papers/rjet_l.uu (full resolution); mpeg movies available at http://www.camk.edu.pl/~tomek/RJET/index.html#Bend . To appear in Proceedings of "Relativistic jets in AGNs", Cracow, 27-30 May 1997. Revised version: slightly modified text, corrected parameter set of model A