BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Olaf Steinbach

Olaf Steinbach appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
math.NANumerical Analysismath.OCmath.APBiological Physicsphysics.med-ph

Catalog footprint

What is connected

11works
6topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

11 published item(s)

preprint2022arXiv

Robust finite element discretization and solvers for distributed elliptic optimal control problems

We consider standard tracking-type, distributed elliptic optimal control problems with $L^2$ regularization, and their finite element discretization. We are investigating the $L^2$ error between the finite element approximation $u_{\varrho h}$ of the state $u_\varrho$ and the desired state (target) $\bar{u}$ in terms of the regularization parameter $\varrho$ and the mesh size $h$ that leads to the optimal choice $\varrho = h^4$. It turns out that, for this choice of the regularization parameter, we can devise simple Jacobi-like preconditioned MINRES or Bramble-Pasciak CG methods that allow us to solve the reduced discrete optimality system in asymptotically optimal complexity with respect to the arithmetical operations and memory demand. The theoretical results are confirmed by several benchmark problems with targets of various regularities including discontinuous targets.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

A generalized inf-sup stable variational formulation for the wave equation

In this paper, we consider a variational formulation for the Dirichlet problem of the wave equation with zero boundary and initial conditions, where we use integration by parts in space and time. To prove unique solvability in a subspace of $H^1(Q$) with $Q$ being the space-time domain, the classical assumption is to consider the right-hand side $f$ in $L^2(Q)$. Here, we analyze a generalized setting of this variational formulation, which allows us to prove unique solvability also for $f$ being in the dual space of the test space, i.e., the solution operator is an isomorphism between the ansatz space and the dual of the test space. This new approach is based on a suitable extension of the ansatz space to include the information of the differential operator of the wave equation at the initial time $t=0$. These results are of utmost importance for the formulation and numerical analysis of unconditionally stable space-time finite element methods, and for the numerical analysis of boundary element methods to overcome the well-known norm gap in the analysis of boundary integral operators.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Numerical results for an unconditionally stable space-time finite element method for the wave equation

In this work, we introduce a new space-time variational formulation of the second-order wave equation, where integration by parts is also applied with respect to the time variable, and a modified Hilbert transformation is used. For this resulting variational setting, ansatz and test spaces are equal. Thus, conforming finite element discretizations lead to Galerkin--Bubnov schemes. We consider a conforming tensor-product approach with piecewise polynomial, continuous basis functions, which results in an unconditionally stable method, i.e., no CFL condition is required. We give numerical examples for a one- and a two-dimensional spatial domain, where the unconditional stability and optimal convergence rates in space-time norms are illustrated.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Robust discretization and solvers for elliptic optimal control problems with energy regularization

We consider the finite element discretization and the iterative solution of singularly perturbed elliptic reaction-diffusion equations in three-dimensional computational domains. These equations arise from the optimality conditions for elliptic distributed optimal control problems with energy regularization that were recently studied by M.~Neumüller and O.~Steinbach (2020). We provide quasi-optimal a priori finite element error estimates which depend both on the mesh size $h$ and on the regularization parameter $\varrho$. The choice $\varrho = h^2$ ensures optimal convergence which only depends on the regularity of the target function. For the iterative solution, we employ an algebraic multigrid preconditioner and a balancing domain decomposition by constraints (BDDC) preconditioner. We numerically study robustness and efficiency of the proposed algebraic preconditioners with respect to the mesh size $h$, the regularization parameter $\varrho$, and the number of subdomains (cores) $p$. Furthermore, we investigate the parallel performance of the BDDC preconditioned conjugate gradient solver.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Space-time finite element discretization of parabolic optimal control problems with energy regularization

We analyze space-time finite element methods for the numerical solution of distributed parabolic optimal control problems with energy regularization in the Bochner space $L^2(0,T;H^{-1}(Ω))$. By duality, the related norm can be evaluated by means of the solution of an elliptic quasi-stationary boundary value problem. When eliminating the control, we end up with the reduced optimality system that is nothing but the variational formulation of the coupled forward-backward primal and adjoint equations. Using Babuška's theorem, we prove unique solvability in the continuous case. Furthermore, we establish the discrete inf-sup condition for any conforming space-time finite element discretization yielding quasi-optimal discretization error estimates. Various numerical examples confirm the theoretical findings. We emphasize that the energy regularization results in a more localized control with sharper contours for discontinuous target functions, which is demonstrated by a comparison with an $L^2$ regularization and with a sparse optimal control approach.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Unstructured space-time finite element methods for optimal control of parabolic equations

