Research
In order to bridge the gaps between observation, theory, and prediction, modern science has come to heavily rely on sophisticated computational methods. From simulating star formation to proton transport, from climate change and sea level rise to turbulence and combustion (not to mention turbulent combustion), and from structural design to the design of personalized treatment pathways in state-of-the-art medical practice, complex computer simulations are indispensable, regardless of scale. By studying the mathematics behind computer simulations, I try to bring clarity and insight into otherwise intractable problems like these.
My recent research has resulted in discovering a new class of finite element methods, called DPG* methods, and classifying a large category of related methods, called discrete least-squares methods (DLS methods). In addition, I have written general and portable software for high-order and adaptive simulations and applied my work to several structural, fluid, and wave mechanics problems. I am currently interested in applying my work to problems at the frontiers of science and engineering research by expanding on state-of-the-art techniques in optimal control, uncertainty quantification, and high performance computing.
Numerical Methods
Much of my PhD research, supervised by Leszek Demkowicz, was spent focusing on discontinuous Petrov-Galerkin finite element methods (DPG methods) which have been used extensively in engineering and scientific applications. Since beginning my postdoctoral research, I have been focused on an emerging approach to finite element analysis which has been dubbed the "surrogate matrix methodology." In additional, as part of my role in the EU's Horizon 2020 research and innovation program ExaQUte, my research has also included shape optimization under uncertainty, multilevel Monte Carlo methods, and turbulent inlet generation for large eddy simulation. The following summaries only pertain to my contributions to the field of finite element methods.
DPG* Methods

DPG* methods are dual to DPG methods. In a well-defined sense, DPG, as a methodology, can be viewed as a practical means to solve overdetermined Petrov–Galerkin discretizations of boundary value problems. In a similar way, DPG* delivers a methodology for underdetermined discretizations.
Demkowicz, L., Gopalakrishnan, J., and Keith, B. (2020). The DPG-star method. Comput. Math. Appl. [preprint] [doi]
Surrogate Matrix Methods

It can be shown that the majority of entries in the coefficient matrices generated by many Galerkin finite element methods can be computed by evaluating a small number of smooth so-called "stencil functions" at evenly-spaced points throughout the computational domain. Surrogate matrix methods provide low-cost coefficient matrices for such methods by simply evaluating approximations of these stencil functions. The methodology has been applied to classical finite element methods with lowest order simplicial elements as well as isogeometric methods, which employ high-order B-spline or NURBS bases.
In an ongoing series of papers, surrogate matrix methods have been successfully used in a large number of problems with steady or transient characteristics. This includes Poisson’s equation and p-Laplacian diffusion, membrane vibration, plate bending, Stokes’ flow, linear and nonlinear elasticity, and high frequency wave propagation problems. Using surrogate methods can significantly reduce the time to solution in isogeometric analysis (IGA) at a few million degrees of freedom or in large-scale matrix-free applications with billions of degrees of freedom. This ongoing work is part of Daniel Drzisga's PhD research.
Drzisga, D., Keith, B., and Wohlmuth, B. (2020). The surrogate matrix methodology: Accelerating isogeometric analysis of waves. Comput. Methods Appl. Mech. Engrg., 372:113322 [preprint] [doi]
Drzisga, D., Keith, B., and Wohlmuth, B. (2020). The surrogate matrix methodology: A reference implementation for low-cost assembly in isogeometric analysis. MethodsX 7:100813 [preprint] [doi] [code]
Drzisga, D., Keith, B., and Wohlmuth, B. (2020). The surrogate matrix methodology: Low-cost assembly for isogeometric analysis. Comput. Methods Appl. Mech. Engrg., 361:112776 [preprint] [doi]
Drzisga, D., Keith, B., and Wohlmuth, B. (2019). The surrogate matrix methodology: a priori error estimation. SIAM J. Sci. Comput., 41(6):A3806-A3838 [preprint] [doi]
DLS Methods

From extensive research on DPG methods, the term discrete least-squares finite element method (DLS method) was coined. The discrete linear systems in DLS methods have a very specific algebraic structure which can be exploited to accelerate computation or to reduce round-off error.
Keith, B., Petrides, S., Fuentes, F., and Demkowicz, L. (2017). Discrete least-squares finite element methods. Comput. Methods Appl. Mech. Engrg., 327:226-255. [preprint] [doi]
Goal-Oriented Methods

Goal-oriented methods are tailored for accuracy in a specific output. With goal-oriented methods, great efficiency improvements can be achieved because a globally high-quality solution is often not necessary. Once the output is determined, there are many different ways in which to seek efficiency improvements. One way is through goal-oriented adaptive mesh refinement, which can be rigorously formulated in non-symmetric functional settings and applied to DPG methods.
Keith, B., Vaziri Astaneh, A., and Demkowicz, L. (2019). Goal-oriented adaptive mesh refinement for discontinuous Petrov–Galerkin methods. SIAM J. Numer. Anal., 57(4):1649-1676 [preprint] [doi]
Applications
Solids

