Ultra-weak least squares discretizations for unique continuation and Cauchy problems

In this paper, conditional stability estimates are derived for unique continuation and Cauchy problems associated to the Poisson equation in ultra-weak variational form. Numerical approximations are obtained as minima of regularized least squares functionals. The arising dual norms are replaced by discretized dual norms, which leads to a mixed formulation in terms of trial- and test-spaces. For stable pairs of such spaces, and a proper choice of the regularization parameter, the L₂-error on a subdomain in the obtained numerical approximation can be bounded by the best possible fractional power of the sum of the data error and the error of best approximation. Compared to the use of a standard variational formulation, the latter two errors are measured in weaker norms. To avoid the use of C¹-finite element test spaces, nonconforming finite element test spaces can be applied as well. They either lead to the qualitatively same error bound, or in a simplified version, to such an error bound modulo an additional data oscillation term. Numerical results illustrate our theoretical findings.