BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

G. Bodo

G. Bodo contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
astro-ph.HEastro-ph.SRastro-phastro-ph.EPastro-ph.COastro-ph.IMphysics.comp-ph

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
20works
0followers
7topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

20 published item(s)

preprint2025arXiv

How do recollimation-induced instabilities shape the propagation of hydrodynamic relativistic jets?

Recollimation is a phenomenon of particular importance in the dynamic evolution of jets and in the emission of high-energy radiation. Additionally, the full comprehension of this phenomenon provides insights into fundamental properties of jets in the vicinity of the Active Galactic Nucleus (AGN). Three-dimensional (magneto-)hydrodynamic simulations revealed that the jet conditions at recollimation favor the growth of strong instabilities, challenging the traditional view-supported from two-dimensional simulations-of confined jets undergoing a series of recollimation and reflection shocks. To investigate the stability of relativistic jets in AGNs at recollimation sites, we perform a set of long duration three-dimensional relativistic hydrodynamic simulations with the state-of-the-art PLUTO code, to focus on the development of hydrodynamical instabilities. We explore the non-linear growth of the instabilities and their effects on the physical jet properties as a function of the initial jet parameters: jet Lorentz factor, temperature, opening angle and jet-environment density-contrast. The parameter space is designed to describe low-power, weakly magnetized jets at small distances from the core (around the parsec scale). All collimating jets we simulated develop instabilities. Recollimation instabilities decelerate the jet, heat it, entrain external material, and move the recollimation point to shorter distances from the core. This is true for both conical and cylindrical jets. The instabilities, that are first triggered by the centrifugal instability, appear to be less disruptive in the case of narrower, denser, more relativistic, and warmer jets. These results provide valuable insights into the complex processes governing AGN jets and could be used to model the properties of low-power, weakly magnetized jetted AGNs.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Making Fanaroff-Riley I radio sources III. The effects of the magnetic field on relativistic jets' propagation and source morphologies

Extragalactic radio sources appear under different morphologies, the most frequent ones are classified as Fanaroff-Riley type I (FR I), typically with lower luminosities, and Fanaroff-Riley type II, (FR II), typically more luminous. This simple classification, however, has many exceptions that we intend to investigate. Following previous analyses in the three-dimensional Hydrodynamic and Magneto-Hydrodynamic limits, we extend the numerical investigation to the Relativistic Magneto-Hydrodynamic regime, to include sources whose jet kinetic power sets in the range that separates FR Is from FR IIs. We consider weakly and mildly relativistic, underdense, supersonic jets that propagate in a stratified medium. In the model, the ambient temperature increases with distance from the jet origin maintaining constant pressure. We present three cases with low, high and intermediate kinetic luminosity that evolve into different morphologies. We find that the resulting morphology can be highly time dependent and that, apart from the jet power, the jet-to-ambient density ratio and the magnetization parameter play a crucial role in the jet evolution as well.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Linear stability analysis of magnetized jets: the rotating case

We perform a linear stability analysis of magnetized rotating cylindrical jet flows in the approximation of zero thermal pressure. We focus our analysis on the effect of rotation on the current driven mode and on the unstable modes introduced by rotation. We find that rotation has a stabilizing effect on the current driven mode only for rotation velocities of the order of the Alfvén velocity. Rotation introduces also a new unstable centrifugal buoyancy mode and the "cold" magnetorotational instability. The first mode is analogous to the Parker instability with the centrifugal force playing the role of effective gravity. The magnetorotational instability can be present, but only in a very limited region of the parameter space and is never dominant. The current driven mode is characterized by large wavelenghts and is dominant at small values of the rotational velocity, while the buoyancy mode becomes dominant as rotation is increased and is characterized by small wavelenghts.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Making Fanaroff-Riley I radio sources. Numerical Hydrodynamic 3D Simulations of Low Power Jets

