Max Winkler

Thomas Apel, Ariel L. Lombardi, Max Winkler

The paper is concerned with the finite element solution of the Poisson equation with homogeneous Dirichlet boundary condition in a three-dimensional domain. Anisotropic, graded meshes from a former paper are reused for dealing with the singular behaviour of the solution in the vicinity of the non-smooth parts of the boundary. The discretization error is analyzed for the piecewise linear approximation in the H^1(Ω)- and L^2(Ω)-norms by using a new quasi-interpolation operator. This new interpolant is introduced in order to prove the estimates for L^2(Ω)-data in the differential equation which is not possible for the standard nodal interpolant. These new estimates allow for the extension of certain error estimates for optimal control problems with elliptic partial differential equation and for a simpler proof of the discrete compactness property for edge elements of any order on this kind of finite element meshes.

J. Pfefferer, M. Winkler

This paper is concerned with approximations and related discretization error estimates for the normal derivatives of solutions of linear elliptic partial differential equations. In order to illustrate the ideas, we consider the Poisson equation with homogeneous Dirichlet boundary conditions and use standard linear finite elements for its discretization. The underlying domain is assumed to be polygonal but not necessarily convex. Approximations of the normal derivatives are introduced in a standard way as well as in a variational sense. On general quasi-uniform meshes, one can show that these approximate normal derivatives possess a convergence rate close to one in $L^2$ as long as the singularities due to the corners are mild enough. Using boundary concentrated meshes, we show that the order of convergence can even be doubled in terms of the mesh parameter while increasing the complexity of the discrete problems only by a small factor. As an application, we use these results for the numerical analysis of Dirichlet boundary control problems, where the control variable corresponds to the normal derivative of some adjoint variable.

Thomas Apel, Johannes Pfefferer, Sergejs Rogovs et al.

This paper deals with a priori pointwise error estimates for the finite element solution of boundary value problems with Neumann boundary conditions in polygonal domains. Due to the corners of the domain, the convergence rate of the numerical solutions can be lower than in case of smooth domains. As a remedy the use of local mesh refinement near the corners is investigated. In order to prove quasi-optimal a priori error estimates regularity results in weighted Sobolev spaces are exploited. This is the first work on the Neumann boundary value problem where both the regularity of the data is exactly specified and the sharp convergence order $h^2 \lvert \ln h \rvert $ in the case of piecewise linear finite element approximations is obtained. As an extension we show the same rate for the approximate solution of a semilinear boundary value problem. The proof relies in this case on the supercloseness between the Ritz projection to the continuous solution and the finite element solution.

Jan Blechschmidt, Tom-Christian Riemer, Max Winkler et al.

We develop a novel physics informed deep learning approach for solving nonlinear drift-diffusion equations on metric graphs. These models represent an important model class with a large number of applications in areas ranging from transport in biological cells to the motion of human crowds. While traditional numerical schemes require a large amount of tailoring, especially in the case of model design or parameter identification problems, physics informed deep operator networks (DeepONet) have emerged as a versatile tool for the solution of partial differential equations with the particular advantage that they easily incorporate parameter identification questions. We here present an approach where we first learn three DeepONet models for representative inflow, inner and outflow edges, resp., and then subsequently couple these models for the solution of the drift-diffusion metric graph problem by relying on an edge-based domain decomposition approach. We illustrate that our framework is applicable for the accurate evaluation of graph-coupled physics models and is well suited for solving optimization or inverse problems on these coupled networks.