Olaf Steinbach

Olaf Steinbach, Barbara Wohlmuth, Linus Wunderlich

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.

Peter Gangl, Mario Gobrial, Olaf Steinbach

In this paper we formulate and analyze a space-time finite element method for the numerical simulation of rotating electric machines where the finite element mesh is fixed in space-time domain. Based on the Babuška--Nečas theory we prove unique solvability both for the continuous variational formulation and for a standard Galerkin finite element discretization in the space-time domain. This approach allows for an adaptive resolution of the solution both in space and time, but it requires the solution of the overall system of algebraic equations. While the use of parallel solution algorithms seems to be mandatory, this also allows for a parallelization simultaneously in space and time. This approach is used for the eddy current approximation of the Maxwell equations which results in an elliptic-parabolic interface problem. Numerical results for linear and nonlinear constitutive material relations confirm the applicability, efficiency and accuracy of the proposed approach.

Stefan Dohr, Michal Merta, Günther Of et al.

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.

Stefan Dohr, Olaf Steinbach

In this note we describe a space-time boundary element discretization of the heat equation and an efficient and robust preconditioning strategy which is based on the use of boundary integral operators of opposite orders, but which requires a suitable stability condition for the boundary element spaces used for the discretization. We demonstrate the method for the simple spatially one-dimensional case. However, the presented results, particularly the stability analysis of the boundary element spaces, can be used to extend the method to the two- and three-dimensional problem.

Michael Reichelt, Michael Wiesheu, Melina Merkel et al.

Space-time methods promise more efficient time-domain simulations, in particular of electrical machines. However, most approaches require the motion to be known in advance so that it can be included in the space-time mesh. To overcome this problem, this paper proposes to use the well-known air-gap element for the rotor-stator coupling of an isogeometric machine model. First, we derive the solution in the air-gap region and then employ it to couple the rotor and stator. This coupling is angle dependent and we show how to efficiently update the coupling matrices to a different angle, avoiding expensive quadrature. Finally, the resulting time-dependent problem is solved in a space-time setting. The spatial discretization using isogeometric analysis is particularly suitable for coupling via the air-gap element, as NURBS can exactly represent the geometry of the air-gap. Furthermore, the model including the air-gap element can be seamlessly transferred to the space-time setting. However, the air-gap element is well known in the literature. The originality of this work is the application to isogeometric analysis and space-time.

Christoph M. Augustin, Gerhard A. Holzapfel, Olaf Steinbach

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.