Extragalactic radio sources have been classified into two classes, Fanaroff-Riley I and II, which differ in morphology and radio power. Strongly emitting sources belong to the edge-brightened FR II class, and weakly emitting sources to the edge-darkened FR I class. The origin of this dichotomy is not yet fully understood. Numerical simulations are successful in generating FR II morphologies, but they fail to reproduce the diffuse structure of FR Is. By means of hydro-dynamical 3D simulations of supersonic jets, we investigate how the displayed morphologies depend on the jet parameters. Bow shocks and Mach disks at the jet head, which are probably responsible for the hot spots in the FR II sources, disappear for a jet kinetic power L_kin < 10^43 erg/s. This threshold compares favorably with the luminosity at which the FR I/FR II transition is observed. The problem is addressed by numerical means carrying out 3D HD simulations of supersonic jets that propagate in a non-homogeneous medium with the ambient temperature that increases with distance from the jet origin, which maintains constant pressure. The jet energy in the lower power sources, instead of being deposited at the terminal shock, is gradually dissipated by the turbulence. The jets spread out while propagating, and they smoothly decelerate while mixing with the ambient medium and produce the plumes characteristic of FR I objects. Three-dimensionality is an essential ingredient to explore the FR I evolution because the properties of turbulence in two and three dimensions are very different, since there is no energy cascade to small scales in two dimensions, and two-dimensional simulations with the same parameters lead to FRII-like behavior.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Astrophysical fluid simulations of thermally ideal gases with non-constant adiabatic index: numerical implementation

An Equation of State (\textit{EoS}) closes the set of fluid equations. Although an ideal EoS with a constant \textit{adiabatic index} $Γ$ is the preferred choice due to its simplistic implementation, many astrophysical fluid simulations may benefit from a more sophisticated treatment that can account for diverse chemical processes. Here, we first review the basic thermodynamic principles of a gas mixture in terms of its thermal and caloric EoS by including effects like ionization, dissociation as well as temperature dependent degrees of freedom such as molecular vibrations and rotations. The formulation is revisited in the context of plasmas that are either in equilibrium conditions (local thermodynamic- or collisional excitation- equilibria) or described by non-equilibrium chemistry coupled to optically thin radiative cooling. We then present a numerical implementation of thermally ideal gases obeying a more general caloric EoS with non-constant adiabatic index in Godunov-type numerical schemes.We discuss the necessary modifications to the Riemann solver and to the conversion between total energy and pressure (or vice-versa) routinely invoked in Godunov-type schemes. We then present two different approaches for computing the EoS.The first one employs root-finder methods and it is best suited for EoS in analytical form. The second one leans on lookup table and interpolation and results in a more computationally efficient approach although care must be taken to ensure thermodynamic consistency. A number of selected benchmarks demonstrate that the employment of a non-ideal EoS can lead to important differences in the solution when the temperature range is $500-10^4$ K where dissociation and ionization occur. The implementation of selected EoS introduces additional computational costs although using lookup table methods can significantly reduce the overhead by a factor $3\sim 4$.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Global Properties of Fully Convective Accretion Disks from Local Simulations

We present an approach to deriving global properties of accretion disks from the knowledge of local solutions derived from numerical simulations based on the shearing box approximation. The approach consists of a two-step procedure. First a local solution valid for all values of the disk height is constructed by piecing together an interior solution obtained numerically with an analytical exterior radiative solution. The matching is obtained by assuming hydrostatic balance and radiative equilibrium. Although in principle the procedure can be carried out in general, it simplifies considerably when the interior solution is fully convective. In these cases, the construction is analogous to the derivation of the Hayashi tracks for protostars. The second step consists of piecing together the local solutions at different radii to obtain a global solution. Here we use the symmetry of the solutions with respect to the defining dimensionless numbers--in a way similar to the use of homology relations in stellar structure theory--to obtain the scaling properties of the various disk quantities with radius.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Linear and nonlinear evolution of current-carrying highly magnetized jets

