An efficient GPU-based h-adaptation framework via linear trees for the flux reconstruction method

DOI10.1016/J.JCP.2024.112823MaRDI QIDQ6129902zbMATH OpenWikidataFDO

Authors Lai Wang, F. D. Witherden, Antony Jameson

Publication date 17 April 2024

Published in Journal of Computational Physics (Search for Journal in Brave)

Full work available at URL https://arxiv.org/abs/2306.06744

zbMATH Keywords

AMRflux reconstructionGPUlinear octreelinear quadtree

Mathematics Subject Classification ID

Numerical methods for partial differential equations, initial value and time-dependent initial-boundary value problems (65Mxx) Computer aspects of numerical algorithms (65Yxx) Basic methods in fluid mechanics (76Mxx)

Abstract: In this paper, we develop the first entirely graphic processing unit (GPU) based h-adaptive flux reconstruction (FR) method with linear trees. The adaptive solver fully operates on the GPU hardware, using a linear quadtree for two dimensional (2D) problems and a linear octree for three dimensional (3D) problems. We articulate how to efficiently perform tree construction, 2:1 balancing, connectivity query, and how to perform adaptation for the flux reconstruction method on the GPU hardware. As a proof of concept, we apply the adaptive flux reconstruction method to solve the inviscid isentropic vortex propagation problem on 2D and 3D meshes to demonstrate the efficiency of the developed adaptive FR method on a single GPU card. Depending on the computational domain size, acceleration of one or two orders of magnitude can be achieved compared to uniform meshing. The total computational cost of adaption, including tree manipulations, connectivity query and data transfer, compared to that of the numerical solver, is insignificant. It can be less than 2% of the total wall clock time for 3D problems even if we perform adaptation as frequent as every 10 time steps with an explicit 3-stage Runge--Kutta time integrator.

Recommendations

• An h-adaptive implementation of the discontinuous Galerkin method for nonlinear hyperbolic conservation laws on unstructured meshes for graphics processing units
• GPU accelerated cell-based adaptive mesh refinement on unstructured quadrilateral grid
• PyFR: an open source framework for solving advection-diffusion type problems on streaming architectures using the flux reconstruction approach
• Adaptive mesh refinement on graphics processing units for applications in gas dynamics
• \(p\)-adaptive hybridized flux reconstruction schemes

Cited in

(1)

• A freestream preserving algorithm for the flux reconstruction method implemented on tetrahedral element

This page was built for publication: An efficient GPU-based \(h\)-adaptation framework via linear trees for the flux reconstruction method

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q6129902)