Return to search

High-performance implementation of H(div)-conforming elements for incompressible flows

In this thesis, evaluation of H(div)-conforming finite elements is implemented in a high-performance setting and used to solve the incompressible Navier-Stokes equation, obtaining an exactly point-wise divergence-free velocity field. In particular, the anisotropic Raviart-Thomas tensor-product polynomial space is considered, where the finite element operators are evaluated with quadrature in a matrix-free fashion using sum-factorization on tensor-product meshes. The implementation includes evaluation over elements and faces in two- and three-dimensional space, supporting non-conforming meshes with hanging nodes, and using the contravariant Piola transformation to preserve normal components on element boundaries. In terms of throughput, the implementation achieves up to an order of magnitude faster evaluation of finite element operators compared to a matrix-based evaluation. Correctness is demonstrated with optimal convergence rates for various polynomial degrees, as well as exactly divergence-free solutions for the velocity field.

Identiferoai:union.ndltd.org:UPSALLA1/oai:DiVA.org:uu-478600
Date January 2022
CreatorsWik, Niklas
PublisherUppsala universitet, Avdelningen för beräkningsvetenskap
Source SetsDiVA Archive at Upsalla University
LanguageEnglish
Detected LanguageEnglish
TypeStudent thesis, info:eu-repo/semantics/bachelorThesis, text
Formatapplication/pdf
Rightsinfo:eu-repo/semantics/openAccess
RelationUPTEC F, 1401-5757 ; 22027

Page generated in 0.0018 seconds