carrier image

Anisotropic mesh adaptation in Firedrake with PETSc DMPlex

Barral, Nicolas, Matthew G. Knepley, Michael Lange, Matthew D. Piggott, Gerard J. Gorman

Research Note, 25th International Meshing Roundtable, Sandia National Laboratories, September 26-30 2016

INTERNATIONAL
MESHING
ROUNTABLE

25th International Meshing Roundtable
Washington DC, U.S.A.
September 26-30, 2016

Nicolas Barral, Imperial College London, GB, n.barral@imperial.ac.uk
Matthew G. Knepley, Rice University, US, knepley@rice.edu
Michael Lange, Imperial College London, GB, michael.lange@imperial.ac.uk
Matthew D. Piggott, Imperial College London, GB, m.d.piggott@imperial.ac.uk
Gerard J. Gorman, Imperial College London, GB, g.gorman@imperial.ac.uk

Research Note Abstract
Despite decades of research in this area, mesh adaptation capabilities are still rarely found in numerical simulation software. We postulate that the primary reason for this is lack of usability. Integrating mesh adaptation into existing software is difficult as non-trivial operators, such as error metrics and interpolation operators, are required, and integrating available adaptive remeshers is not straightforward. Our approach presented here is to first integrate Pragmatic, an anisotropic mesh adaptation library, into DMPlex, a PETSc object that manages unstructured meshes and their interactions with PETSc’s solvers and I/O routines. As PETSc is already widely used, this will make anisotropic mesh adaptation available to a much larger community. As a demonstration of this we describe the integration of anisotropic mesh adaptation into Firedrake, an automated Finite Element based system for the portable solution of partial differential equations which already uses PETSc solvers and I/O via DMPlex. We present a proof of concept of this integration with a three-dimensional advection test case.

Download Full Paper (PDF)


Contact author(s) or publisher for availability and copyright information on above referenced article