The University of Southampton
University of Southampton Institutional Repository

Resistive and multi-fluid RMHD on graphics processing units

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
0067-0049
1-12
Wright, A.J.
4960f51d-7e48-4b59-91d9-359af6d559c1
Hawke, I.
fc964672-c794-4260-a972-eaf818e7c9f4
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. The Astrophysical Journal Supplement Series, 240 (1), 1-12, [8]. (doi:10.3847/1538-4365/aaf1b0).

Record type: Article

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
Available under License Creative Commons Attribution.
Download (1MB)

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
ORCID for I. Hawke: ORCID iD orcid.org/0000-0003-4805-0309

Catalogue record

Date deposited: 22 Mar 2019 17:30
Last modified: 20 Jul 2019 00:56

Export record

Altmetrics

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

Atom RSS 1.0 RSS 2.0

Contact ePrints Soton: eprints@soton.ac.uk

ePrints Soton supports OAI 2.0 with a base URL of http://eprints.soton.ac.uk/cgi/oai2

This repository has been built using EPrints software, developed at the University of Southampton, but available to everyone to use.

We use cookies to ensure that we give you the best experience on our website. If you continue without changing your settings, we will assume that you are happy to receive cookies on the University of Southampton website.

×