A GPU-enabled implicit finite volume solver for the ideal two-fluid plasma model on unstructured grids
From MaRDI portal
Publication:6100271
DOI10.1016/J.CPC.2019.01.019arXiv1809.02532OpenAlexW2890919586MaRDI QIDQ6100271FDOQ6100271
Authors: Isaac Alonso Asensio, Alejandro Alvarez Laguna, Mohamed Hassanine Aissa, S. Poedts, Nataly Ozak, Andrea Lani 
Publication date: 12 May 2023
Published in: Computer Physics Communications (Search for Journal in Brave)
Abstract: This paper describes the main features of a pioneering unsteady solver for simulating ideal two-fluid plasmas on unstructured grids, taking profit of GPGPU (General-purpose computing on graphics processing units). The code, which has been implemented within the open source COOLFluiD platform, is implicit, second-order in time and space, relying upon a Finite Volume method for the spatial discretization and a three-point backward Euler for the time integration. In particular, the convective fluxes are computed by a multi-fluid version of the AUSM+up scheme for the plasma equations, in combination with a modified Rusanov scheme with tunable dissipation for the Maxwell equations. Source terms are integrated with a one-point rule, using the cell-centered value. Some critical aspects of the porting to GPU's are discussed, as well as the performance of two open source linear system solvers (i.e. PETSc, PARALUTION). The code design allows for computing both flux and source terms on the GPU along with their Jacobian, giving a noticeable decrease in the computational time in comparison with the original CPU-based solver. The code has been tested in a wide range of mesh sizes and in three different systems, each one with a different GPU. The increased performance (up to 14x) is demonstrated in two representative 2D benchmarks: propagation of circularly polarized waves and the more challenging Geospace Environmental Modeling (GEM) magnetic reconnection challenge.
Full work available at URL: https://arxiv.org/abs/1809.02532
Cites Work
- A sequel to AUSM II: AUSM\(^+\)-up for all speeds
- A three-dimensional finite-volume solver for the Maxwell equations with divergence cleaning on unstructured meshes
- General purpose molecular dynamics simulations fully implemented on graphics processing units
- An object-oriented implementation of a parallel Monte Carlo code for radiation transport
- A GPU-enabled finite volume solver for global magnetospheric simulations on unstructured grids
- A finite volume implicit time integration method for solving the equations of ideal magnetohydrodynamics for the hyperbolic divergence cleaning approach
- Computational Science – ICCS 2005
- A discontinuous Galerkin method for ideal two-fluid plasma equations
- Toward a GPU-aware comparison of explicit and implicit CFD simulations on structured meshes
- Fully-implicit finite volume method for the ideal two-fluid plasma model
- A fully-implicit finite-volume method for multi-fluid reactive and collisional magnetized plasmas on unstructured meshes
- Divergence-free approximate Riemann solver for the quasi-neutral two-fluid plasma model
- Numerical Methods for Ordinary Differential Equations
Cited In (5)
- An asymptotic preserving well-balanced scheme for the isothermal fluid equations in low-temperature plasmas at low-pressure
- A fully-implicit finite-volume method for multi-fluid reactive and collisional magnetized plasmas on unstructured meshes
- Development of an unstructured mesh gyrokinetic particle-in-cell code for exascale fusion plasma simulations on GPUs
- A simple GPU-accelerated two-dimensional MUSCL-Hancock solver for ideal magnetohydrodynamics
- Development of an implicit high-order flux reconstruction solver for the langtry-menter laminar-turbulent transition RANS model
This page was built for publication: A GPU-enabled implicit finite volume solver for the ideal two-fluid plasma model on unstructured grids
Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q6100271)
