Source author record

Thomas Führer

Thomas Führer appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

19works
3topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

19 published item(s)

preprint2022arXiv

A DPG method for Reissner-Mindlin plates

We present a discontinuous Petrov-Galerkin (DPG) method with optimal test functions for the Reissner-Mindlin plate bending model. Our method is based on a variational formulation that utilizes a Helmholtz decomposition of the shear force. It produces approximations of the primitive variables and the bending moments. For any canonical selection of boundary conditions the method converges quasi-optimally. In the case of hard-clamped convex plates, we prove that the lowest-order scheme is locking free. Several numerical experiments confirm our results.

preprint2022arXiv

Space-time finite element methods for parabolic distributed optimal control problems

We present a method for the numerical approximation of distributed optimal control problems constrained by parabolic partial differential equations. We complement the first-order optimality condition by a recently developed space-time variational formulation of parabolic equations which is coercive in the energy norm, and a Lagrangian multiplier. Our final formulation fulfills the Babuška-Brezzi conditions on the continuous as well as discrete level, without restrictions. Consequently, we can allow for final-time desired states, and obtain an a-posteriori error estimator which is efficient and reliable. Numerical experiments confirm our theoretical findings.

preprint2021arXiv

DPG methods for a fourth-order div problem

We study a fourth-order div problem and its approximation by the discontinuous Petrov-Galerkin method with optimal test functions. We present two variants, based on first and second-order systems. In both cases we prove well-posedness of the formulation and quasi-optimal convergence of the approximation. Our analysis includes the fully-discrete schemes with approximated test functions, for general dimension and polynomial degree in the first-order case, and for two dimensions and lowest-order approximation in the second-order case. Numerical results illustrate the performance for quasi-uniform and adaptively refined meshes.

preprint2020arXiv

A locking-free DPG scheme for Timoshenko beams

We develop a discontinuous Petrov-Galerkin scheme with optimal test functions (DPG method) for the Timoshenko beam bending model with various boundary conditions, combining clamped, supported, and free ends. Our scheme approximates the transverse deflection and bending moment. It converges quasi-optimally in $L_2$ and is locking free. In particular, it behaves well (converges quasi-optimally) in the limit case of the Euler-Bernoulli model. Several numerical results illustrate the performance of our method.

preprint2020arXiv

A short note on plain convergence of adaptive least-squares finite element methods

We show that adaptive least-squares finite element methods driven by the canonical least-squares functional converge under weak conditions on PDE operator, mesh-refinement, and marking strategy. Contrary to prior works, our plain convergence does neither rely on sufficiently fine initial meshes nor on severe restrictions on marking parameters. Finally, we prove that convergence is still valid if a contractive iterative solver is used to obtain the approximate solutions (e.g., the preconditioned conjugate gradient method with optimal preconditioner). The results apply within a fairly abstract framework which covers a variety of model problems.

preprint2016arXiv

A time-stepping DPG scheme for the heat equation

We introduce and analyze a discontinuous Petrov-Galerkin method with optimal test functions for the heat equation. The scheme is based on the backward Euler time stepping and uses an ultra-weak variational formulation at each time step. We prove the stability of the method for the field variables (the original unknown and its gradient weighted by the square root of the time step) and derive a Céa-type error estimate. For low-order approximation spaces this implies certain convergence orders when time steps are not too small in comparison with mesh sizes. Some numerical experiments are reported to support our theoretical results.

preprint2016arXiv

On the DPG method for Signorini problems

We derive and analyze discontinuous Petrov-Galerkin methods with optimal test functions for Signorini-type problems as a prototype of a variational inequality of the first kind. We present different symmetric and non-symmetric formulations where optimal test functions are only used for the PDE part of the problem, not the boundary conditions. For the symmetric case and lowest order approximations, we provide a simple a posteriori error estimate. In a second part, we apply our technique to the singularly perturbed case of reaction dominated diffusion. Numerical results show the performance of our method and, in particular, its robustness in the singularly perturbed case.

preprint2016arXiv

Robust coupling of DPG and BEM for a singularly perturbed transmission problem

We consider a transmission problem consisting of a singularly perturbed reaction diffusion equation on a bounded domain and the Laplacian in the exterior, connected through standard transmission conditions. We establish a DPG scheme coupled with Galerkin boundary elements for its discretization, and prove its robustness for the field variables in so-called balanced norms. Our coupling scheme is the one from [Führer, Heuer, Karkulik: On the coupling of DPG and BEM, arXiv:1508.00630], adapted to the singularly perturbed case by using the scheme from [Heuer, Karkulik: A robust DPG method for singularly perturbed reaction diffusion problems, arXiv:1509.07560]. Essential feature of our method is that optimal test functions have to be computed only locally. We report on various numerical experiments in two dimensions.

preprint2015arXiv

A wirebasket preconditioner for the mortar boundary element method

We present and analyze a preconditioner of the additive Schwarz type for the mortar boundary element method. As a basic splitting, on each subdomain we separate the degrees of freedom related to its boundary from the inner degrees of freedom. The corresponding wirebasket-type space decomposition is stable up to logarithmic terms. For the blocks that correspond to the inner degrees of freedom standard preconditioners for the hypersingular integral operator on open boundaries can be used. For the boundary and interface parts as well as the Lagrangian multiplier space, simple diagonal preconditioners are optimal. Our technique applies to quasi-uniform and non-uniform meshes of shape-regular elements. Numerical experiments on triangular and quadrilateral meshes confirm theoretical bounds for condition and MINRES iteration numbers.

preprint2015arXiv

DPG method with optimal test functions for a fractional advection diffusion equation

