A stabilized finite element method for inverse problems subject to the convection-diffusion equation. I: diffusion-dominated regime

• Published in: arXiv.org
• Main Authors: Burman, Erik, Nechita, Mihai, Oksanen, Lauri
• Format: Paper Journal Article
• Language: English
• Published: Ithaca Cornell University Library, arXiv.org 05.10.2019
• Subjects: Approximation; Computer Science - Numerical Analysis; Convection-diffusion equation; Dependence; Diffusion; Divergence; Estimates; Finite element analysis; Finite element method; Inverse problems; Mathematical analysis; Mathematics - Analysis of PDEs; Mathematics - Numerical Analysis; Nonlinear programming; Norms; Numerical analysis; Numerical stability; Upper bounds
• ISSN: 2331-8422
• DOI: 10.48550/arxiv.1811.00431

The numerical approximation of an inverse problem subject to the convection--diffusion equation when diffusion dominates is studied. We derive Carleman estimates that are on a form suitable for use in numerical analysis and with explicit dependence on the Péclet number. A stabilized finite element method is then proposed and analysed. An upper bound on the condition number is first derived. Combining the stability estimates on the continuous problem with the numerical stability of the method, we then obtain error estimates in local \(H^1\)- or \(L^2\)-norms that are optimal with respect to the approximation order, the problem's stability and perturbations in data. The convergence order is the same for both norms, but the \(H^1\)-estimate requires an additional divergence assumption for the convective field. The theory is illustrated in some computational examples.