2019
DOI: 10.3847/1538-4365/aaf1b0
|View full text |Cite
|
Sign up to set email alerts
|

Resistive and Multi-fluid RMHD on Graphics Processing Units

Abstract: In this work we present a proof of concept of CUDA-capable, resistive, multi-fluid models of relativistic magnetohydrodynamics (RMHD). Resistive and multi-fluid codes for simulating models of RMHD suffer from stiff source terms, so it is common to implement a set of semi-implicit time integrators to maintain numerical stability. We show, for the first time, that finite volume IMEX schemes for resistive and two-fluid models of RMHD can be accelerated by execution on graphics processing units, significantly redu… Show more

Help me understand this report
View preprint versions

Search citation statements

Order By: Relevance

Paper Sections

Select...
1
1
1
1

Citation Types

0
5
0

Year Published

2019
2019
2024
2024

Publication Types

Select...
3
3

Relationship

1
5

Authors

Journals

citations
Cited by 6 publications
(5 citation statements)
references
References 39 publications
0
5
0
Order By: Relevance
“…All results have been generated using the METHOD code 2 . Details of the numerical schemes used in METHOD can be found in more detail in [13]. Briefly, the numerical flux is computed using flux splitting, with a third order WENO Larger resistivities lead to a greater rate of magnetic diffusion and thus greater smearing out compared to the ideal MHD solution.…”
Section: Methodsmentioning
confidence: 99%
See 1 more Smart Citation
“…All results have been generated using the METHOD code 2 . Details of the numerical schemes used in METHOD can be found in more detail in [13]. Briefly, the numerical flux is computed using flux splitting, with a third order WENO Larger resistivities lead to a greater rate of magnetic diffusion and thus greater smearing out compared to the ideal MHD solution.…”
Section: Methodsmentioning
confidence: 99%
“…One issue with more physically complex models that prevents the wide spread adoption of the resistive GRMHD equations is the additional computation required in their evolution. Due to the possibly stiff source term of resistive GRMHD [2], an implicit integrator, often the IMEX schemes of [10], is required to keep execution times rea-sonable [2,3,11,12,13]. Whilst these schemes allow for a faster evolution than would be possible using traditional explicit schemes for the majority of the parameter space, they still result in at least a factor 5× slow-down compared to conventional, ideal GRMHD models [5].…”
Section: Motivationmentioning
confidence: 99%
“…The restrictions on 𝜎 and 𝑏/𝑎 reduce the sample to 912 galaxies. The parameters 𝑎 deV , 𝑏/𝑎, and 𝜎 fiber in 𝑟-band are taken from the catalogs PhotObjAll and SpecObjAll of DR15 (Aguado et al 2019).…”
Section: 𝑀 * − 𝜎 Relationmentioning
confidence: 99%
“…In RHD and RMHD codes, this task requires solving a highly nonlinear algebraic system of equations, which constitutes one of the most challenging tasks from the point of view of computational efficiency. The recovery of primitive variables is usually considered as one of the bottlenecks for the speed-up of the code (Wright & Hawke 2019). Of the different strategies that exist in the literature (see, e.g, Martí & Müller 2015, and references therein), we chose to recover the set of primitive variables following the inversion scheme of Mignone & McKinney (2007).…”
Section: A22 Spatial Cell Reconstructionmentioning
confidence: 99%