Repository logo
  • English
  • Italiano
Log In
Have you forgotten your password?
  1. Home
  2. PRODOTTI RICERCA INAF
  3. 1 CONTRIBUTI IN RIVISTE (Journal articles)
  4. 1.01 Articoli in rivista
  5. A Constrained Transport Method for the Solution of the Resistive Relativistic MHD Equations
 

A Constrained Transport Method for the Solution of the Resistive Relativistic MHD Equations

Journal
MONTHLY NOTICES OF THE ROYAL ASTRONOMICAL SOCIETY  
Date Issued
2019
Author(s)
Mignone, A.
•
Mattia, G.
•
BODO, Gianluigi  
•
Del Zanna, L.
DOI
10.1093/mnras/stz1015
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 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.
Volume
486
Issue
3
Start page
4252
Uri
http://hdl.handle.net/20.500.12386/29339
Url
https://academic.oup.com/mnras/article/486/3/4252/5479252
http://arxiv.org/abs/1904.01530v1
Issn Identifier
0035-8711
Rights
open.access
File(s)
Loading...
Thumbnail Image
Name

29339-stz1015_P01.pdf

Description
Miur
Size

5.42 MB

Format

Adobe PDF

Checksum (MD5)

f1f1e475e40f0dcac4d2a784b1dde3e8

Loading...
Thumbnail Image
Name

stz1015.pdf

Description
pdf editoriale
Size

10.86 MB

Format

Adobe PDF

Checksum (MD5)

f363ced2e86738cc0598e2dc16a850d6

Explore By
  • Communities and Collection
  • Research Outputs
  • Researchers
  • Organizations
  • Projects
Information and guides for authors
  • https://openaccess-info.inaf.it: all about open access in INAF
  • How to enter a product: guides to OA@INAF
  • The INAF Policy on Open Access
  • Downloadable documents and templates

Built with DSpace-CRIS software - Extension maintained and optimized by 4Science

  • Privacy policy
  • End User Agreement
  • Send Feedback