Adaptive least-squares space-time finite element methods
Abstract
We consider the numerical solution of an abstract operator equation Bu=f by using a least-squares approach. We assume that B: X Y* is an isomorphism, and that A : Y Y* implies a norm in Y, where X and Y are Hilbert spaces. The minimizer of the least-squares functional 12 \, \| Bu-f \|A-12, i.e., the solution of the operator equation, is then characterized by the gradient equation Su=B* A-1f with an elliptic and self-adjoint operator S:=B* A-1 B : X X*. When introducing the adjoint p = A-1(f-Bu) we end up with a saddle point formulation to be solved numerically by using a mixed finite element method. Based on a discrete inf-sup stability condition we derive related a priori error estimates. While the adjoint p is zero by construction, its approximation ph serves as a posteriori error indicator to drive an adaptive scheme when discretized appropriately. While this approach can be applied to rather general equations, here we consider second order linear partial differential equations, including the Poisson equation, the heat equation, and the wave equation, in order to demonstrate its potential, which allows to use almost arbitrary space-time finite element methods for the adaptive solution of time-dependent partial differential equations.
Turn this paper into a lesson
ArcXiv compiles a structured reading guide from this paper's metadata: plain-English importance, contributions, prerequisite concepts, which sections to read first, flashcards, and a quiz. Grounded in the abstract, never invented.