|
In this talk, we discuss the least-squares discretization of operator equations Bu=f where B:X→Y* is assumed to be invertible when considering suitable Hilbert spaces Xand Y. Choosing a self-adjoint and elliptic operatorA:Y→Y* we consider the minimization of 1/2||Bu−f||A−12, where the minimizer u=B−1f solves the gradient equation B*A−1Bu=B*A−1f. When introducing the adjoint variable p=A−1(f−Bu), we conclude the systemAp+Bu=f,B*p=0 to be solved by some mixed discretization scheme. While for the solution (u,p)∈X×Y we conclude Bu=f and p≡0, we can use ph∈Yh as a posteriori error indicator to drive adaptive schemes.
The abstract theory is given in our recent publication [2], where we discuss the elliptic diffusion equation, the parabolic heat equation, and the hyperbolic wave equation. In addition, [3] covers the case of the advection-diffusion equation. A more prominent example is the pressure recovery in fluid mechanics, when the velocity is known [4]. Other applications involve boundary element methods in the case of less regular data [5], and least-squares space-time boundary element methods for the wave equation [1].
References
[1] D. Hoonhout, R. Löscher, O. Steinbach, C. Urzua-Torres: Stable least-squares space-time boundary element methods for the wave equation.Adv. Comput. Math.52, 7 (2026).
[2] C. Köthe, R. Löscher, O. Steinbach: Adaptive least-squares space-time finite element methods.Electron. Trans. Numer. Anal., to appear, 2026.
[3] C. Köthe, O. Steinbach: Adaptive least-squares space-time finite element methods for convection-diffusion problems.Comput. Methods Appl. Math.26(2026) 211-239.
[4] D. R. Q. Pacheco, O. Steinbach: Optimal pressure recovery using an ultra-weak finite element method for the pressure Poisson equation and a least-squares approach for the gradient equation.Comput. Methods Appl. Math.24(2024) 921-934.
[5] O. Steinbach: An adaptive least squares boundary element method for elliptic boundary value problems.Numer. Math.158(2026) 653-669. |