We investigate the linear and nonlinear evolution of current-carrying jets in a periodic configuration by means of high resolution three-dimensional numerical simulations. The jets under consideration are strongly magnetized with a variable pitch profile and initially in equilibrium under the action of a force-free magnetic field. The growth of current-driven (CDI) and Kelvin-Helmholtz (KHI) instabilities is quantified using three selected cases corresponding to static, Alfvenic and super-Alfvenic jets. During the early stages, we observe large-scale helical deformations of the jet corresponding to the growth of the initially excited CDI mode. A direct comparison between our simulation results and the analytical growth rates obtained from linear theory reveals good agreement on condition that high-resolution and accurate discretization algorithms are employed. After the initial linear phase, the jet structure is significantly altered and, while slowly-moving jets show increasing helical deformations, larger velocity shear are violently disrupted on a few Alfven crossing time leaving a turbulent flow structure. Overall, kinetic and magnetic energies are quickly dissipated into heat and during the saturated regime the jet momentum is redistributed on a larger surface area with most of the jet mass travelling at smaller velocities. The effectiveness of this process is regulated by the onset of KHI instabilities taking place at the jet/ambient interface and can be held responsible for vigorous jet braking and entrainment.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

On the convergence of Magnetorotational turbulence in stratified isothermal shearing boxes

We consider the problem of convergence in stratified isothermal shearing boxes with zero net magnetic flux. We present results with the highest resolution to-date--up to 200 grid-point per pressure scale height--that show no clear evidence of convergence. Rather, the Maxwell stresses continue to decrease with increasing resolution. We propose some possible scenarios to explain the lack of convergence based on multi-layer dynamo systems.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Fully Convective Magnetorotational Turbulence in Stratified Shearing Boxes

We present a numerical study of turbulence and dynamo action in stratified shearing boxes with zero magnetic flux. We assume that the fluid obeys the perfect gas law and has finite (constant) thermal diffusivity. We choose radiative boundary conditions at the vertical boundaries in which the heat flux is propor- tional to the fourth power of the temperature. We compare the results with the corresponding cases in which fixed temperature boundary conditions are applied. The most notable result is that the formation of a fully convective state in which the density is nearly constant as a function of height and the heat is transported to the upper and lower boundaries by overturning motions is robust and persists even in cases with radiative boundary conditions. Interestingly, in the convective regime, although the diffusive transport is negligible the mean stratification does not relax to an adiabatic state.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Linear stability analysis of magnetized relativistic jets: the nonrotating case

We perform a linear analysis of the stability of a magnetized relativistic non-rotating cylindrical flow in the aproximation of zero thermal pressure, considering only the m = 1 mode. We find that there are two modes of instability: Kelvin-Helmholtz and current driven. The Kelvin-Helmholtz mode is found at low magnetizations and its growth rate depends very weakly on the pitch parameter. The current driven modes are found at high magnetizations and the value of the growth rate and the wavenumber of the maximum increase as we decrease the pitch parameter. In the relativistic regime the current driven mode is splitted in two branches, the branch at high wavenumbers is characterized by the eigenfunction concentrated in the jet core, the branch at low wavenumbers is instead characterized by the eigenfunction that extends outside the jet velocity shear region.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Revisiting linear dynamics of non-axisymmetric perturbations in weakly magnetized accretion discs

