Resistive and multi-fluid RMHD on graphics processing units
Resistive and multi-fluid RMHD on graphics processing units
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 reducing the demand set by these kinds of problems. We report parallel speed-ups of over 21× using double-precision floating-point accuracy, and highlight the optimization strategies required for these schemes, and how they differ from ideal RMHD models. The impact of these results is discussed in the context of the next-generation simulations of neutron star mergers.
magnetohydrodynamics (MHD), methods: numerical, plasmas, relativistic processes
1-12
Wright, A.J.
4960f51d-7e48-4b59-91d9-359af6d559c1
Hawke, I.
fc964672-c794-4260-a972-eaf818e7c9f4
January 2019
Wright, A.J.
4960f51d-7e48-4b59-91d9-359af6d559c1
Hawke, I.
fc964672-c794-4260-a972-eaf818e7c9f4
Wright, A.J. and Hawke, I.
(2019)
Resistive and multi-fluid RMHD on graphics processing units.
Astrophysical Journal, Supplement Series, 240 (1), , [8].
(doi:10.3847/1538-4365/aaf1b0).
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 reducing the demand set by these kinds of problems. We report parallel speed-ups of over 21× using double-precision floating-point accuracy, and highlight the optimization strategies required for these schemes, and how they differ from ideal RMHD models. The impact of these results is discussed in the context of the next-generation simulations of neutron star mergers.
Text
Wright_2019_ApJS_240_8
- Version of Record
More information
Accepted/In Press date: 15 November 2018
e-pub ahead of print date: 4 January 2019
Published date: January 2019
Keywords:
magnetohydrodynamics (MHD), methods: numerical, plasmas, relativistic processes
Identifiers
Local EPrints ID: 429147
URI: http://eprints.soton.ac.uk/id/eprint/429147
ISSN: 0067-0049
PURE UUID: b81ce379-c939-46c8-83cf-f391109be739
Catalogue record
Date deposited: 22 Mar 2019 17:30
Last modified: 18 Mar 2024 03:01
Export record
Altmetrics
Contributors
Author:
A.J. Wright
Download statistics
Downloads from ePrints over the past year. Other digital versions may also be available to download e.g. from the publisher's website.
View more statistics