We develop an ultra-weak variational formulation of a fractional advection diffusion problem in one space dimension and prove its well-posedness. Based on this formulation, we define a DPG approximation with optimal test functions and show its quasi-optimal convergence. Numerical experiments confirm expected convergence properties, for uniform and adaptively refined meshes.

preprint2015arXiv

On the coupling of DPG and BEM

We develop and analyze strategies to couple the discontinuous Petrov-Galerkin method with optimal test functions to (i) least-squares boundary elements and (ii) various variants of standard Galerkin boundary elements. Essential feature of our methods is that, despite the use of boundary integral equations, optimal test functions have to be computed only locally. We apply our findings to a standard transmission problem in full space and present numerical experiments to validate our theory.

preprint2015arXiv

Optimal additive Schwarz methods for the $hp$-BEM: the hypersingular integral operator in 3D on locally refined meshes

We propose and analyze an overlapping Schwarz preconditioner for the $p$ and $hp$ boundary element method for the hypersingular integral equation in 3D. We consider surface triangulations consisting of triangles. The condition number is bounded uniformly in the mesh size $h$ and the polynomial order $p$. The preconditioner handles adaptively refined meshes and is based on a local multilevel preconditioner for the lowest order space. Numerical experiments on different geometries illustrate its robustness.

preprint2014arXiv

Convergence of adaptive BEM and adaptive FEM-BEM coupling for estimators without h-weighting factor

We analyze adaptive mesh-refining algorithms in the frame of boundary element methods (BEM) and the coupling of finite elements and boundary elements (FEM-BEM). Adaptivity is driven by the two-level error estimator proposed by Ernst P. Stephan, Norbert Heuer, and coworkers in the frame of BEM and FEM-BEM or by the residual error estimator introduced by Birgit Faermann for BEM for weakly-singular integral equations. We prove that in either case the usual adaptive algorithm drives the associated error estimator to zero. Emphasis is put on the fact that the error estimators considered are {not even globally equivalent to weighted-residual error estimators for which recently convergence with quasi-optimal algebraic rates has been derived.

preprint2014arXiv

Multiscale modeling in micromagnetics: existence of solutions and numerical integration

Various applications ranging from spintronic devices, giant magnetoresistance sensors, and magnetic storage devices, include magnetic parts on very different length scales. Since the consideration of the Landau-Lifshitz-Gilbert equation (LLG) constrains the maximum element size to the exchange length within the media, it is numerically not attractive to simulate macroscopic parts with this approach. On the other hand, the magnetostatic Maxwell equations do not constrain the element size, but cannot describe the short-range exchange interaction accurately. A combination of both methods allows to describe magnetic domains within the micromagnetic regime by use of LLG and also considers the macroscopic parts by a non-linear material law using the Maxwell equations. In our work, we prove that under certain assumptions on the non-linear material law, this multiscale version of LLG admits weak solutions. Our proof is constructive in the sense that we provide a linear-implicit numerical integrator for the multiscale model such that the numerically computable finite element solutions admit weak $H^1$-convergence (at least for a subsequence) towards a weak solution.

preprint2013arXiv

Adaptive FEM with optimal convergence rates for a certain class of non-symmetric and possibly non-linear problems

We analyze adaptive mesh-refining algorithms for conforming finite element discretizations of certain non-linear second-order partial differential equations. We allow continuous polynomials of arbitrary, but fixed polynomial order. The adaptivity is driven by the residual error estimator. We prove convergence even with optimal algebraic convergence rates. In particular, our analysis covers general linear second-order elliptic operators. Unlike prior works for linear non-symmetric operators, our analysis avoids the interior node property for the refinement, and the differential operator has to satisfy a Gårding inequality only. If the differential operator is uniformly elliptic, no additional assumption on the initial mesh is posed.

preprint2013arXiv

ZZ-type aposteriori error estimators for adaptive boundary element methods on a curve

In the context of the adaptive finite element method (FEM), ZZ-error estimators named after Zienkiewicz and Zhu are mathematically well-established and widely used in practice. In this work, we propose and analyze ZZ-type error estimators for the adaptive boundary element method (BEM). We consider weakly-singular and hyper-singular integral equations and prove, in particular, convergence of the related adaptive mesh-refining algorithms.

preprint2012arXiv

Classical FEM-BEM coupling methods: nonlinearities, well-posedness, and adaptivity

We consider a (possibly) nonlinear interface problem in 2D and 3D, which is solved by use of various adaptive FEM-BEM coupling strategies, namely the Johnson-Nédélec coupling, the Bielak-MacCamy coupling, and Costabel's symmetric coupling. We provide a framework to prove that the continuous as well as the discrete Galerkin solutions of these coupling methods additionally solve an appropriate operator equation with a Lipschitz continuous and strongly monotone operator. Therefore, the coupling formulations are well-defined, and the Galerkin solutions are quasi-optimal in the sense of a Céa-type lemma. For the respective Galerkin discretizations with lowest-order polynomials, we provide reliable residual-based error estimators. Together with an estimator reduction property, we prove convergence of the adaptive FEM-BEM coupling methods. A key point for the proof of the estimator reduction are novel inverse-type estimates for the involved boundary integral operators which are advertized. Numerical experiments conclude the work and compare performance and effectivity of the three adaptive coupling procedures in the presence of generic singularities.

preprint2012arXiv

Inverse estimates for elliptic boundary integral operators and their application to the adaptive coupling of FEM and BEM

We prove inverse-type estimates for the four classical boundary integral operators associated with the Laplace operator. These estimates are used to show convergence of an h-adaptive algorithm for the coupling of a finite element method with a boundary element method which is driven by a weighted residual error estimator.