We investigate linear dynamics of non-axisymmetric perturbations in incompressible, vertically stratified Keplerian discs with a weak vertical magnetic field in the shearing box approximation. Perturbations are decomposed into shearing waves whose evolution is followed via numerical integration of the linearized ideal MHD equations. There are two basic modes in the system -- inertia-gravity waves and magnetic mode, which displays the magnetorotational instability (MRI). As distinct from previous studies, we introduce `eigenvariables&#39; characterizing each (counter-propagating) component of the inertia-gravity and magnetic modes, which are governed by a set of four first order coupled ordinary differential equations. This allowed us to identify a new process of linear coupling of the two above non-axisymmetric modes due to the disc&#39;s differential rotation. We did a comparative analysis of the dynamics of non-axisymmetric and axisymmetric magnetic mode perturbations. It is shown that the growth of optimal and close-to-optimal non-axisymmetric harmonics of this mode, having transient nature, can prevail over the exponential growth of axisymmetric ones (i.e., over the axisymmetric MRI) during dynamical time. A possible implication of this result for axisymmetric channel solutions is discussed. Specifically, the formation of the channel may be affected/impeded by non-axisymmetric modes already at the early linear stage leading to its untimely disruption -- the outcome strongly depends on the amplitude and spectrum of initial perturbation. So, this competition may result in an uncertainty in the magnetic mode&#39;s non-linear dynamics. It is also shown that a maximum non-axisymmetric growth is at vertical wavelengths close to the scale-height for which compressibility effects are important. This indirectly suggests that compressibility plays a role in the dynamics of the non-axisymmetric MRI.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Magnetorotational turbulence in stratified shearing boxes with perfect gas equation of state and finite thermal diffusivity

We present a numerical study of turbulence and dynamo action in stratified shearing boxes with zero mean magnetic flux. We assume that the fluid obeys the perfect gas law and has finite (constant) thermal diffusivity. The calculations begin from an isothermal state spanning three scale heights above and below the mid-plane. After a long transient the layers settle to a stationary state in which thermal losses out of the boundaries are balanced by dissipative heating. We identify two regimes. A conductive regime in which the heat is transported mostly by conduction and the density decreases with height. In the limit of large thermal diffusivity this regime resembles the more familiar isothermal case. Another, the convective regime, observed at smaller values of the thermal diffusivity, in which the layer becomes unstable to overturning motions, the heat is carried mostly by advection and the density becomes nearly constant throughout the layer. In this latter constant-density regime we observe evidence for large-scale dynamo action leading to a substantial increase in transport efficiency relative to the conductive cases.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The PLUTO Code for Adaptive Mesh Computations in Astrophysical Fluid Dynamics

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

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

High-order conservative finite difference GLM-MHD schemes for cell-centered MHD

We present and compare third- as well as fifth-order accurate finite difference schemes for the numerical solution of the compressible ideal MHD equations in multiple spatial dimensions. The selected methods lean on four different reconstruction techniques based on recently improved versions of the weighted essentially non-oscillatory (WENO) schemes, monotonicity preserving (MP) schemes as well as slope-limited polynomial reconstruction. The proposed numerical methods are highly accurate in smooth regions of the flow, avoid loss of accuracy in proximity of smooth extrema and provide sharp non-oscillatory transitions at discontinuities. We suggest a numerical formulation based on a cell-centered approach where all of the primary flow variables are discretized at the zone center. The divergence-free condition is enforced by augmenting the MHD equations with a generalized Lagrange multiplier yielding a mixed hyperbolic/parabolic correction, as in Dedner et al. (J. Comput. Phys. 175 (2002) 645-673). The resulting family of schemes is robust, cost-effective and straightforward to implement. Compared to previous existing approaches, it completely avoids the CPU intensive workload associated with an elliptic divergence cleaning step and the additional complexities required by staggered mesh algorithms. Extensive numerical testing demonstrate the robustness and reliability of the proposed framework for computations involving both smooth and discontinuous features.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

High Resolution 3D Relativistic MHD Simulations of Jets

Relativistic magnetized jets are key elements in Active Galactic Nuclei and in other astrophysical environments. Their structure and evolution involves a complex nonlinear physics that can be approached by numerical studies only. Still, owing to a number of challenging computational aspects, only a few numerical investigations have been undertaken so far. In this paper, we present high-resolution three dimensional numerical simulations of relativistic magnetized jets carrying an initially toroidal magnetic field. The presence of a substantial toroidal component of the field is nowadays commonly invoked and held responsible for the process of jet acceleration and collimation. We find that the typical nose cone structures, commonly observed in axisymmetric two-dimensional simulations, are not produced in the 3D case. Rather, the toroidal field gives rise to strong current driven kink instabilities leading to jet wiggling. However, it appears to be able to maintain an highly relativistic spine along its full length. By comparing low and high resolution simulations, we emphasize the impact of resolution on the jet dynamical properties.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Linear coupling of modes in 2D radially stratified astrophysical discs

We investigate mode coupling in a two dimensional compressible disc with radial stratification and differential rotation. We employ the global radial scaling of linear perturbations and study the linear modes in the local shearing sheet approximation. We employ a three-mode formalism and study the vorticity (W), entropy (S) and compressional (P) modes and their coupling properties. The system exhibits asymmetric three-mode coupling: these include mutual coupling of S and P-modes, S and W-modes, and asymmetric coupling between the W and P-modes. P-mode perturbations are able to generate potential vorticity through indirect three-mode coupling. This process indicates that compressional perturbations can lead to the development of vortical structures and influence the dynamics of radially stratified hydrodynamic accretion and protoplanetary discs.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Time-dependent MHD shocks and line intensity ratios in the HH 30 jet: A focus on cooling function and numerical resolution

The coupling between time-dependent, multidimensional MHD numerical codes and radiative line emission is of utmost importance in the studies of the interplay between dynamical and radiative processes in many astrophysical environments, with particular interest for problems involving radiative shocks. There is a widespread consensus that line emitting knots observed in Herbig-Haro jets can be interpreted as radiative shocks. In this paper we address two different aspects relevant to the time-dependent calculations of the line intensity ratios of forbidden transitions, resulting from the excitation by planar, time-dependent radiative shocks traveling in a stratified medium. The first one concerns the impact of the radiation and ionization processes included in the cooling model, and the second one the effects of the numerical grid resolution. In this paper we apply the AMR methodology to the treatment of radiating shocks and show how this method is able to vastly reduce the integration time. The technique is applied to the knots of the HH 30 jet to obtain the observed line intensity ratios and derive the physical parameters, such as density, temperature and ionization fraction. We consider the impact of two different cooling functions and different grid resolutions on the results. We conclude that the use of different cooling routines has effects on results whose weight depends upon the line ratio considered. Moreover, we find the minimum numerical resolution of the simulation grid behind the shock to achieve convergence in the results. This is crucial for the forthcoming 2D calculations of radiative shocks.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

A five-wave HLL Riemann solver for relativistic MHD

We present a five-wave Riemann solver for the equations of ideal relativistic magnetohydrodynamics. Our solver can be regarded as a relativistic extension of the five-wave HLLD Riemann solver initially developed by Miyoshi and Kusano for the equations of ideal MHD. The solution to the Riemann problem is approximated by a five wave pattern, comprised of two outermost fast shocks, two rotational discontinuities and a contact surface in the middle. The proposed scheme is considerably more elaborate than in the classical case since the normal velocity is no longer constant across the rotational modes. Still, proper closure to the Rankine-Hugoniot jump conditions can be attained by solving a nonlinear scalar equation in the total pressure variable which, for the chosen configuration, has to be constant over the whole Riemann fan. The accuracy of the new Riemann solver is validated against one dimensional tests and multidimensional applications. It is shown that our new solver considerably improves over the popular HLL solver or the recently proposed HLLC schemes.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

An HLLC Riemann Solver for Relativistic Flows: I. Hydrodynamics

We present an extension of the HLLC approximate Riemann solver by Toro, Spruce and Speares to the relativistic equations of fluid dynamics. The solver retains the simplicity of the original two-wave formulation proposed by Harten, Lax and van Leer (HLL) but it restores the missing contact wave in the solution of the Riemann problem. The resulting numerical scheme is computationally efficient, robust and positively conservative. The performance of the new solver is evaluated through numerical testing in one and two dimensions.

0 citations0 reviewsNo code linkedOpen access
preprint1997arXiv

BeppoSAX observations of low power radio galaxies: possible detection of obscured nuclei

We present the first results of BeppoSAX observations of a small sample of low brightness FRI radio galaxies. The flux of all the targets is consistent with a thermal spectrum, as due to the presence of hot intracluster gas or galactic corona. Moreover in three sources a non thermal absorbed spectrum can be present in the MECS spectrum at energies larger than 7 keV, while for a fourth object a high energy flux has been detected in the PDS instrument at energies larger than 15 keV. This component could be related to the inner AGN surrounded by an obscuring torus.

0 citations0 reviewsNo code linkedOpen access