-
A multifidelity Bayesian optimization method for inertial confinement fusion design
Authors:
J. Wang,
N. Chiang,
A. Gillette,
J. L. Peterson
Abstract:
Due to their cost, experiments for inertial confinement fusion (ICF) heavily rely on numerical simulations to guide design. As simulation technology progresses, so too can the fidelity of models used to plan for new experiments. However, these high-fidelity models are by themselves insufficient for optimal experimental design, because their computational cost remains too high to efficiently and ef…
▽ More
Due to their cost, experiments for inertial confinement fusion (ICF) heavily rely on numerical simulations to guide design. As simulation technology progresses, so too can the fidelity of models used to plan for new experiments. However, these high-fidelity models are by themselves insufficient for optimal experimental design, because their computational cost remains too high to efficiently and effectively explore the numerous parameters required to describe a typical experiment. Traditionally, ICF design has relied on low-fidelity modeling to initially identify potentially interesting design regions, which are then subsequently explored via selected high-fidelity modeling. In this paper, we demonstrate that this two-step approach can be insufficient: even for simple design problems, a two-step optimization strategy can lead high-fidelity searching towards incorrect regions and consequently waste computational resources on parameter regimes far away from the true optimal solution. We reveal that a primary cause of this behavior in ICF design problems is the presence of low-fidelity optima in distinct regions of the parameter space from high-fidelity optima. To address this issue, we propose an iterative multifidelity Bayesian optimization method based on Gaussian Process Regression that leverages both low- and high-fidelity modelings. We demonstrate, using both two- and eight-dimensional ICF test problems, that our algorithm can effectively utilize low-fidelity modeling for exploration, while automatically refining promising designs with high-fidelity models. This approach proves to be more efficient than relying solely on high-fidelity modeling for optimization.
△ Less
Submitted 15 December, 2023;
originally announced December 2023.
-
Learning robust marking policies for adaptive mesh refinement
Authors:
Andrew Gillette,
Brendan Keith,
Socratis Petrides
Abstract:
In this work, we revisit the marking decisions made in the standard adaptive finite element method (AFEM). Experience shows that a naïve marking policy leads to inefficient use of computational resources for adaptive mesh refinement (AMR). Consequently, using AFEM in practice often involves ad-hoc or time-consuming offline parameter tuning to set appropriate parameters for the marking subroutine.…
▽ More
In this work, we revisit the marking decisions made in the standard adaptive finite element method (AFEM). Experience shows that a naïve marking policy leads to inefficient use of computational resources for adaptive mesh refinement (AMR). Consequently, using AFEM in practice often involves ad-hoc or time-consuming offline parameter tuning to set appropriate parameters for the marking subroutine. To address these practical concerns, we recast AMR as a Markov decision process in which refinement parameters can be selected on-the-fly at run time, without the need for pre-tuning by expert users. In this new paradigm, the refinement parameters are also chosen adaptively via a marking policy that can be optimized using methods from reinforcement learning. We use the Poisson equation to demonstrate our techniques on $h$- and $hp$-refinement benchmark problems, and our experiments suggest that superior marking policies remain undiscovered for many classical AFEM applications. Furthermore, an unexpected observation from this work is that marking policies trained on one family of PDEs are sometimes robust enough to perform well on problems far outside the training family. For illustration, we show that a simple $hp$-refinement policy trained on 2D domains with only a single re-entrant corner can be deployed on far more complicated 2D domains, and even 3D domains, without significant performance loss. For reproduction and broader adoption, we accompany this work with an open-source implementation of our methods.
△ Less
Submitted 26 December, 2023; v1 submitted 13 July, 2022;
originally announced July 2022.
-
Data-driven geometric scale detection via Delaunay interpolation
Authors:
Andrew Gillette,
Eugene Kur
Abstract:
Accurate approximation of a real-valued function depends on two aspects of the available data: the density of inputs within the domain of interest and the variation of the outputs over that domain. There are few methods for assessing whether the density of inputs is \textit{sufficient} to identify the relevant variations in outputs -- i.e., the "geometric scale" of the function -- despite the fact…
▽ More
Accurate approximation of a real-valued function depends on two aspects of the available data: the density of inputs within the domain of interest and the variation of the outputs over that domain. There are few methods for assessing whether the density of inputs is \textit{sufficient} to identify the relevant variations in outputs -- i.e., the "geometric scale" of the function -- despite the fact that sampling density is closely tied to the success or failure of an approximation method. In this paper, we introduce a general purpose, computational approach to detecting the geometric scale of real-valued functions over a fixed domain using a deterministic interpolation technique from computational geometry. Our algorithm is based on the observation that a sequence of piecewise linear interpolants will converge to a continuous function at a quadratic rate (in $L^2$ norm) if and only if the data are sampled densely enough to distinguish the feature from noise. We present numerical experiments demonstrating how our method can identify feature scale, estimate uncertainty in feature scale, and assess the sampling density for fixed (i.e. static) datasets of input-output pairs. In addition, we include analytical results in support of our numerical findings and will release lightweight code that can be adapted for use in a variety of data science settings.
△ Less
Submitted 10 March, 2022;
originally announced March 2022.
-
Bringing Trimmed Serendipity Methods to Computational Practice in Firedrake
Authors:
Justin Crum,
Cyrus Cheng,
David A. Ham,
Lawrence Mitchell,
Robert C. Kirby,
Joshua A. Levine,
Andrew Gillette
Abstract:
We present an implementation of the trimmed serendipity finite element family, using the open source finite element package Firedrake. The new elements can be used seamlessly within the software suite for problems requiring $H^1$, \hcurl, or \hdiv-conforming elements on meshes of squares or cubes. To test how well trimmed serendipity elements perform in comparison to traditional tensor product ele…
▽ More
We present an implementation of the trimmed serendipity finite element family, using the open source finite element package Firedrake. The new elements can be used seamlessly within the software suite for problems requiring $H^1$, \hcurl, or \hdiv-conforming elements on meshes of squares or cubes. To test how well trimmed serendipity elements perform in comparison to traditional tensor product elements, we perform a sequence of numerical experiments including the primal Poisson, mixed Poisson, and Maxwell cavity eigenvalue problems. Overall, we find that the trimmed serendipity elements converge, as expected, at the same rate as the respective tensor product elements while being able to offer significant savings in the time or memory required to solve certain problems.
△ Less
Submitted 8 October, 2021; v1 submitted 27 April, 2021;
originally announced April 2021.
-
Extending discrete exterior calculus to a fractional derivative
Authors:
Justin Crum,
Joshua A. Levine,
Andrew Gillette
Abstract:
Fractional partial differential equations (FDEs) are used to describe phenomena that involve a "non-local" or "long-range" interaction of some kind. Accurate and practical numerical approximation of their solutions is challenging due to the dense matrices arising from standard discretization procedures. In this paper, we begin to extend the well-established computational toolkit of Discrete Exteri…
▽ More
Fractional partial differential equations (FDEs) are used to describe phenomena that involve a "non-local" or "long-range" interaction of some kind. Accurate and practical numerical approximation of their solutions is challenging due to the dense matrices arising from standard discretization procedures. In this paper, we begin to extend the well-established computational toolkit of Discrete Exterior Calculus (DEC) to the fractional setting, focusing on proper discretization of the fractional derivative. We define a Caputo-like fractional discrete derivative, in terms of the standard discrete exterior derivative operator from DEC, weighted by a measure of distance between $p$-simplices in a simplicial complex. We discuss key theoretical properties of the fractional discrete derivative and compare it to the continuous fractional derivative via a series of numerical experiments.
△ Less
Submitted 16 July, 2019; v1 submitted 2 May, 2019;
originally announced May 2019.
-
Computational Serendipity and Tensor Product Finite Element Differential Forms
Authors:
Andrew Gillette,
Tyler Kloefkorn,
Victoria Sanders
Abstract:
Many conforming finite elements on squares and cubes are elegantly classified into families by the language of finite element exterior calculus and presented in the Periodic Table of the Finite Elements. Use of these elements varies, based principally on the ease or difficulty in finding a "computational basis" of shape functions for element families. The tensor product family, $Q^-_rΛ^k$, is most…
▽ More
Many conforming finite elements on squares and cubes are elegantly classified into families by the language of finite element exterior calculus and presented in the Periodic Table of the Finite Elements. Use of these elements varies, based principally on the ease or difficulty in finding a "computational basis" of shape functions for element families. The tensor product family, $Q^-_rΛ^k$, is most commonly used because computational basis functions are easy to state and implement. The trimmed and non-trimmed serendipity families, $S^-_rΛ^k$ and $S_rΛ^k$ respectively, are used less frequently because they are newer to the community and, until now, lacked a straightforward technique for computational basis construction. This represents a missed opportunity for computational efficiency as the serendipity elements in general have fewer degrees of freedom than elements of equivalent accuracy from the tensor product family. Accordingly, in pursuit of easy adoption of the serendipity families, we present complete lists of computational bases for both serendipity families, for any order $r\geq 1$ and for any differential form order $0\leq k\leq n$, for problems in dimension $n=2$ or $3$. The bases are defined via shared subspace structures, allowing easy comparison of elements across families. We use and include code in SageMath to find, list, and verify these computational basis functions.
△ Less
Submitted 31 May, 2018;
originally announced June 2018.
-
Nonstandard finite element de Rham complexes on cubical meshes
Authors:
Andrew Gillette,
Kaibo Hu,
Shuo Zhang
Abstract:
We propose two general operations on finite element differential complexes on cubical meshes that can be used to construct and analyze sequences of "nonstandard" convergent methods. The first operation, called DoF-transfer, moves edge degrees of freedom to vertices in a way that reduces global degrees of freedom while increasing continuity order at vertices. The second operation, called serendipit…
▽ More
We propose two general operations on finite element differential complexes on cubical meshes that can be used to construct and analyze sequences of "nonstandard" convergent methods. The first operation, called DoF-transfer, moves edge degrees of freedom to vertices in a way that reduces global degrees of freedom while increasing continuity order at vertices. The second operation, called serendipity, eliminates interior bubble functions and degrees of freedom locally on each element without affecting edge degrees of freedom. These operations can be used independently or in tandem to create nonstandard complexes that incorporate Hermite, Adini and "trimmed-Adini" elements. We show that the resulting elements provide convergent, non-conforming methods for problems requiring stronger regularity and satisfy a discrete Korn inequality. We discuss potential benefits of applying these elements to Stokes, biharmonic and elasticity problems.
△ Less
Submitted 12 April, 2018;
originally announced April 2018.
-
Numerical studies of serendipity and tensor product elements for eigenvalue problems
Authors:
Andrew Gillette,
Craig Gross,
Ken Plackowski
Abstract:
While the use of finite element methods for the numerical approximation of eigenvalues is a well-studied problem, the use of serendipity elements for this purpose has received little attention in the literature. We show by numerical experiments that serendipity elements, which are defined on a square reference geometry, can attain the same order of accuracy as their tensor product counterparts whi…
▽ More
While the use of finite element methods for the numerical approximation of eigenvalues is a well-studied problem, the use of serendipity elements for this purpose has received little attention in the literature. We show by numerical experiments that serendipity elements, which are defined on a square reference geometry, can attain the same order of accuracy as their tensor product counterparts while using dramatically fewer degrees of freedom. In some cases, the serendipity method uses only half as many basis functions as the tensor product method while still producing the same numerical approximation of an eigenvalue. To encourage the further use and study of serendipity elements, we provide a table of serendipity basis functions for low order cases and a Mathematica file that can be used to generate the basis functions for higher order cases.
△ Less
Submitted 30 November, 2016;
originally announced December 2016.
-
Trimmed Serendipity Finite Element Differential Forms
Authors:
Andrew Gillette,
Tyler Kloefkorn
Abstract:
We introduce the family of trimmed serendipity finite element differential form spaces, defined on cubical meshes in any number of dimensions, for any polynomial degree, and for any form order. The relation between the trimmed serendipity family and the (non-trimmed) serendipity family developed by Arnold and Awanou [Math. Comp. 83(288) 2014] is analogous to the relation between the trimmed and (n…
▽ More
We introduce the family of trimmed serendipity finite element differential form spaces, defined on cubical meshes in any number of dimensions, for any polynomial degree, and for any form order. The relation between the trimmed serendipity family and the (non-trimmed) serendipity family developed by Arnold and Awanou [Math. Comp. 83(288) 2014] is analogous to the relation between the trimmed and (non-trimmed) polynomial finite element differential form families on simplicial meshes from finite element exterior calculus. We provide degrees of freedom in the general setting and prove that they are unisolvent for the trimmed serendipity spaces. The sequence of trimmed serendipity spaces with a fixed polynomial order r provides an explicit example of a system described by Christiansen and Gillette [ESAIM:M2AN 50(3) 2016], namely, a minimal compatible finite element system on squares or cubes containing order r-1 polynomial differential forms.
△ Less
Submitted 5 January, 2018; v1 submitted 2 July, 2016;
originally announced July 2016.
-
Serendipity and Tensor Product Affine Pyramid Finite Elements
Authors:
Andrew Gillette
Abstract:
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 robu…
▽ More
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.
△ Less
Submitted 9 September, 2016; v1 submitted 22 December, 2015;
originally announced December 2015.
-
Constructions of some minimal finite element systems
Authors:
Snorre H. Christiansen,
Andrew Gillette
Abstract:
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 el…
▽ More
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.
△ Less
Submitted 20 October, 2015; v1 submitted 17 April, 2015;
originally announced April 2015.
-
Interpolation Error Estimates for Harmonic Coordinates On Polytopes
Authors:
Andrew Gillette,
Alexander Rand
Abstract:
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…
▽ More
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.
△ Less
Submitted 2 October, 2015; v1 submitted 2 April, 2015;
originally announced April 2015.
-
Construction of scalar and vector finite element families on polygonal and polyhedral meshes
Authors:
Andrew Gillette,
Alexander Rand,
Chandrajit Bajaj
Abstract:
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…
▽ More
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.
△ Less
Submitted 25 April, 2016; v1 submitted 27 May, 2014;
originally announced May 2014.
-
Nodal bases for the serendipity family of finite elements
Authors:
Michael S. Floater,
Andrew Gillette
Abstract:
Using the notion of multivariate lower set interpolation, we construct nodal basis functions for the serendipity family of finite elements, of any order and any dimension. For the purpose of computation, we also show how to express these functions as linear combinations of tensor-product polynomials.
Using the notion of multivariate lower set interpolation, we construct nodal basis functions for the serendipity family of finite elements, of any order and any dimension. For the purpose of computation, we also show how to express these functions as linear combinations of tensor-product polynomials.
△ Less
Submitted 15 February, 2016; v1 submitted 24 April, 2014;
originally announced April 2014.
-
Gradient bounds for Wachspress coordinates on polytopes
Authors:
Michael Floater,
Andrew Gillette,
N. Sukumar
Abstract:
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 hyper…
▽ More
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.
△ Less
Submitted 5 December, 2013; v1 submitted 18 June, 2013;
originally announced June 2013.
-
Hermite and Bernstein Style Basis Functions for Cubic Serendipity Spaces on Squares and Cubes
Authors:
Andrew Gillette
Abstract:
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 ins…
▽ More
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.
△ Less
Submitted 12 February, 2014; v1 submitted 29 August, 2012;
originally announced August 2012.
-
Finite Element Exterior Calculus for Evolution Problems
Authors:
Andrew Gillette,
Michael Holst,
Yunrong Zhu
Abstract:
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…
▽ More
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.
△ Less
Submitted 12 October, 2016; v1 submitted 7 February, 2012;
originally announced February 2012.
-
Interpolation Error Estimates for Mean Value Coordinates over Convex Polygons
Authors:
Alexander Rand,
Andrew Gillette,
Chandrajit Bajaj
Abstract:
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 polygo…
▽ More
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.
△ Less
Submitted 18 September, 2012; v1 submitted 23 November, 2011;
originally announced November 2011.
-
Quadratic Serendipity Finite Elements on Polygons Using Generalized Barycentric Coordinates
Authors:
Alexander Rand,
Andrew Gillette,
Chandrajit Bajaj
Abstract:
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…
▽ More
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.
△ Less
Submitted 20 July, 2012; v1 submitted 15 September, 2011;
originally announced September 2011.
-
Dual Formulations of Mixed Finite Element Methods with Applications
Authors:
Andrew Gillette,
Chandrajit Bajaj
Abstract:
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 esse…
▽ More
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.
△ Less
Submitted 11 May, 2011; v1 submitted 17 December, 2010;
originally announced December 2010.
-
Error Estimates for Generalized Barycentric Interpolation
Authors:
Andrew Gillette,
Alexander Rand,
Chandrajit Bajaj
Abstract:
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 Harmo…
▽ More
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.
△ Less
Submitted 15 April, 2011; v1 submitted 24 October, 2010;
originally announced October 2010.
-
Stable Mesh Decimation
Authors:
Chandrajit Bajaj,
Andrew Gillette,
Qin Zhang
Abstract:
Current mesh reduction techniques, while numerous, all primarily reduce mesh size by successive element deletion (e.g. edge collapses) with the goal of geometric and topological feature preservation. The choice of geometric error used to guide the reduction process is chosen independent of the function the end user aims to calculate, analyze, or adaptively refine. In this paper, we argue that su…
▽ More
Current mesh reduction techniques, while numerous, all primarily reduce mesh size by successive element deletion (e.g. edge collapses) with the goal of geometric and topological feature preservation. The choice of geometric error used to guide the reduction process is chosen independent of the function the end user aims to calculate, analyze, or adaptively refine. In this paper, we argue that such a decoupling of structure from function modeling is often unwise as small changes in geometry may cause large changes in the associated function. A stable approach to mesh decimation, therefore, ought to be guided primarily by an analysis of functional sensitivity, a property dependent on both the particular application and the equations used for computation (e.g. integrals, derivatives, or integral/partial differential equations). We present a methodology to elucidate the geometric sensitivity of functionals via two major functional discretization techniques: Galerkin finite element and discrete exterior calculus. A number of examples are given to illustrate the methodology and provide numerical examples to further substantiate our choices.
△ Less
Submitted 8 October, 2009;
originally announced October 2009.