Partial differential equations (PDEs) model a wide variety of phenomena, ranging from how an airplane wing deforms in response to turbulence, to how radio waves travel through and around objects, to how black holes generate gravitational waves when they merge. Numerical analysts develop algorithms for simulating these systems by solving PDEs on a computer; these simulations enable engineers and scientists to develop prototypes and to interpret data from sensors. For example, the NSF-funded Nobel-winning detection of gravitational waves would not have been possible without advances in numerical analysis. In recent decades, numerical analysts discovered that ideas from differential geometry, an area of pure mathematics, can be used to develop good algorithms for solving PDEs. In fact, these ideas help not only for geometric problems in fields of study like computer vision and general relativity, but also for fields like electromagnetism that have little to do with geometry. Although applying differential geometry to numerical analysis has been very successful, thus far this link has been explored only for a small number of differential geometry ideas. In this project, the investigators will continue exploring this link, taking more ideas from differential geometry and applying them to develop new numerical algorithms. These algorithms could then be used both in applied areas, by solving PDEs in science and engineering, and in pure areas, by solving PDEs in differential geometry itself. The project will also support the training of graduate student researchers.<br/><br/>This project focuses on problems at the cusp of numerical analysis and differential geometry. It deals specifically with the design of finite element methods for PDEs that involve vector fields and tensor fields on Riemannian manifolds. In the long term, these efforts have the potential to lead to robust numerical methods for solving geometric PDEs like the Einstein field equations, which are useful for studying gravitational wave signals, as well as PDEs like the elasticity equations, which model how objects deform under stress. This project has three main goals. The first is to develop a new family of finite elements for discretizing algebraic curvature tensors and other bi-forms---tensor products of differential forms---on simplicial triangulations. The second goal is to develop an intrinsic finite element discretization of the Bochner Laplacian, which is a basic differential operator in Riemannian geometry that differs from the familiar Hodge Laplacian from finite element exterior calculus. The third goal is to leverage what we learn to design numerical methods for a wide range of geometric problems, such as computing spectra of elliptic operators on manifolds, simulating intrinsic geometric flows, and solving prescribed curvature problems.<br/><br/>This award reflects NSF's statutory mission and has been deemed worthy of support through evaluation using the Foundation's intellectual merit and broader impacts review criteria.