Energy Correction Schemes in Space and Time

Barbara Wohlmuth*

In this talk, we address the role of the energy within numerical schemes for partial differential equations (pdes). We focus on two dierent aspects: stability of time integration and a priori estimates for spatial discretizations.\newline Lagrange multiplier based space discretizations for time dependent constrained pdes result quite often in DAE systems of index three. Thus typically spurious oscillations can be observed in case of standard time integrations schemes. Here we show how simple variational crimes based on suitable quadrature rules can reduce the index from three to one and yield stable energy preserving time integration schemes. As model example we consider a dynamic contact problem with Coulomb friction and illustrate numerically the influence of a local mass redistribution.\newline In the second part, we focus on stationary elliptic problems where the solution is globally non-smooth. Typical examples are domains with re-entrant corners, heterogeneous coefficients and right hand sides being not in the $H^1$-dual space. These situations occur quite often in the mathematical modelling of multi-physics applications. Standard remedies to handle these type of problems are graded meshes or enrichment, both techniques result in extra implementational work and computational cost. In the case of point sources, we provide optimality results and show that globally no pollution occurs. Although the solution is globally not in $H^1$, we observe on a sequence of uniformly refined meshes optimal $L^2$-norm a priori convergence on subdomains excluding the point sources. In the case of heterogeneous coefficients or re-entrant corners, the situation is more complex. Then the well-known pollution effect is observed, and the convergence can be extremely poor. This also holds true not only for the $L^2$ norm on subdomains excluding the cross-points and corners but also for other quantities of interests such as the stress intensity factors or eigenvalues. Here we introduce a purely local energy correction function and modify locally the bilinear form. Numerical examples in 2D illustrate that we can recover full optimality in case of uniform meshes for re-entrant corners, heterogeneous coeffcients and linear elasticity. As quantities of interest we select the convergence of eigenvalues, the ux across an interface and the stress intensity factor.\newline Partly joint work with M. Melenk (TU Wien), H. Egger (TU Darmstadt), U. Rüde (FAU) and with F. Benesch, T. Horger, M. Huber, T. Köppl, C. Waluga (TU München)

Mathematics Subject Classification:


Plenary Lecture