Source author record

Andrew Gillette

Andrew Gillette appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

12works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

12 published item(s)

preprint2016arXiv

Construction of scalar and vector finite element families on polygonal and polyhedral meshes

We combine theoretical results from polytope domain meshing, generalized barycentric coordinates, and finite element exterior calculus to construct scalar- and vector-valued basis functions for conforming finite element methods on generic convex polytope meshes in dimensions 2 and 3. Our construction recovers well-known bases for the lowest order Nédélec, Raviart-Thomas, and Brezzi-Douglas-Marini elements on simplicial meshes and generalizes the notion of Whitney forms to non-simplicial convex polygons and polyhedra. We show that our basis functions lie in the correct function space with regards to global continuity and that they reproduce the requisite polynomial differential forms described by finite element exterior calculus. We present a method to count the number of basis functions required to ensure these two key properties.

preprint2016arXiv

Finite Element Exterior Calculus for Evolution Problems

Arnold, Falk, and Winther [Bull. Amer. Math. Soc. 47 (2010), 281--354] showed that mixed variational problems, and their numerical approximation by mixed methods, could be most completely understood using the ideas and tools of Hilbert complexes. This led to the development of the Finite Element Exterior Calculus (FEEC) for a large class of linear elliptic problems. More recently, Holst and Stern [Found. Comp. Math. 12:3 (2012), 263--293 and 363--387] extended the FEEC framework to semi-linear problems, and to problems containing variational crimes, allowing for the analysis and numerical approximation of linear and nonlinear geometric elliptic partial differential equations on Riemannian manifolds of arbitrary spatial dimension, generalizing surface finite element approximation theory. In this article, we develop another distinct extension to the FEEC, namely to parabolic and hyperbolic evolution systems, allowing for the treatment of geometric and other evolution problems. Our approach is to combine the recent work on the FEEC for elliptic problems with a classical approach to solving evolution problems via semi-discrete finite element methods, by viewing solutions to the evolution problem as lying in time-parameterized Hilbert spaces (or Bochner spaces). Building on classical approaches by Thomee for parabolic problems and Geveci for hyperbolic problems, we establish a priori error estimates for Galerkin FEM approximation in the natural parametrized Hilbert space norms. In particular, we recover the results of Thomee and Geveci for two-dimensional domains and lowest-order mixed methods as special cases, effectively extending their results to arbitrary spatial dimension and to an entire family of mixed methods. We also show how the Holst and Stern framework allows for extensions of these results to certain semi-linear evolution problems.

preprint2016arXiv

Serendipity and Tensor Product Affine Pyramid Finite Elements

Using the language of finite element exterior calculus, we define two families of $H^1$-conforming finite element spaces over pyramids with a parallelogram base. The first family has matching polynomial traces with tensor product elements on the base while the second has matching polynomial traces with serendipity elements on the base. The second family is new to the literature and provides a robust approach for linking between Lagrange elements on tetrahedra and serendipity elements on affinely-mapped cubes while preserving continuity and approximation properties. We define shape functions and degrees of freedom for each family and prove unisolvence and polynomial reproduction results.

preprint2015arXiv

Constructions of some minimal finite element systems

Within the framework of finite element systems, we show how spaces of differential forms may be constructed, in such a way that they are equipped with commuting interpolators and contain prescribed functions, and are minimal under these constraints. We show how various known mixed finite element spaces fulfill such a design principle, including trimmed polynomial differential forms, serendipity elements and TNT elements. We also comment on virtual element methods and provide a dimension formula for minimal compatible finite element systems containing polynomials of a given degree on hypercubes.

preprint2015arXiv

Interpolation Error Estimates for Harmonic Coordinates On Polytopes

Interpolation error estimates in terms of geometric quality measures are established for harmonic coordinates on polytopes in two and three dimensions. First we derive interpolation error estimates over convex polygons that depend on the geometric quality of the triangles in the constrained Delaunay triangulation of the polygon. This characterization is sharp in the sense that families of polygons with poor quality triangles in their constrained Delaunay triangulations are shown to produce large error when interpolating a basic quadratic function. Non-convex polygons exhibit a similar limitation: large constrained Delaunay triangles caused by vertices approaching a non-adjacent edge also lead to large interpolation error. While this relationship is generalized to convex polyhedra in three dimensions, the possibility of sliver tetrahedra in the constrained Delaunay triangulation prevent the analogous estimate from sharply reflecting the actual interpolation error. Non-convex polyhedra are shown to be fundamentally different through an example of a family of polyhedra containing vertices which are arbitrarily close to non-adjacent faces yet the interpolation error remains bounded.

preprint2014arXiv

Hermite and Bernstein Style Basis Functions for Cubic Serendipity Spaces on Squares and Cubes

