Stencil scaling for vector-valued PDEs on hybrid grids with applications to generalized Newtonian fluids

08/23/2019
by   Daniel Drzisga, et al.
0

Matrix-free finite element implementations for large applications provide an attractive alternative to standard sparse matrix data formats due to the significantly reduced memory consumption. Here, we show that they are also competitive with respect to the run time in the low order case if combined with suitable stencil scaling techniques. We focus on variable coefficient vector-valued partial differential equations as they arise in many physical applications. The presented method is based on scaling constant reference stencils originating from a linear finite element discretization instead of evaluating the bilinear forms on-the-fly. This method assumes the usage of hierarchical hybrid grids, and it may be applied to vector-valued second-order elliptic partial differential equations directly or as a part of more complicated problems. We provide theoretical and experimental performance estimates showing the advantages of this new approach compared to the traditional on-the-fly integration and stored matrix approaches. In our numerical experiments, we consider two specific mathematical models. Namely, linear elastostatics and incompressible Stokes flow. The final example considers a non-linear shear-thinning generalized Newtonian fluid. For this type of non-linearity, we present an efficient approach to compute a regularized strain rate which is then used to define the node-wise viscosity. Depending on the compute architecture, we could observe maximum speedups of 64 and 122 involved solving a Stokes problem with 12288 compute cores on the state of the art supercomputer SuperMUC-NG.

READ FULL TEXT
research
08/23/2019

Stencil scaling for vector-valued PDEs with applications to generalized Newtonian fluids

Matrix-free finite element implementations for large applications provid...
research
09/20/2017

A stencil scaling approach for accelerating matrix-free finite element implementations

We present a novel approach to fast on-the-fly low order finite element ...
research
10/27/2022

A matrix-free ILU realization based on surrogates

Matrix-free techniques play an increasingly important role in large-scal...
research
09/27/2019

Error estimation for second-order PDEs in non-variational form

Second-order partial differential equations in non-divergence form are c...
research
07/14/2016

Finite Element Integration with Quadrature on the GPU

We present a novel, quadrature-based finite element integration method f...
research
06/23/2023

On particular solutions of linear partial differential equations with polynomial right-hand-sides

This paper introduces general methodologies for constructing closed-form...
research
05/21/2019

Parallel memory-efficient all-at-once algorithms for the sparse matrix triple products in multigrid methods

Multilevel/multigrid methods is one of the most popular approaches for s...

Please sign up or login with your details

Forgot password? Click here to reset