High-order matrix-free incompressible flow solvers with GPU acceleration and low-order refined preconditioners

10/07/2019
by   Michael Franco, et al.
0

We present a matrix-free flow solver for high-order finite element discretizations of the incompressible Navier-Stokes and Stokes equations with GPU acceleration. For high polynomial degrees, assembling the matrix for the linear systems resulting from the finite element discretization can be prohibitively expensive, both in terms of computational complexity and memory. For this reason, it is necessary to develop matrix-free operators and preconditioners, which can be used to efficiently solve these linear systems without access to the matrix entries themselves. The matrix-free operator evaluations utilize GPU-accelerated sum-factorization techniques to minimize memory movement and maximize throughput. The preconditioners developed in this work are based on a low-order refined methodology with parallel subspace corrections. The saddle-point Stokes system is solved using block-preconditioning techniques, which are robust in mesh size, polynomial degree, time step, and viscosity. For the incompressible Navier-Stokes equations, we make use of projection (fractional step) methods, which require Helmholtz and Poisson solves at each time step. The performance of our flow solvers is assessed on several benchmark problems in two and three spatial dimensions.

READ FULL TEXT
research
05/31/2022

Matrix-free Monolithic Multigrid Methods for Stokes and Generalized Stokes Problems

We consider the widely used continuous 𝒬_k-𝒬_k-1 quadrilateral or hexahe...
research
04/24/2023

Matrix-free GPU-accelerated saddle-point solvers for high-order problems in H(div)

This work describes the development of matrix-free GPU-accelerated solve...
research
10/21/2022

End-to-end GPU acceleration of low-order-refined preconditioning for high-order finite element discretizations

In this paper, we present algorithms and implementations for the end-to-...
research
09/23/2020

Initial Guesses for Sequences of Linear Systems in a GPU-Accelerated Incompressible Flow Solver

We consider several methods for generating initial guesses when iterativ...
research
08/18/2023

Performant low-order matrix-free finite element kernels on GPU architectures

Numerical methods such as the Finite Element Method (FEM) have been succ...
research
08/16/2021

A fast direct solver for integral equations on locally refined boundary discretizations and its application to multiphase flow simulations

In transient simulations of particulate Stokes flow, to accurately captu...
research
11/06/2015

A quantitative performance analysis for Stokes solvers at the extreme scale

This article presents a systematic quantitative performance analysis for...

Please sign up or login with your details

Forgot password? Click here to reset