-
A Fourth-Order Finite Volume Scheme for Resistive Relativistic Magnetohydrodynamics
Authors:
Andrea Mignone,
Vittoria Berta,
Marco Rossazza,
Matteo Bugli,
Giancarlo Mattia,
Luca Del Zanna,
Lorenzo Pareschi
Abstract:
We present a finite-volume, genuinely 4th-order accurate numerical method for solving the equations of resistive relativistic magnetohydrodynamics (Res-RMHD) in Cartesian coordinates. In our formulation, the magnetic field is evolved in time in terms of face-average values via the constrained-transport method while the remaining variables (density, momentum, energy and electric fields) are advance…
▽ More
We present a finite-volume, genuinely 4th-order accurate numerical method for solving the equations of resistive relativistic magnetohydrodynamics (Res-RMHD) in Cartesian coordinates. In our formulation, the magnetic field is evolved in time in terms of face-average values via the constrained-transport method while the remaining variables (density, momentum, energy and electric fields) are advanced as cell volume-averages. Spatial accuracy employs 5th-order accurate WENO-Z reconstruction from point values (as described in a companion paper) to obtain left and right states at zone interfaces. Explicit flux evaluation is carried out by solving a Riemann problem at cell interfaces, using the Maxwell-Harten-Lax-van Leer with contact wave resolution (MHLLC). Time stepping is based on the implicit-explicit (IMEX) Runge-Kutta (RK) methods, of which we consider both the 3rd-order strong stability preserving SSP3(4,3,3) and a recent 4th-order additive RK scheme, to cope with the stiffness introduced by the source term in Ampere's law. Numerical benchmarks are presented in order to assess the accuracy and robustness of our implementation.
△ Less
Submitted 11 July, 2024;
originally announced July 2024.
-
A $4^{\rm th}$-order accurate finite volume method for ideal classical and special relativistic MHD based on pointwise reconstructions
Authors:
Vittoria Berta,
Andrea Mignone,
Matteo Bugli,
Giancarlo Mattia
Abstract:
We present a novel implementation of a genuinely $4^{\rm th}$-order accurate finite volume scheme for multidimensional classical and special relativistic magnetohydrodynamics (MHD) based on the constrained transport (CT) formalism. The scheme introduces several novel aspects when compared to its predecessors yielding a more efficient computational tool. Among the most relevant ones, our scheme exp…
▽ More
We present a novel implementation of a genuinely $4^{\rm th}$-order accurate finite volume scheme for multidimensional classical and special relativistic magnetohydrodynamics (MHD) based on the constrained transport (CT) formalism. The scheme introduces several novel aspects when compared to its predecessors yielding a more efficient computational tool. Among the most relevant ones, our scheme exploits pointwise to pointwise reconstructions (rather than one-dimensional finite volume ones), employs the generic upwind constrained transport averaging and sophisticated limiting strategies that include both a discontinuity detector and an order reduction procedure. Selected numerical benchmarks demonstrate the accuracy and robustness of the method.
△ Less
Submitted 18 December, 2023; v1 submitted 18 October, 2023;
originally announced October 2023.
-
Resistive relativistic MHD simulations of astrophysical jets
Authors:
Giancarlo Mattia,
Luca Del Zanna,
Matteo Bugli,
Andrea Pavan,
Riccardo Ciolfi,
Gianluigi Bodo,
Andrea Mignone
Abstract:
Aims. The main goal of the present paper is to provide the first systematic numerical study of the propagation of astrophysical relativistic jets, in the context of high-resolution shock-capturing resistive relativistic magnetohydrodynamics (RRMHD) simulations. We aim at investigating different values and models for the plasma resistivity coefficient, and at assessing their impact on the level of…
▽ More
Aims. The main goal of the present paper is to provide the first systematic numerical study of the propagation of astrophysical relativistic jets, in the context of high-resolution shock-capturing resistive relativistic magnetohydrodynamics (RRMHD) simulations. We aim at investigating different values and models for the plasma resistivity coefficient, and at assessing their impact on the level of turbulence, the formation of current sheets and reconnection plasmoids, the electromagnetic energy content, and the dissipated power. Methods. We use the PLUTO code for simulations and we assume an axisymmetric setup for jets, endowed with both poloidal and toroidal magnetic fields, and propagating in a uniform magnetized medium. The gas is assumed to be characterized by a realistic Synge-like equation of state (Taub equation), appropriate for such type of astrophysical jets. The Taub equation is combined here for the first time with the Implicit-Explicit Runge-Kutta time-stepping procedure, as required in RRMHD simulations. Results. The main result is that turbulence is clearly suppressed for the highest values of resistivity (low Lundquist numbers), current sheets are broader, and plasmoids are barely present, while for low values of resistivity results are very similar to ideal runs, where dissipation is purely numerical. We find that recipes employing a variable resistivity based on the advection of a jet tracer or on the assumption of a uniform Lundquist number improve on the use of a constant coefficient and are probably more realistic, preserving the development of turbulence and of sharp current sheets, possible sites for the acceleration of the non-thermal particles producing the observed high-energy emission.
△ Less
Submitted 12 September, 2023; v1 submitted 18 August, 2023;
originally announced August 2023.
-
Jets from accretion disk dynamos: consistent quenching modes for dynamo and resistivity
Authors:
Giancarlo Mattia,
Christian Fendt
Abstract:
Astrophysical jets are launched from strongly magnetized systems that host an accretion disk surrounding a central object. The origin of the magnetic field, which is a key component of the launching process, is still an open question. Here we address the question of how the magnetic field required for jet launching is generated and maintained by a dynamo process. By carrying out non-ideal MHD simu…
▽ More
Astrophysical jets are launched from strongly magnetized systems that host an accretion disk surrounding a central object. The origin of the magnetic field, which is a key component of the launching process, is still an open question. Here we address the question of how the magnetic field required for jet launching is generated and maintained by a dynamo process. By carrying out non-ideal MHD simulations (PLUTO code), we investigate how the feedback of the generated magnetic field on the mean-field dynamo affects the disk and jet properties. We find that a stronger quenching of the dynamo leads to a saturation of the magnetic field at a lower disk magnetization. Nevertheless, we find that, while applying different dynamo feedback models, the overall jet properties remain unaffected. We then investigate a feedback model which encompasses a quenching of the magnetic diffusivity. Our modeling considers a more consistent approach for mean-field dynamo modeling simulations, as the magnetic quenching of turbulence should be considered for both, a turbulent dynamo and turbulent magnetic diffusivity. We find that, after the magnetic field is saturated, the Blandford-Payne mechanism can work efficiently, leading to more collimated jets, that move, however, with slower speed. We find strong intermittent periods of flaring and knot ejection for low Coriolis numbers. In particular, flux ropes are built up and advected towards the inner disk thereby cutting off of the inner disk wind, leading to magnetic field reversals, reconnection and the emergence of intermittent flares.
△ Less
Submitted 2 July, 2022;
originally announced July 2022.
-
Exploring the disk-jet connection in NGC 315
Authors:
L. Ricci,
B. Boccardi,
E. Nokhrina,
M. Perucho,
N. MacDonald,
G. Mattia,
P. Grandi,
E. Madika,
T. P. Krichbaum,
J. A. Zensus
Abstract:
Aims. Hot accretion flows are thought to be able to power the relativistic jets observed in Active Galactic Nuclei. They can present themselves as SANE (Standard And Normal Evolution) disks or MAD (Magnetically Arrested Disks), two states implying profound differences in the physical properties of the disks themselves and of the outflows they produce. Methods. In this paper we use a multi-frequenc…
▽ More
Aims. Hot accretion flows are thought to be able to power the relativistic jets observed in Active Galactic Nuclei. They can present themselves as SANE (Standard And Normal Evolution) disks or MAD (Magnetically Arrested Disks), two states implying profound differences in the physical properties of the disks themselves and of the outflows they produce. Methods. In this paper we use a multi-frequency and multi-epoch data set to study the giant radio galaxy NGC 315, with the goal to explore the properties of its accretion disk and sub-parsec jet. We analyze the source maps with a pixel-based analysis and we use theoretical models to link the observational properties of the jet to the physical state of the accretion disk. Results. We propose that the bulk flow in NGC 315 accelerates on sub-pc scales, concurrently with the parabolic expansion. We show that this fast acceleration can be theoretically reconciled with a magnetically driven acceleration. Along the acceleration and collimation zone, we observe an unexpected spectral behavior, with very steep spectral index values $α\sim -1.5$ ($S_ν\propto ν^α$) between 22 GHz and 43 GHz. Based on the properties of this region, we predict the black hole of NGC 315 to be fast rotating and the magnetic flux threading the accretion disk to be in excellent agreement with that expected in the case of a MAD. Using a new formalism based on the core-shift effect, we model the magnetic field downstream a quasi-parabolic accelerating jet and we reconstruct it up to the event horizon radius. In the MAD scenario, we compare it with the expected magnetic saturation strengths in the disk, finding a good agreement.
△ Less
Submitted 24 June, 2022;
originally announced June 2022.
-
A comparison of approximate non-linear Riemann solvers for Relativistic MHD
Authors:
Giancarlo Mattia,
Andrea Mignone
Abstract:
We compare a particular selection of approximate solutions of the Riemann problem in the context of ideal relativistic magnetohydrodynamics. In particular, we focus on Riemann solvers not requiring a full eigenvector structure. Such solvers recover the solution of the Riemann problem by solving a simplified or reduced set of jump conditions, whose level of complexity depends on the intermediate mo…
▽ More
We compare a particular selection of approximate solutions of the Riemann problem in the context of ideal relativistic magnetohydrodynamics. In particular, we focus on Riemann solvers not requiring a full eigenvector structure. Such solvers recover the solution of the Riemann problem by solving a simplified or reduced set of jump conditions, whose level of complexity depends on the intermediate modes that are included. Five different approaches - namely the HLL, HLLC, HLLD, HLLEM and GFORCE schemes - are compared in terms of accuracy and robustness against one- and multi-dimensional standard numerical benchmarks. Our results demonstrate that - for weak or moderate magnetizations - the HLLD Riemann solver yields the most accurate results, followed by HLLC solver(s). The GFORCE approach provides a valid alternative to the HLL solver being less dissipative and equally robust for strongly magnetized environments. Finally, our tests show that the HLLEM Riemann solver is not cost-effective in improving the accuracy of the solution and reducing the numerical dissipation.
△ Less
Submitted 23 November, 2021; v1 submitted 17 November, 2021;
originally announced November 2021.
-
MHD accretion-ejection: jets launched by a non-isotropic accretion disk dynamo. II. A dynamo tensor defined by the disk Coriolis number
Authors:
Giancarlo Mattia,
Christian Fendt
Abstract:
Astrophysical jets are launched from strongly magnetized systems that host an accretion disk surrounding a central object. Here we address the question how to generate the accretion disk magnetization and field structure required for jet launching. We continue our work from Paper I (Mattia & Fendt 2020a), considering a non-scalar accretion disk mean-field $α^2Ω$-dynamo in the context of large scal…
▽ More
Astrophysical jets are launched from strongly magnetized systems that host an accretion disk surrounding a central object. Here we address the question how to generate the accretion disk magnetization and field structure required for jet launching. We continue our work from Paper I (Mattia & Fendt 2020a), considering a non-scalar accretion disk mean-field $α^2Ω$-dynamo in the context of large scale disk-jet simulations. We now investigate a disk dynamo that follows analytical solutions of mean-field dynamo theory, essentially based only on a single parameter, the Coriolis number. We thereby confirm the anisotropy of the dynamo tensor acting in accretion disks, allowing to relate both the resistivity and mean-field dynamo to the disk turbulence. Our new model recovers previous simulations applying a purely radial initial field, while allowing for a more stable evolution for seed fields with a vertical component. We also present correlations between the strength of the disk dynamo coefficients and the dynamical parameters of the jet that is launched, and discuss their implication for observed jet quantities.
△ Less
Submitted 3 August, 2020;
originally announced August 2020.
-
MHD accretion-ejection: jets launched by a non-isotropic accretion disk dynamo. I. Validation and application of selected dynamo tensorial components
Authors:
Giancarlo Mattia,
Christian Fendt
Abstract:
Astrophysical jets are launched from strongly magnetized systems that host an accretion disk surrounding a central object. The origin of the jet launching magnetic field is one of the open questions for modeling the accretion-ejection process. Here we address the question how to generate the accretion disk magnetization and field structure required for jet launching. Applying the PLUTO code, we pr…
▽ More
Astrophysical jets are launched from strongly magnetized systems that host an accretion disk surrounding a central object. The origin of the jet launching magnetic field is one of the open questions for modeling the accretion-ejection process. Here we address the question how to generate the accretion disk magnetization and field structure required for jet launching. Applying the PLUTO code, we present the first resistive MHD simulations of jet launching including a non-scalar accretion disk mean-field $α^2Ω$-dynamo in the context of large scale disk-jet simulations. Essentially, we find the $α_φ$-dynamo component determining the amplification of the poloidal magnetic field, which is strictly related to the disk magnetization (and, as a consequence, to the jet speed, mass and collimation), while the $α_R$ and $α_θ$-dynamo components trigger the formation of multiple, anti-aligned magnetic loops in the disk, with strong consequences on the stability and dynamics of the disk-jet system. In particular, such loops trigger the formation of dynamo inefficient zones, which are characterized by a weak magnetic field, and therefore a lower value of the magnetic diffusivity. The jet mass, speed and collimation are strongly affected by the formation of the dynamo inefficient zones. Moreover, the $θ$-component of the $α$-dynamo plays a key role when interacting with a non-radial component of the seed magnetic field. We also present correlations between the strength of the disk toy dynamo coefficients and the dynamical parameters of the jet that is launched.
△ Less
Submitted 30 July, 2020;
originally announced July 2020.
-
A Constrained Transport Method for the Solution of the Resistive Relativistic MHD Equations
Authors:
A. Mignone,
G. Mattia,
G. Bodo,
L. Del Zanna
Abstract:
We describe a novel Godunov-type numerical method for solving the equations of resistive relativistic magnetohydrodynamics. In the proposed approach, the spatial components of both magnetic and electric fields are located at zone interfaces and are evolved using the constrained transport formalism. Direct application of Stokes' theorem to Faraday's and Ampere's laws ensures that the resulting disc…
▽ More
We describe a novel Godunov-type numerical method for solving the equations of resistive relativistic magnetohydrodynamics. In the proposed approach, the spatial components of both magnetic and electric fields are located at zone interfaces and are evolved using the constrained transport formalism. Direct application of Stokes' theorem to Faraday's and Ampere's laws ensures that the resulting discretization is divergence-free for the magnetic field and charge-conserving for the electric field. Hydrodynamic variables retain, instead, the usual zone-centred representation commonly adopted in finite-volume schemes. Temporal discretization is based on Runge-Kutta implicit-explicit (IMEX) schemes in order to resolve the temporal scale disparity introduced by the stiff source term in Ampere's law. The implicit step is accomplished by means of an improved and more efficient Newton-Broyden multidimensional root-finding algorithm. The explicit step relies on a multidimensional Riemann solver to compute the line-averaged electric and magnetic fields at zone edges and it employs a one-dimensional Riemann solver at zone interfaces to update zone-centred hydrodynamic quantities. For the latter, we introduce a five-wave solver based on the frozen limit of the relaxation system whereby the solution to the Riemann problem can be decomposed into an outer Maxwell solver and an inner hydrodynamic solver. A number of numerical benchmarks demonstrate that our method is superior in stability and robustness to the more popular charge-conserving divergence cleaning approach where both primary electric and magnetic fields are zone-centered. In addition, the employment of a less diffusive Riemann solver noticeably improves the accuracy of the computations.
△ Less
Submitted 2 April, 2019;
originally announced April 2019.
-
Linear Wave Propagation for Resistive Relativistic Magnetohydrodynamics
Authors:
A. Mignone,
G. Mattia,
G. Bodo
Abstract:
We present a linear mode analysis of the relativistic MHD equations in the presence of finite electrical conductivity. Starting from the fully relativistic covariant formulation, we derive the dispersion relation in the limit of small linear perturbations. It is found that the system supports ten wave modes which can be easily identified in the limits of small or large conductivities. In the resis…
▽ More
We present a linear mode analysis of the relativistic MHD equations in the presence of finite electrical conductivity. Starting from the fully relativistic covariant formulation, we derive the dispersion relation in the limit of small linear perturbations. It is found that the system supports ten wave modes which can be easily identified in the limits of small or large conductivities. In the resistive limit, matter and electromagnetic fields decouple and solution modes approach pairs of light and acoustic waves as well as a number of purely damped (non-propagating) modes. In the opposite (ideal) limit, the frozen-in condition applies and the modes of propagation coincide with a pair of fast magnetosonic, a pair of slow and Alfvén modes, as expected. In addition, the contact mode is always present and it is unaffected by the conductivity. For finite values of the conductivity, the dispersion relation gives rise to either pairs of opposite complex conjugate roots or purely imaginary (damped) modes. In all cases, the system is dissipative and also dispersive as the phase velocity depends nonlineary on the wavenumber. Occasionally, the group velocity may exceed the speed of light although this does not lead to superluminal signal propagation.
△ Less
Submitted 4 September, 2018;
originally announced September 2018.
-
A Particle Module for the PLUTO Code: I - an implementation of the MHD-PIC equations
Authors:
A. Mignone,
G. Bodo,
B. Vaidya,
G. Mattia
Abstract:
We describe an implementation of a particle physics module available for the PLUTO code, appropriate for the dynamical evolution of a plasma consisting of a thermal fluid and a non-thermal component represented by relativistic charged particles, or cosmic rays (CR). While the fluid is approached using standard numerical schemes for magnetohydrodynamics, CR particles are treated kinetically using c…
▽ More
We describe an implementation of a particle physics module available for the PLUTO code, appropriate for the dynamical evolution of a plasma consisting of a thermal fluid and a non-thermal component represented by relativistic charged particles, or cosmic rays (CR). While the fluid is approached using standard numerical schemes for magnetohydrodynamics, CR particles are treated kinetically using conventional Particle-In-Cell (PIC) techniques.
The module can be used to describe either test particles motion in the fluid electromagnetic field or to solve the fully coupled MHD-PIC system of equations with particle backreaction on the fluid as originally introduced by \cite{Bai_etal.2015}. Particle backreaction on the fluid is included in the form of momentum-energy feedback and by introducing the CR-induced Hall term in Ohm's law. The hybrid MHD-PIC module can be employed to study CR kinetic effects on scales larger than the (ion) skin depth provided the Larmor gyration scale is properly resolved. When applicable, this formulation avoids to resolve microscopic scales offering a substantial computational saving with respect to PIC simulations.
We present a fully-conservative formulation which is second-order accurate in time and space and extends to either Runge-Kutta (RK) or corner-transport-upwind (CTU) time-stepping schemes (for the fluid) while a standard Boris integrator is employed for the particles. For highly-energetic relativistic CRs and in order to overcome the time step restriction a novel sub-cycling strategy that retains second-order accuracy in time is presented. Numerical benchmarks and applications including Bell instability, diffusive shock acceleration and test particle acceleration in reconnecting layers are discussed.
△ Less
Submitted 5 April, 2018;
originally announced April 2018.