Primal-Dual Mixed Finite Element Methods for the Elliptic Cauchy Problem

• consider primal-dual mixed finite element methods for the solution of the elliptic Cauchy problem, or other related data assimilation problems. The method has a local conservation property. We derive a priori error estimates using known conditional stability estimates and determine the minimal amount of weakly consistent stabilization and Tikhonov regularization that yields optimal convergence for smooth exact solutions. The effect of perturbations in data is also accounted for. A reduced version of the method, obtained by choosing a special stabilization of the dual variable, can be viewed as a variant of the least squares mixed finite element method introduced by Darde, Hannukainen, and Hyvonen in [SIAM T. Numer. Anal., 51 (2013), pp. 2123-2148]. The main difference is that our choice of regularization does not depend on auxiliary parameters, the mesh size being the only asymptotic parameter. Finally, we show that the reduced method can be used for defect correction iteration to determine the solution of the full method. The theory is illustrated by some numerical examples.

Ämnesord

NATURVETENSKAP  -- Matematik -- Beräkningsmatematik (hsv//swe)

NATURAL SCIENCES  -- Mathematics -- Computational Mathematics (hsv//eng)

Nyckelord

inverse problem

elliptic Cauchy problem

mixed finite element method

primal-dual method

stabilized methods

Publikations- och innehållstyp

ref (ämneskategori)

art (ämneskategori)