In a sequence of two papers, Federico Fuentes and I analyzed four different formulations — denoted strong, mixed, primal, and ultraweak, respectively (see figure; left) — of a common linearized elasticity model. This project culminated in the analysis of a very difficult high material contrast coupled rubber and steel model (see figure; center and right). Here, due to the incompressibility of the rubber, standard displacement-only methods will often break down or ''lock.'' Additionally, the steel in the model is very thin and so, likewise, some alternative formulations will often break down there. Our solution was to develop a special method coupling two different formulations together at the material interface; hence, greatly improving the accuracy possible solely with either individual formulation.
Fuentes, F., Keith, B., Demkowicz, L., and Le Tallec, P. (2017). Coupled variational formulations of linear elasticity and the DPG methodology. J. Comput. Phys., 348:715-731. [preprint] [doi]
Keith, B., Fuentes, F., and Demkowicz, L. (2016). The DPG methodology applied to different variational formulations of linear elasticity. Comput. Methods Appl. Mech. Engrg., 309:579-609. [preprint] [doi]
Fluids

Viscoelastic fluid models are commonly used in engineering to simulate blood and polymer melts. These models are well-known to very challenging both in simulation and in mathematical analysis.
I developed a new DPG finite element method of the Oldroyd-B viscoelastic fluid model which is intrinsically stable throughout a broad range of model parameters. I then studied parameter-dependent adaptive mesh refinement with the method. The simulations for this project were completed using the Camellia software library.
Keith, B., Knechtges, P., Roberts, N.V., Elgeti, S., Behr, M., and Demkowicz, L (2017). An ultraweak DPG method for viscoelastic fluids. J. Non-Newton. Fluid Mech., 247:107-122. [preprint] [doi]
Waves

Perfectly matched layers (PMLs) are a very widespread type of artificial absorbing boundary layer used in numerical methods for wave propagation problems defined on unbounded domains. In order to lay the groundwork for some of my colleagues to begin their dissertation research, Ali Vaziri Astaneh and I developed PMLs for high-order DPG methods for acoustic, elastodynamic, and electromagnetic models.
Vaziri Astaneh, A., Keith, B., and Demkowicz, L. (2019). On perfectly matched layers for discontinuous Petrov–Galerkin methods. Comput. Mech., 63(6):1131-1145 [preprint] [doi]
Mixing patterns

A Lagrangian coherent structure (LCS) is an influential material surface which acts as a skeleton of observed mixing patterns in a dynamical system. During my Master's studies, supervised by George Haller, I analyzed several problems via new variational characterizations of LCSs.
Keith, B. (2014). Lagrangian coherent structures in three-dimensional steady flows. Master’s Thesis, McGill University, Montreal, Quebec, Canada. [link]
Software
I have written a significant amount of software in Fortran, Matlab, Python, and C++. Some of the software projects which I have led or contributed to are described below.
hp2D/hp3D
Throughout my PhD studies, I maintained and contributed to hp2D and hp3D. These two finite element libraries, written in Fortran and developed by Leszek Demkowicz, have complete 2D/3D support for local hierarchical anisotropic h- and p-refinement with one level of hanging nodes and shape functions for all standard elements conforming in each of the canonical de Rahm complex of Hilbert spaces:

The software is not open-source due to limited documentation, but it is well-used by the Electromagnetics and Acoustics Group at the Oden Institute and by many of their collaborators.
ESEAS

The ESEAS (exact sequence elements of all shapes) software library supports a complete set of features for evaluating shape functions in high-order finite element methods written in the Fortran 77 standard. The library is innovative for its support for all standard elements in one, two, and three spatial dimensions for arbitrary polynomial order (see figure above; left). It relies on a minimal set of hierarchical routines to construct all the shape functions. This hierarchical construction was the breakthrough which lead to the first explicit construction of arbitrary order pyramid elements (see figure above; right).
ESEAS is currently being used inside hp2D and hp3D, HOFEM, and Sierra. It is being rewritten and optimized in C++ for the Trilinos and MFEM libraries.
Fuentes, F., Keith, B., Demkowicz, L., and Nagaraj, S. (2015). Orientation embedded high order shape functions for the exact sequence elements of all shapes. Comput. Math. Appl., 70(4):353-458. [preprint] [doi] [code]
Camellia
The Camellia software library is a C++ toolbox developed by Nathan Roberts which uses Sandia's Trilinos library of packages. It is a publicly available software with many tools for rapid implementation of finite element methods including discontinuous Galerkin (DG), discontinuous Petrov-Galerkin (DPG), DPG*, hybridizable discontinuous Galerkin (HDG), and first-order system least-squares methods (FOSLS).