We introduce new Hermite-style and Bernstein-style geometric decompositions of the cubic order serendipity finite element spaces $S_3(I^2)$ and $S_3(I^3)$, as defined in the recent work of Arnold and Awanou [Found. Comput. Math. 11 (2011), 337--344]. The serendipity spaces are substantially smaller in dimension than the more commonly used bicubic and tricubic Hermite tensor product spaces - 12 instead of 16 for the square and 32 instead of 64 for the cube - yet are still guaranteed to obtain cubic order \textit{a priori} error estimates in $H^1$ norm when used in finite element methods. The basis functions we define have a canonical relationship both to the finite element degrees of freedom as well as to the geometry of their graphs; this means the bases may be suitable for applications employing isogeometric analysis where domain geometry and functions supported on the domain are described by the same basis functions. Moreover, the basis functions are linear combinations of the commonly used bicubic and tricubic polynomial Bernstein or Hermite basis functions, allowing their rapid incorporation into existing finite element codes.

preprint2013arXiv

Gradient bounds for Wachspress coordinates on polytopes

We derive upper and lower bounds on the gradients of Wachspress coordinates defined over any simple convex d-dimensional polytope P. The bounds are in terms of a single geometric quantity h_*, which denotes the minimum distance between a vertex of P and any hyperplane containing a non-incident face. We prove that the upper bound is sharp for d=2 and analyze the bounds in the special cases of hypercubes and simplices. Additionally, we provide an implementation of the Wachspress coordinates on convex polyhedra using Matlab and employ them in a 3D finite element solution of the Poisson equation on a non-trivial polyhedral mesh. As expected from the upper bound derivation, the H^1-norm of the error in the method converges at a linear rate with respect to the size of the mesh elements.

preprint2012arXiv

Interpolation Error Estimates for Mean Value Coordinates over Convex Polygons

In a similar fashion to estimates shown for Harmonic, Wachspress, and Sibson coordinates in [Gillette et al., AiCM, doi:10.1007/s10444-011-9218-z], we prove interpolation error estimates for the mean value coordinates on convex polygons suitable for standard finite element analysis. Our analysis is based on providing a uniform bound on the gradient of the mean value functions for all convex polygons of diameter one satisfying certain simple geometric restrictions. This work makes rigorous an observed practical advantage of the mean value coordinates: unlike Wachspress coordinates, the gradients of the mean value coordinates do not become large as interior angles of the polygon approach pi.

preprint2012arXiv

Quadratic Serendipity Finite Elements on Polygons Using Generalized Barycentric Coordinates

We introduce a finite element construction for use on the class of convex, planar polygons and show it obtains a quadratic error convergence estimate. On a convex n-gon satisfying simple geometric criteria, our construction produces 2n basis functions, associated in a Lagrange-like fashion to each vertex and each edge midpoint, by transforming and combining a set of n(n+1)/2 basis functions known to obtain quadratic convergence. The technique broadens the scope of the so-called `serendipity' elements, previously studied only for quadrilateral and regular hexahedral meshes, by employing the theory of generalized barycentric coordinates. Uniform `a priori' error estimates are established over the class of convex quadrilaterals with bounded aspect ratio as well as over the class of generic convex planar polygons satisfying additional shape regularity conditions to exclude large interior angles and short edges. Numerical evidence is provided on a trapezoidal quadrilateral mesh, previously not amenable to serendipity constructions, and applications to adaptive meshing are discussed.

preprint2011arXiv

Dual Formulations of Mixed Finite Element Methods with Applications

Mixed finite element methods solve a PDE using two or more variables. The theory of Discrete Exterior Calculus explains why the degrees of freedom associated to the different variables should be stored on both primal and dual domain meshes with a discrete Hodge star used to transfer information between the meshes. We show through analysis and examples that the choice of discrete Hodge star is essential to the numerical stability of the method. Additionally, we define interpolation functions and discrete Hodge stars on dual meshes which can be used to create previously unconsidered mixed methods. Examples from magnetostatics and Darcy flow are examined in detail.

preprint2011arXiv

Error Estimates for Generalized Barycentric Interpolation

We prove the optimal convergence estimate for first order interpolants used in finite element methods based on three major approaches for generalizing barycentric interpolation functions to convex planar polygonal domains. The Wachspress approach explicitly constructs rational functions, the Sibson approach uses Voronoi diagrams on the vertices of the polygon to define the functions, and the Harmonic approach defines the functions as the solution of a PDE. We show that given certain conditions on the geometry of the polygon, each of these constructions can obtain the optimal convergence estimate. In particular, we show that the well-known maximum interior angle condition required for interpolants over triangles is still required for Wachspress functions but not for Sibson functions.