This work presents and analyzes space-time finite element methods on fully unstructured simplicial space-time meshes for the numerical solution of parabolic optimal control problems. Using Babuška's theorem, we show well-posedness of the first-order optimality systems for a typical model problem with linear state equations, but without control constraints. This is done for both continuous and discrete levels. Based on these results, we derive discretization error estimates. Then we consider a semilinear parabolic optimal control problem arising from the Schlögl model. The associated nonlinear optimality system is solved by Newton's method, where a linear system, that is similar to the first-order optimality systems considered for the linear model problems, has to be solved at each Newton step. We present various numerical experiments including results for adaptive space-time finite element discretizations based on residual-type error indicators. In the last two examples, we also consider semilinear parabolic optimal control problems with box constraints imposed on the control.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Unstructured Space-Time Finite Element Methods for Optimal Sparse Control of Parabolic Equations

We consider a space-time finite element method on fully unstructured simplicial meshes for optimal sparse control of semilinear parabolic equations. The objective is a combination of a standard quadratic tracking-type functional including a Tikhonov regularization term and of the $L^1$-norm of the control that accounts for its spatio-temporal sparsity. We use a space-time Petrov-Galerkin finite element discretization for the first-order necessary optimality system of the associated discrete optimal sparse control problem. The discretization is based on a variational formulation that employs piecewise linear finite elements simultaneously in space and time. Finally, the discrete nonlinear optimality system that consists of coupled forward-backward state and adjoint state equations is solved by a semismooth Newton method.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

A parallel solver for a preconditioned space-time boundary element method for the heat equation

We describe a parallel solver for the discretized weakly singular space-time boundary integral equation of the spatially two-dimensional heat equation. The global space-time nature of the system matrices leads to improved parallel scalability in distributed memory systems in contrast to time-stepping methods where the parallelization is usually limited to spatial dimensions. We present a parallelization technique which is based on a decomposition of the input mesh into submeshes and a distribution of the corresponding blocks of the system matrices among processors. To ensure load balancing, the distribution is based on a cylic decomposition of complete graphs. In addition, the solution of the global linear system requires the use of an efficient preconditioner. We present a robust preconditioning strategy which is based on boundary integral operators of opposite order, and extend the introduced parallel solver to the preconditioned system.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Classical and all-floating FETI methods for the simulation of arterial tissues

High-resolution and anatomically realistic computer models of biological soft tissues play a significant role in the understanding of the function of cardiovascular components in health and disease. However, the computational effort to handle fine grids to resolve the geometries as well as sophisticated tissue models is very challenging. One possibility to derive a strongly scalable parallel solution algorithm is to consider finite element tearing and interconnecting (FETI) methods. In this study we propose and investigate the application of FETI methods to simulate the elastic behavior of biological soft tissues. As one particular example we choose the artery which is - as most other biological tissues - characterized by anisotropic and nonlinear material properties. We compare two specific approaches of FETI methods, classical and all-floating, and investigate the numerical behavior of different preconditioning techniques. In comparison to classical FETI, the all-floating approach has not only advantages concerning the implementation but in many cases also concerning the convergence of the global iterative solution method. This behavior is illustrated with numerical examples. We present results of linear elastic simulations to show convergence rates, as expected from the theory, and results from the more sophisticated nonlinear case where we apply a well-known anisotropic model to the realistic geometry of an artery. Although the FETI methods have a great applicability on artery simulations we will also discuss some limitations concerning the dependence on material parameters.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Trace and flux a priori error estimates in finite element approximations of Signorni-type problems

Variational inequalities play in many applications an important role and are an active research area. Optimal a priori error estimates in the natural energy norm do exist but only very few results in other norms exist. Here we consider as prototype a simple Signorini problem and provide new optimal order a priori error estimates for the trace and the flux on the Signorini boundary. The a priori analysis is based on the exact and a mesh-dependent Steklov-Poincaré operator as well as on duality in Aubin-Nitsche type arguments. Numerical results illustrate the convergence rates of the finite element approach.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Boundary Integral Equations for the Laplace-Beltrami Operator

We present a boundary integral method, and an accompanying boundary element discretization, for solving boundary-value problems for the Laplace-Beltrami operator on the surface of the unit sphere $§$ in $\mathbb{R}^3$. We consider a closed curve ${\cal C}$ on ${\cal S}$ which divides ${\cal S}$ into two parts ${\cal S}_1$ and ${\cal S}_2$. In particular, ${\cal C} = \partial {\cal S}_1$ is the boundary curve of ${\cal S}_1$. We are interested in solving a boundary value problem for the Laplace-Beltrami operator in $§_2$, with boundary data prescribed on $\C$.

0 citations0 reviewsNo code linkedOpen access